Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB143 44,301 (T)9→8 intergenic (+132/‑359) oprD_1 → / → czcR_2 Porin D/Transcriptional activator protein CzcR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB14344,2930T.98.5% 302.5 / ‑3.2 65intergenic (+124/‑367)oprD_1/czcR_2Porin D/Transcriptional activator protein CzcR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (29/35);  minor base C (1/0);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTCTGTAATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGGAACGAAGCCA  >  ACB143/44160‑44422
                                                                                                                                     |                                                                                                                                 
gATTCTGTAATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCttt                                                                                                                                 >  7:251294/1‑136 (MQ=255)
                 ccagtgCCAGTTATCCATTGTGGGGTCCTGTTAGTAAAACGATAATTAAGCATTCAGCTTGAGTTGCTTAATTGCTCCTTTGGTAAAAGTGGGATATTTGGCGCCCATGGGGCGCCTTTTTTTTGCCTAAAATTc                                                                                                                <  8:122406/130‑1 (MQ=255)
                     ggCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAg                                                                                                             <  5:353768/134‑1 (MQ=255)
                      gCCATTTATCCATTGGCGTGTCCTGTTAGTAAAACTAGCATTAGGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGTATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAg                                                                                                          <  4:215035/136‑1 (MQ=255)
                      gCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAg                                                                                                          <  3:289786/136‑1 (MQ=255)
                      gCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAg                                                                                                          <  6:277801/136‑1 (MQ=255)
                        cAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                                        <  4:265819/136‑1 (MQ=255)
                        cAGTTATACATTGGCGTGTGCTGTTAGTAAAACGAGGATTACTCATTCACCTTGAGTGGCTTTAGTGCTCCTTTGGTAAAAGGTGGTTATTTGGCGCCCTTTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                                        <  8:195043/136‑1 (MQ=255)
                           ttATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGtt                                                                                                     <  4:306143/136‑1 (MQ=255)
                                cATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGAt                                                                                                 <  8:79955/135‑1 (MQ=255)
                                  ttGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGt                                                                                              <  3:110116/136‑1 (MQ=255)
                                  ttGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGACTAAAATTCAAGTAGGCGTTCGATTg                                                                                               >  5:160203/1‑135 (MQ=255)
                                       gggTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgca                                                                                         <  3:52821/136‑1 (MQ=255)
                                        ggTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTAAAGTAGGCGTTCGATTGTAcccat                                                                                        >  6:6099/1‑136 (MQ=255)
                                         gTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatc                                                                                       <  1:100200/136‑1 (MQ=255)
                                         gTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatc                                                                                       <  2:101471/136‑1 (MQ=255)
                                           ccTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcg                                                                                      >  6:256747/1‑135 (MQ=255)
                                           ccTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcg                                                                                      >  3:256485/1‑135 (MQ=255)
                                            cTGTTAGTAAAACGAGAATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGTCTCCCATTGGGCGACTTTTTTTTGCCTAATATTCAAGTAGGCGTTCGATTGTAcgcatcgca                                                                                    >  6:90776/1‑136 (MQ=255)
                                             tGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcgcat                                                                                   >  4:99103/1‑136 (MQ=255)
                                                  gTAAAACGAGCATTACGCATTCACCTTGAGTTGCTTAAGTTCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                               <  3:363905/135‑1 (MQ=255)
                                                  gTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                               <  5:194862/135‑1 (MQ=255)
                                                     aaaCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTg                                                                            <  4:143998/135‑1 (MQ=255)
                                                      aaCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGaa                                                                          <  1:124971/136‑1 (MQ=255)
                                                      aaCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGaa                                                                          <  1:43588/136‑1 (MQ=255)
                                                         gAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGtt                                                                       <  1:271248/136‑1 (MQ=255)
                                                         gAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGtt                                                                       <  8:155760/136‑1 (MQ=255)
                                                            cATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCAt                                                                    <  8:52029/136‑1 (MQ=255)
                                                              ttACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGa                                                                  <  7:46023/136‑1 (MQ=255)
                                                              ttACGCATTCAACTTGAGTGGCTTCAGTGCTCCTTTGGTACAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGa                                                                  <  2:303581/136‑1 (MQ=255)
                                                                  gCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTac                                                              <  8:218348/136‑1 (MQ=255)
                                                                  gCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTac                                                              <  3:203015/136‑1 (MQ=255)
                                                                  gCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCACGCCGGTGAAGTTCATGACTac                                                              <  1:239674/136‑1 (MQ=255)
                                                                       cACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCaa                                                         >  2:127153/1‑136 (MQ=255)
                                                                        aCCTTGAGTGGCTTTAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCTCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCaa                                                         <  2:203935/135‑1 (MQ=255)
                                                                        aCCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCaa                                                         <  5:60711/135‑1 (MQ=255)
                                                                        aCCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTAATGACTACACCAAt                                                        >  5:285901/1‑136 (MQ=255)
                                                                         ccTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAAt                                                        >  7:72701/1‑135 (MQ=255)
                                                                             gAGTTGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCTCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAcc                                                    <  6:188099/135‑1 (MQ=255)
                                                                             gAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAcc                                                    <  5:230108/135‑1 (MQ=255)
                                                                             gAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                                     <  5:45467/134‑1 (MQ=255)
                                                                             gAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                                     <  6:255310/134‑1 (MQ=255)
                                                                               gTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGt                                                 >  2:272239/1‑136 (MQ=255)
                                                                                 ggCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCTCCCATTGGGCGCCTTTTTTTTGCCTAAAATTAAAGTAGGCGTTCGATTGTACGCATCGCATGCCCGGGAATTTCAAGACTACAACAATGACCACtc                                                >  4:284145/1‑135 (MQ=255)
                                                                                 ggCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGGCGAAAATTCAAGTAGGCGTTCGATTGTACGCATAGCATGCCGGTGAAGTTCATGACTTCACCAAGGACCAGtcg                                               >  4:311425/1‑135 (MQ=255)
                                                                                 ggCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGtct                                               >  1:280187/1‑136 (MQ=255)
                                                                                 ggCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGtc                                                >  8:228712/1‑135 (MQ=255)
                                                                                   cTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTctct                                             >  7:226832/1‑136 (MQ=255)
                                                                                   cTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTctct                                             >  7:231573/1‑136 (MQ=255)
                                                                                     tAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCCATTGTACGAATCGCAAGCCGGTGAAGTTCATTACTAAACCAATGACCAGtctcat                                            >  8:254237/1‑133 (MQ=255)
                                                                                     tAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTAAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCtt                                            >  2:95671/1‑135 (MQ=255)
                                                                                     tAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTCGGCCGCTTTTTTTTGGCTAAATTTCAAGTAGTCGTTGGATTGTACCGATAACATGCCGGTGAAGTTCATGAATTAACACAAGCCCAGTCTCtt                                            >  8:343123/1‑135 (MQ=255)
                                                                                       aGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGTCTTTCGAATGTACGCAAAGCATGCCGGTGAAGTTCATGTATACACCAATGACCAGTCTCATCaa                                         >  7:331476/1‑136 (MQ=255)
                                                                                           cTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCTCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAAtgt                                      <  4:106978/135‑1 (MQ=255)
                                                                                                ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGa                                >  3:361318/1‑136 (MQ=255)
                                                                                                ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTAAATGTGAGCGa                                >  8:326459/1‑136 (MQ=255)
                                                                                                ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTAAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGTCTACACCAATGACCAGTCTATTCAATGTGAGCGa                                >  1:271619/1‑136 (MQ=255)
                                                                                                 tGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTATGCGTTCGATTGTACGCATCGCATGCCGGTGAATTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGAc                               >  4:6192/1‑136 (MQ=255)
                                                                                                              atatTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTAGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCt                  <  1:215991/136‑1 (MQ=255)
                                                                                                               tatTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTg                 <  8:120084/136‑1 (MQ=255)
                                                                                                                atTTGGCGCCCATTGGGCGCCTTTTTTTTGGCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGTCTACACCAATGAGCAGTCTCTTCAATGTGAGCGCCCGCGCCGTcgactgc                >  8:345431/1‑130 (MQ=255)
                                                                                                                atTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCt                  >  5:125710/1‑134 (MQ=255)
                                                                                                                       gcCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGCCCTCGCCGTCATCTGAGCGGAAc         >  8:329824/1‑136 (MQ=255)
                                                                                                                       gcCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGAATAAACCAATGACAAGTCTCTTAAATGTGAGCGCCCTCGCCGTCATCTGAGCGGAAc         >  7:138970/1‑136 (MQ=255)
                                                                                                                        cccATTGGGCGCCTTTTTTTTGGCTAAAATTAAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGGAAc         >  5:34224/1‑135 (MQ=255)
                                                                                                                              gggCGCCTTTTTTTTGGCTAAAATTCAAGTAGGCGTTGGATTGTACGGATCGCATGCCGGTGAAGTTGATGACTACCCCAATGACCAGTGTCTTCAATGTGAGCGACCTCGCCGTGATCTGCGCGGAACGAAggcc  >  4:6933/1‑133 (MQ=255)
                                                                                                                                     |                                                                                                                                 
GATTCTGTAATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGGAACGAAGCCA  >  ACB143/44160‑44422

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: