Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913144,1840TG37.8% 145.2 / 24.0 92G161G (GGT→GGGyadHputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (13/43);  new base G (34/0);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 3.27e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTGGCGCGTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGTTTATATGATC  >  NC_000913/144038‑144334
                                                                                                                                                   |                                                                                                                                                       
ggcgttgcgCGTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTTTTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGgtgtgt                                                                                                                                                     <  2:593897/151‑1 (MQ=255)
ggcgtggcgCGTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGgtgtgt                                                                                                                                                     <  1:1053969/151‑1 (MQ=255)
  cgtggcgCGTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGttt                                                                                                                                                   <  1:502172/151‑1 (MQ=255)
   gtgccgCGTGGTCAGTTTGTTGTAATTCTGGTGAAGGCAATTTAACTGTTTTTTGTGCCATTTCAGGTGCATTAGAGGGTATTCGCTGCCTTAAACCTGGTGTTCACGGCGGCGGTGTTTCCCCATCCCGGGCTGCTGAACGttgactttt                                                                                                                                                  >  2:1039640/1‑145 (MQ=255)
      gcgcGTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCa                                                                                                                                              <  2:828207/151‑1 (MQ=255)
      gcgcGTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCa                                                                                                                                              <  1:805325/151‑1 (MQ=255)
       cgggTGGTCTGTTTGTTGGCAGTCTGGTGACGGCAATTTCCCTGTTTTTTGTGCCATTTCAGGTTCATTCGTGGGTATTCGTTGCCTTACCGCTGGTGCTCACGGCCGTGTTGTTCTCCCTTGCGGGTCTGCTGAACGGTGTGTTTGCCaa                                                                                                                                             <  2:84195/148‑1 (MQ=255)
        gcgTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCaaa                                                                                                                                            <  1:588285/151‑1 (MQ=255)
             gTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGtt                                                                                                                                       <  1:708208/151‑1 (MQ=255)
             gTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGtt                                                                                                                                       <  2:948337/151‑1 (MQ=255)
             gTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGtt                                                                                                                                       <  1:719323/151‑1 (MQ=255)
              tCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGAGTTGTTGTCCCTCGCGGGGTTGCTGGACGGGGCGTTTGCCgaaatcgtc                                                                                                                                      >  1:225945/1‑146 (MQ=255)
                   ttgttgGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTTCGGGGTTGGTGAAAGGGGTGTTTGGCAAAAAGTTCGAt                                                                                                                                   >  1:537853/1‑149 (MQ=255)
                    tgttgGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGGGTTGTTTTCCCTTGCGGGGTTTGTGGACGGGGGGTTTGTCCAAAAGtt                                                                                                                                       >  2:6946/1‑144 (MQ=255)
                       tgGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTTCCTTTCCGGGTTTTTGGAAGGGGTGTTTGTCCAAAAGGTCGATGACa                                                                                                                               >  1:380088/1‑149 (MQ=255)
                       tgGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCTCCTGCGGGGTTGCTGGACGGGGTGGTTTGCAAAACGGTCGATGAAATc                                                                                                                             >  2:1418942/1‑151 (MQ=255)
                        ggCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAAcaaatcggaagagc                                                                                                                             <  1:444055/150‑14 (MQ=255)
                            tgctagTGACGGCATTTTCACTGTTGTTTCTTCTCTATCGGGGTCATTCGTGGGTATTCGTTGGTTTAACGGTTGTGCTCACGGCGGTGTTGTTCTCGCTTGCAGGTTTGCTGAACGGTGTGTTTTCTAAAACGTTCGATGa                                                                                                                                 <  2:1131469/137‑1 (MQ=255)
                             tCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCt                                                                                                                        <  1:903651/150‑1 (MQ=255)
                              cTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAACGGGGGGTTTGCCAAAACGTTCGATGACATCAGCCTgg                                                                                                                      >  1:646977/1‑151 (MQ=255)
                               tGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGGTGTTCTCCCCTTCGGGGTTTCTGAACGGCGGGGTTGCCAAAACGGTCGATGACaaccccc                                                                                                                         >  1:1054487/1‑141 (MQ=255)
                                ggTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGt                                                                                                                     <  1:1293557/150‑1 (MQ=255)
                                ggTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTg                                                                                                                    <  2:1143130/151‑1 (MQ=255)
                                    aCGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCaa                                                                                                                <  1:1195659/151‑1 (MQ=255)
                                    aCGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCCTTCGGGGTTTCTGAACGGTGTGTTTGTCAAAACGTTCGATGACATCAGCCTGGTGcc                                                                                                                  >  1:3823/1‑149 (MQ=255)
                                     cGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGGTTTTTTCCCCTTCGGGGTTTGTGAACGGGGTGTTTGCCAAAAAGTTCGATAAAATCAGGCCGGTGccca                                                                                                                >  1:306488/1‑148 (MQ=255)
                                        cAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTTCGGGGTTGCTGAACGGGGGGTTTGCCAAAACGGTCGATTACATCAACCCGGTGCCAAcact                                                                                                            >  2:155935/1‑148 (MQ=255)
                                        cAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTGCGGGGTTGCCGAACGGGGGGGTTGCCCAAACGTTCGATGACATCAGCCTGGGGCCAACCtt                                                                                                            >  1:670220/1‑151 (MQ=255)
                                        cAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGCTGCCAACCtt                                                                                                            <  2:407962/151‑1 (MQ=255)
                                        cAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGCGTTATTCTCCCTTGCGGGGTTGCTGAACGGGGGGGTTTGCCAAAAGTTCGATGACAACAGACCGGTGCCAACCtt                                                                                                            >  1:954537/1‑151 (MQ=255)
                                         aaTTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTGCGGGGTTGGTGAACGGGGTGTTTTCCCAAAAGTTCGATGACATCAGCCTGGTGCCCACCttt                                                                                                           >  1:504190/1‑151 (MQ=255)
                                         aaTTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCCTTCGGGTTTGCTGAACGGGGTGTTTTCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCttt                                                                                                           >  1:793864/1‑151 (MQ=255)
                                          aTTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTtg                                                                                                          <  1:1049856/151‑1 (MQ=255)
                                              cACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTTCCTTTCGGGGTTGCTGAAAGGGGTGTTTGCCAAAACGTTCGATGACATCAGCCTTGTGCCAACCTTTTTGtt                                                                                                      >  1:814248/1‑151 (MQ=255)
                                               aCTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTTCCTTTCGGGTTTTCTGAAAGGGGTGTTTTCCCAAACGTTTGGTGACATCAGCCCGGTGCCCACCttt                                                                                                           >  2:1426431/1‑145 (MQ=255)
                                                  gTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTaa                                                                                                    <  2:1099483/149‑1 (MQ=255)
                                                  gTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAAc                                                                                                   <  1:1419944/150‑1 (MQ=255)
                                                   ttttttGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGGTGTTTTCCCTTGCGGGGTTGGTGGACCGGGTGTTTGCCCAAACGTTTGATGACATCAGGCCGGTGGCAACCttt                                                                                                           >  1:1017211/1‑141 (MQ=255)
                                                    tatttGTGCCATTTCAGTTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCGTGCGGGTTTGCGGAACGGTGTTTTTGCAAAAACGTTCGATGACATCAGCCTGGTGCCACCCTTTGTGTTaa                                                                                                    <  2:599391/145‑1 (MQ=255)
                                                      tttGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTTTTTTCCCTTGCGGGGTTTCTGAACGGTGTGTTTTCCAAAACGTTTGGTTACATTAGCCTGGGGCCAACCTTTTTGTTTAcaccac                                                                                              >  2:922725/1‑151 (MQ=255)
                                                      tttGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTCTGGGGTTGCTGAAAGGGGTGTTTTCCCAAACGTTCGATGGAATTAGCCCAGTGGCAACCTTTGTGTTaaaccc                                                                                                >  2:1062551/1‑145 (MQ=255)
                                                        tgtgCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTGCGGGTTTGCTGAAAGGGGGGTTTGCCAAAAAGTTCGATGAAATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTc                                                                                            >  1:1136314/1‑151 (MQ=255)
                                                         gtgCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAACGGGGTGTTTTCCCAAACGTTCGATGACATCAGACCGGTGCCCACCTTTGTGTTTACGCCACTc                                                                                            >  1:1250300/1‑150 (MQ=255)
                                                          tgCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGGTTTTCTTCCTTGCGGGTTTGGTGGACGGGGTGTTTGGCCCAAACGTCGATGACATCCGCCTGGGGCCAACCTTTGTGGTTACGCCACTc                                                                                            >  2:1444985/1‑149 (MQ=255)
                                                             cATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCCTGCGGGGTTTGTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCCACCTTTGTGTTAACGCCCCTCAc                                                                                          >  1:158667/1‑148 (MQ=255)
                                                              aTTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAAAGGGGGGTTTGCCCAAACGTTCGATGACATCAGCCTGGGGCCAACCTTTGTGTTTACGCCACTCACGTAt                                                                                      >  1:1566268/1‑151 (MQ=255)
                                                               tttCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTTCGGGGTTGCTGAAAGGGGTGTTTGCCAAAAAGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTAtt                                                                                     >  1:963283/1‑151 (MQ=255)
                                                               tttCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTAt                                                                                      <  2:254759/150‑1 (MQ=255)
                                                                ttcgcgTGCATTAGCGGATATTCGTTGCCTTAACGCTGGGGATAACTGCGGTGTTATTCTGCCTTGCGGGTGGGCTGAACGGAGGGTTTGCCAAAACGCTAGATGACAACAGCCTGGTGCCAACCTTTGCGTCAACGTAACTCGCGTAtt                                                                                     <  2:1151451/145‑1 (MQ=255)
                                                                ttCAGGTGCATTCGTGGGTATTCGTTTCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCGCTCACGTa                                                                                       <  2:907553/148‑1 (MQ=255)
                                                                 tCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGGTCTTCCCTTCGGGGTTGCTGAACCGGGTGTTTGCCCAAACGGTCGCTGACATTAGGCTGGTGGCAACCCTTGTGTTTACGGCACTCACCGAttt                                                                                    >  1:560385/1‑150 (MQ=255)
                                                                  cAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACTGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATGAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTg                                                                                   <  1:1063090/150‑1 (MQ=255)
                                                                  cAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTTTCCCTTTCGGGTTTGCTGAAAGGGGTGTTTGCCCAAAAGTTCGATGACATCAGCCTGGTGCCCACCTTTGTGTTTACGCCACTCACGTAttttg                                                                                  >  1:255247/1‑149 (MQ=255)
                                                                  cAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAACGGGGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGGGCCAACCTTTGTGTTAACGCCACTCACGTAttt                                                                                    >  1:566950/1‑149 (MQ=255)
                                                                   aGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTggg                                                                                 <  2:1142492/151‑1 (MQ=255)
                                                                    ggTGCATTCGTGGGTATTCGTTGCCGTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGTTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGc                                                                                <  2:90080/151‑1 (MQ=255)
                                                                         aTTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTTCGGGGTTGCTGAACGGTGTGTTTTCCCAAACGGTCGATGACATCAGCCTGGTGGCCACCTTTGTGTTTAAGCCACTCACGTATTTTGGCggggg                                                                           >  1:38500/1‑150 (MQ=255)
                                                                           tCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTTCCTTGCGGGGTTGCTGGACGGGGTGTTTGCCAAAACGTTCGATGACCTCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCTCGTATTTGGGCGGGGTc                                                                          >  2:30287/1‑150 (MQ=255)
                                                                           tCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTc                                                                          <  2:276087/150‑1 (MQ=255)
                                                                             gTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTcgt                                                                        <  2:751309/150‑3 (MQ=255)
                                                                                ggTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTAc                                                                    <  1:1295126/151‑1 (MQ=255)
                                                                                    ttCGTTGCCTTAACGCTGGTGCTCTCGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTTGGCGGGGTCTTTTACTCa                                                                 <  2:153971/150‑1 (MQ=255)
                                                                                      cGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTTCCTTGCGGGTTTGTTGGACGGTGTGTTTGCCAAAACGTTCGGTGAACTCAGCCTGGTGCCAACCTTTTTGGTTACGCCCCTCCCGTATTTGGGCGGGGGCTTTTTCTCAc                                                                >  2:188821/1‑149 (MQ=255)
                                                                                        ttGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGa                                                             <  2:634031/150‑1 (MQ=255)
                                                                                          gCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACt                                                           <  1:802426/150‑1 (MQ=255)
                                                                                          gCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACt                                                           <  1:1537744/150‑1 (MQ=255)
                                                                                           ccTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTTCGGGGTTGCTGAACGGGGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTTGGCGGGGGCTTTTACTCACTGACttt                                                         >  1:1573090/1‑151 (MQ=255)
                                                                                           ccTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTTCTGAACGGGGTGTTTGCCAAAAAGTTCGACGGAATCAGGCTGGTGCCAACCCTTGTGTTAACGCCCCTCCCGTATTTGGGCGGGGGCTTTTAATCACTGACttt                                                         >  2:867714/1‑151 (MQ=255)
                                                                                           ccTTAACGCTGGTGCTCACGGCGGGGTTGTTCTCTCCTGCGGGGTTGCTGAACGGGGTGTTTGCCCAAACGTTCTATGACATCCGCCTGGTGCCCACCTTTGTGTTAAAGCCACCCCCGTATTTGGGGGGGGGGTTTTTCTCACTTACttt                                                         >  2:618256/1‑151 (MQ=255)
                                                                                             ttAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAAAGGGGTGTTTGCCAAAAAGTTCGATGACATCAGCCTGGGGCCCACCTTTGTGTTAACGCCACTCACGTATTTTGGGGGGGGCTTTTTCTCACTGACTttgt                                                       >  1:663195/1‑151 (MQ=255)
                                                                                              tAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTTCGGGGTTGCTGAACGGGGTGTTTTCCAAAAAGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAAAGCCACTCACGGATTTGGGCGGGGGCTTTTACTCACTGACtttttt                                                      >  1:476213/1‑151 (MQ=255)
                                                                                               aaCGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTttgtt                                                      <  1:101587/150‑1 (MQ=255)
                                                                                                aCGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGGTTGCTGAACGGGGTGTTTGCCAAAACGTTCGGTGACATCAGCCTGGTGCCAACCTTTTTGTTAACGCCACTCACGTATTTTGGCGGGGGCTTTTTCTCACTGACTttgtt                                                      >  1:1094427/1‑149 (MQ=255)
                                                                                                 cGCTGGTGCTCACGGCGGTGTTTTTCTCCCTTGCGGGGTTGCCGAAAAGGGTGGTTGCCCAAACGTTTGATGACATCAGCCCGGTGCCAACCTTTTTGTTAACGCCCATCACCTATTTGGGGGGGGGCTTTTACTCACTGACTTTTTTtcc                                                   >  2:637382/1‑151 (MQ=255)
                                                                                                 cGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTgcc                                                   <  1:237341/151‑1 (MQ=255)
                                                                                                       tGCTCACGGCGGCGTTGTTCCCCCTTCCGGGTTTGCTGAACGGGGCGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTTTGTTTACGCCACTCACGTATTTGGGCGGGGCCTTTTAATCACTGACTTTGTTGCCGCCGtt                                             >  2:1223198/1‑151 (MQ=255)
                                                                                                           cACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTgg                                         <  1:393043/151‑1 (MQ=255)
                                                                                                            aCGGCGGCGTTATTCTCCCTTGCGGGGTTGCTGTACCGGGTTTTTGCCCAAACGTTCGCTGACATTCGCCCGGTGGCAACCTTTTTGTTAAAGCCCCTCACGTATTTGGGCGGGGGCTTTTTCTCACTGACTTTGTTgccgc                                                 >  2:109963/1‑142 (MQ=255)
                                                                                                             cggcggTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGGCATCAGCCTGGTGCCACCCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGGCTTTGTTGCCGCCGGTCTGGCa                                       <  2:1266507/151‑1 (MQ=255)
                                                                                                              ggcggTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCaa                                      <  1:204361/151‑1 (MQ=255)
                                                                                                               gcggTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAg                                     <  2:432422/151‑1 (MQ=255)
                                                                                                               gcggTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAg                                     <  1:289058/151‑1 (MQ=255)
                                                                                                                cggTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAgg                                    >  1:394719/1‑151 (MQ=255)
                                                                                                                    gttgttCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTg                                <  1:427307/151‑1 (MQ=255)
                                                                                                                              cTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCAACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTg                       >  2:901527/1‑150 (MQ=255)
                                                                                                                               ttGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGa                      >  1:1392140/1‑150 (MQ=255)
                                                                                                                                 gCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAAcc                   >  2:1396756/1‑151 (MQ=255)
                                                                                                                                  cGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAAccc                  >  1:1526863/1‑151 (MQ=255)
                                                                                                                                    ggTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCaa                <  1:1531910/151‑1 (MQ=255)
                                                                                                                                     gTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAAt               >  2:324614/1‑151 (MQ=255)
                                                                                                                                         gCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGtt           <  1:1612635/151‑1 (MQ=255)
                                                                                                                                            gAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGTTtat        >  1:62313/1‑151 (MQ=255)
                                                                                                                                              aCGGTGTGTTTGCCAAAGCGTTCGATTACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGTTTatat      <  2:1515705/151‑1 (MQ=255)
                                                                                                                                               cGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGTTTATATg     >  1:152856/1‑151 (MQ=255)
                                                                                                                                               cGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGTTTATATg     <  1:416203/151‑1 (MQ=255)
                                                                                                                                                  gtgtgtTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGTTTATATGAt   >  1:228693/1‑151 (MQ=255)
                                                                                                                                                   tgtgtTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGTTTATATGATc  >  2:1177736/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
GGCGTGGCGCGTGGTCTGTTTGTTGGCATTCTGGTGACGGCAATTTCACTGTTTTTTGTGCCATTTCAGGTGCATTCGTGGGTATTCGTTGCCTTAACGCTGGTGCTCACGGCGGTGTTGTTCTCCCTTGCGGGTTTGCTGAACGGTGTGTTTGCCAAAACGTTCGATGACATCAGCCTGGTGCCAACCTTTGTGTTAACGCCACTCACGTATTTGGGCGGGGTCTTTTACTCACTGACTTTGTTGCCGCCGTTCTGGCAAGGGCTGTCGCACCTGAACCCAATCGTTTATATGATC  >  NC_000913/144038‑144334

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