Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913385,7970TG29.8% 116.4 / 13.1 69V189G (GTG→GGG) tauAtaurine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (10/37);  new base G (20/0);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 5.18e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTGAAGGGGAAT  >  NC_000913/385558‑386041
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
ggtAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCAGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGAtt                                                                                                                                                                                                                                                   <  1:223246/251‑1 (MQ=255)
        aaaCTATCAGCAAACCGGAAGATCTGATTGGCAACCGCATCTCCGTACCGTTTAGCTCCACCCCCCACTACAGCCTGCTGGCGGCACGGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATAGTGATGGTGCTTATGTCTGGTGACCGGCGGTTAGCTTCCTGGATAAATGCGGCAAGGTGATGACCGATACTGATCAg                                                                                                                                                                                                                                          <  2:401682/251‑1 (MQ=255)
          aCTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGACCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGt                                                                                                                                                                                                                                       <  1:175771/251‑1 (MQ=255)
             aTCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCTGTTAACGCCCTTGAATAAGACGTCAAGGTGTGGACCGATTCTGAACAGGTCg                                                                                                                                                                                                                                     <  2:647461/250‑1 (MQ=255)
             aTCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGCCTGGGCACCGGCGGCTAACCCCCAGGGCACAGACCGGCAGGGGTTGTCCGATTCTGACCAGGCCgg                                                                                                                                                                                                                                    >  1:261815/1‑251 (MQ=255)
                    aaCCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGGTATAACCCCTGGAAAAAACGGCGAGGTGTTTACCGTGTCTGaaaaagtgggggagtgg                                                                                                                                                                                                                              >  2:499248/1‑236 (MQ=255)
                             aTCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTATCACCCCGGGAAAAAACGGGAAGGGGTTGGCCGGTTTTGAACAgggtggggcgtggggggggccc                                                                                                                                                                                                                     >  1:71175/1‑231 (MQ=255)
                                                 gCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCTGCGGTTAACGCTCTGGAGAAAGACGGCAAGGTGTTGACCGATTCTTAACGGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGggtggt                                                                                                                                                                                             <  1:307177/251‑1 (MQ=255)
                                                  ccGTACCGTTTATCTCCCCCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGgtggtg                                                                                                                                                                                            <  2:53228/251‑1 (MQ=255)
                                                   cGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAACCCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACCACCCTGGACAAGACGGCAACGAGTTGCCCGATTTTGCAGCAGGCGGGGCGTGGGGGGCGCCAACGCTGGACGGCTggggggggg                                                                                                                                                                                           >  1:94990/1‑245 (MQ=255)
                                                    gTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTgcg                                                                                                                                                                                          <  2:97945/251‑1 (MQ=255)
                                                      aCCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAAACAGCAAGGGGTTGACCGGTTCTGAACAGGTCGGGCAGTGGGGGGCGCCAACGCGGGGCGGCTggggggggtgcc                                                                                                                                                                                        >  1:528175/1‑242 (MQ=255)
                                                      aCCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTgcgc                                                                                                                                                                                        <  1:425257/250‑1 (MQ=255)
                                                       ccGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGCAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTCACACCCTGGAAAAAGACGGCAAAGGGTTGAACGGTTCTGACCAGGTCGGGCCGGGGGGCGCGCCCAACCTGGACGTCCGggtggggggca                                                                                                                                                                                        >  1:582569/1‑244 (MQ=255)
                                                       ccGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAAACGGCAAGGGGTTGACCGATTCTTGACACGGCAGGCAGTGGGGCGCGCCAACGCTGGACGTCTGTGTGGGGCGCa                                                                                                                                                                                       >  1:450053/1‑250 (MQ=255)
                                                         gTTTATCTCCACCACCCACTACAGACTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGCTAAAGCCCTGGGAAAAGAAGGCAAGGGGTTGACCGATTTTGAAAAGGGTTGGGAGTGGGGGGCGCCACCGCCGGAAGGCTggggggggggcaaa                                                                                                                                                                                      >  1:593143/1‑239 (MQ=255)
                                                            tATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGGGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGAt                                                                                                                                                                                  <  2:388516/250‑1 (MQ=255)
                                                             aTCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGAtt                                                                                                                                                                                 <  1:445096/250‑1 (MQ=255)
                                                                       cccacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGGTTACGCCCCGGGAAAAAAAGGGCAGGGGGTGGACGGTTCTTGACAGGGCGGGCAGGGGGGGGGGGCCACGCTGGACGTCTGGGGGGTGGGGCAAGATTTTTCCGGGaa                                                                                                                                                                       >  1:646305/1‑250 (MQ=255)
                                                                             aCAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCGAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCt                                                                                                                                                                <  1:63643/251‑1 (MQ=255)
                                                                                 ccTGCTGGCGGCACTGAAACACTGGGGCATTAAGCCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTTGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTTGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAgg                                                                                                                                                            <  1:620204/251‑1 (MQ=255)
                                                                                       ggcggcCCTGAAACACTGAGGCATTAAACCCGGGCAAGTGGAGATAGTGGACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCATATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAGCAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGTGAAACATCCTGAGGTCGt                                                                                                                                                        <  2:556234/249‑1 (MQ=255)
                                                                                       ggcggcACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgc                                                                                                                                                                                                                         <  1:238436/187‑1 (MQ=255)
                                                                                          ggcACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAAGCCCCGGGAAAAAACGGCAAGGGGTTGACCGATTCTTGACAGGGCGGGGCGCGGGGGGCGCCAACGCCGGGCGTCTGGGGGGGGCGCAAAGATTTTGCCGGGAAAAATCCTGAGGGCGTGaaa                                                                                                                                                   >  1:226043/1‑250 (MQ=255)
                                                                                            cacggAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCg                                                                                                                                                <  2:584005/247‑1 (MQ=255)
                                                                                               tGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTCAAGCCCCGGGAAAAAAAGGGAAGGGGTTGACCGCTTCTGTACAGGGCGGGCAGTGGGGGGCGCCCACGCCGGCCGTCTGGGGGGGGCGCCAAGATTTTTCCCAGGAACATCCCTAGGGCGTGAAAGCGttt                                                                                                                                             >  1:682441/1‑250 (MQ=255)
                                                                                                gAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCg                                                                                                                                            <  2:149237/251‑1 (MQ=255)
                                                                                                gAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGGTTACGCCCCGGGAAAAGACGGCCAGGGGTTGGCCCATTCTTTACAAGGCGGGCAGTGGGGGGCGCCCACCCTGGGGCGCTGGGGGGGGCGCAAAGAGTTTGCCGGGCAACATccctggggcggggaagcggtcc                                                                                                                                             >  2:544264/1‑231 (MQ=255)
                                                                                                 aaaCACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCCGGAAAAAAAAGGCCAGGGGTTGGGCGATTCTGAACAGGGCGGGCCAGGGGGCGCGCCCACCCTGGGGGGCTGGGGGGGGCGGCAAGGATTTGCCGAGGAACAACCTGAGGGGGTGGAAGCGTTCGc                                                                                                                                           >  1:466476/1‑251 (MQ=255)
                                                                                                  aaCACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAAACGGCAAGGGGTTGGCCGATTTTGAACAGGGCGGGCAGTGGGGGGCGCCAAAGCTGGACGGCTGGGGGGGGCGCAAAGATTTTTCCGAGAAACATCCTGAGGGCGTGAAAACGGTCGCt                                                                                                                                          >  1:225480/1‑251 (MQ=255)
                                                                                                    cacTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTCTGTCTGGGCACCGGCGGATAAAGCCCCGGGAAAAAACGGCGAGGTGTTGACCGATTTTGAACAGGGCGGGCAGGGGGGGGCGCCCACGCTGGGCGTCTGGGGGGTGCGCCAAGATTTTGCCGAGAAACATCCTGAGGGCGTGAAAGCGTTCGCTa                                                                                                                                         >  1:448640/1‑250 (MQ=255)
                                                                                                                  aCCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATg                                                                                                                          <  2:83063/251‑1 (MQ=255)
                                                                                                                     cGGGCAAGTGGAGATTGGGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCTATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTc                                                                                                                       <  1:485428/251‑1 (MQ=255)
                                                                                                                         cAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGGTAACCCCCCTGAAAAAGACGGCCAGGGGGTTACCGGTTCTTCACAAGGCGGGCCGAGGGGGGCGCCAACACTGGGCGGCTGGGTGGGGCGCCAAGATTTTTCCGCGAAACATCCCGAGGGCGGGGAAAGGTTAGCTAAAAGCGCCATCGATgaccagc                                                                                                                    >  1:48454/1‑244 (MQ=255)
                                                                                                                          aaGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGTGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCaa                                                                                                                 <  1:70790/251‑1 (MQ=255)
                                                                                                                              ggAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGAAGGCGAGGTGTTGTCCGATTTCGAACAGGTCGGGCAGGGGGGCGCGCCCACCCTGGAGCGCTTGGGGGGGCGGCAAGAATTTTCCCGAGAACATCCCGAGGGCGTGTAAACGTTCCCTACAACCGCCATCgtagccctgcaacccc                                                                                                              >  2:231677/1‑236 (MQ=255)
                                                                                                                              ggAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGt                                                                                                             <  1:499248/251‑1 (MQ=255)
                                                                                                                                       gAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTa                                                                                                    <  2:466476/251‑1 (MQ=255)
                                                                                                                                                gccgccCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTTGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACg                                                                                           <  1:297807/251‑1 (MQ=255)
                                                                                                                                                               cGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGc                                                                            <  2:682441/251‑1 (MQ=255)
                                                                                                                                                               cGCTGCCTGGCAGCGGGGAGAGATCGATGGTGCTTATGTCTGGGCAGCGGCGGTGAGCGCTCGGGAAAAAGACGGCAATGTGTTGACCGATTCTGAACAGGGCGGGCATTGGGGCGCGCCAACGGTGGAGGTCTGGGTGGTGCGCAAAGATTGTGTCGAGAAACAGTCTGAGGTCGTTAAAGCGTTCGTTAAAAGCGCCATCGATGCTGAGCAATCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGc                                                                            <  1:306815/251‑1 (MQ=255)
                                                                                                                                                                          aGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTTACGCCCCGGAAAAAGACGGCAAGGGGTTGGCCGATTCTGAACAGGTCGGGCAGGGGGGCGCGCCAACGCTGGACGTCTGGGGGGTGCGGAAAGATTTTTCCGAGAAACATCCTGAGGGCGGGAAAGCGGTCGCTAAAAGCGCCCTCGATGCTCAGCAACCGTACATTGCTAACCCACACGTGTGGCTGAAACAGCCGGCAAACAt                                                                  >  1:5726/1‑250 (MQ=255)
                                                                                                                                                                           gCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAAGCCCCGGGAAAAAACGGGAAGGGGTTGACCGGTTCTGAACAGGGCGGGCCGTGGGGGGCGCCCACGCTGGGCGGCTGGGGGGGGGGCAAAGATTTTTCCGGGAAACATCCTGGGGTCGTGAAAGCGTTCGCTAAAAACGCCATCGATGCTCCGCAACCGCACATTTCTAACCCCGACGTGTGGCTGCAACAGCCGGaaaaaaacc                                                                >  1:228518/1‑246 (MQ=255)
                                                                                                                                                                             ggggAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTCGGAAAAGACGGCAAGGGGGTGAACGGTTCTGAACACGTCGGGCAGGGGGGGGCGCCCACGCTGGACGTCTGGGTGGGGGGCAAAGATTTTTGCGGGAAACATCCTGAGGGCGTGGAAGCGTTTGCTTAAAGCGCCCTCGGAGCTCAGCCACCCGCACTTTCTAACCCCGACGTGTGGCTggcaccgaccgaaaacatcacc                                                              >  2:38384/1‑230 (MQ=255)
                                                                                                                                                                                 agaTATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAAc                                                          <  1:544264/251‑1 (MQ=255)
                                                                                                                                                                                         aTGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGCAAAAAACCGCAAAGGGTTGACCGCGTCTGAAAAGGGCGGGCAGGGGGGGGCGCCCAAGCTGGACGTCCCGGGGGGGCGCAAAGAATTTGCCGGGAAACAACCTGAGGTCGGGAAAGCGTTCGCTAAAAGAGCCATCGATGCTCAGCAACCGTCCATTTCTAACCCAAGCGGGTGGGTGAAAAAACCGGAAAACATCCACAAAATGGCGCGt                                                  >  1:677612/1‑251 (MQ=255)
                                                                                                                                                                                         aTGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACAGCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAATAGCCGGAAAACATCAGCAAACTGGCGCGt                                                  <  1:30617/251‑1 (MQ=255)
                                                                                                                                                                                         aTGGTGCTTATGTCTGGGCACCGGCGGATAAAGCCCTGGGAAAAAACGGCAAGGGGTTGACCGATTTTGAACAGGGCGGGCAGTGGGGGGCGCCCACGCTGGACGTCTGGGGGGGGCGCAAAGATTTTTCCCAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTTAAAACGCCATCGATGCTCCGCAACCGTACATTGCTTACCCCGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGcgcgg                                                  >  1:530296/1‑250 (MQ=255)
                                                                                                                                                                                            gTGCTGATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAGCCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTa                                               <  2:634534/251‑1 (MQ=255)
                                                                                                                                                                                              gCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGGGGTGACCCGTTCTGAAAAGGTCGGGCAGGGGGGGGCGCCAACGCCGGACGTCTGGGGGGGGCGCCAAGATTTTTCCGAGAAACATCCTGAGGTCGTGAAAGCGGTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTTCTAACCCAGACGTGTGGCTGAAACAGCCGGGAAACATCAGCAAACTGGCGCCtttta                                              >  1:515386/1‑248 (MQ=255)
                                                                                                                                                                                                    gTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGt                                       <  2:189918/251‑1 (MQ=255)
                                                                                                                                                                                                        gggCACCGGCGGTTAACGCCCGGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGc                                     <  2:71175/249‑1 (MQ=255)
                                                                                                                                                                                                          gCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGGCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGa                                 <  2:582569/251‑1 (MQ=255)
                                                                                                                                                                                                              cggcggTTAACGCCCTGGAAAAAGACGGCAAGGGGTTGACCGGTTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAAATGGCGCGTTTAAGCGGCGGGGCTGAAggg                             >  1:156577/1‑250 (MQ=255)
                                                                                                                                                                                                                  ggTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGGGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAACGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACg                         >  1:523643/1‑251 (MQ=255)
                                                                                                                                                                                                                   gTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGt                        <  2:228518/251‑1 (MQ=255)
                                                                                                                                                                                                                   gTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCGGAACAGGTCTGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGt                        <  2:226043/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaCGCCCTGGCAAAAGACGGCAAGGTGTTGACCGATTCAGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATAGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGAAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTcc                     >  1:205621/1‑251 (MQ=255)
                                                                                                                                                                                                                        aaCGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTACGCGGCGTGCCTGAAGGTGACGTTcc                     >  2:224766/1‑251 (MQ=255)
                                                                                                                                                                                                                         aCGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTTCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAATCTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCg                    <  1:81546/251‑1 (MQ=255)
                                                                                                                                                                                                                            gCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCggg                  >  2:259662/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGAAAACTGGCGCGTTTACGCGGCGTGCCTGAAGGCGACGTTCCGGGGCTg              >  2:376791/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaaaaGCCGGCAAGGTGTTGACCGAGTCGGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTg           <  2:450053/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aaaaaGACGGCAAGGTGTTGACCGATTCTGCACAGGTCCGGCCGCGGGGCGAGCCAACGCTGGACGCCTGGGCGGTGCGACAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCCTCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTg           >  1:92312/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aaaaGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTGa          <  1:25344/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aaaGTCGGCAAGGAGTTGACCGATTCTGAACAGGACGGGCAGTGGGGCGCGGCAACGCTGGACGTCTGGGAGGTGCGCGAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATAGATGCTCAGCAACCGTGCATTGCTAACCCAGACGTGTGGCTGCAACAGCCGGAAAACATCCGCAAACTGGCGCGTTTACGCGGCGTGCCTGACGGGGTCGTTCTGGGGCTGGTGa          >  2:119096/1‑250 (MQ=255)
                                                                                                                                                                                                                                           cGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTGAAGGGGa    <  1:38384/251‑1 (MQ=255)
                                                                                                                                                                                                                                             ggCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGCGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTGAAGGGGa    >  1:592491/1‑250 (MQ=255)
                                                                                                                                                                                                                                              gCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCTTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTGAAGGGGAAt  <  1:670815/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
GGTAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTGAAGGGGAAT  >  NC_000913/385558‑386041

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: