Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913385,7970TG27.9% 160.7 / 12.4 71V189G (GTG→GGG) tauAtaurine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (14/35);  new base G (19/0);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 2.44e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCGCTGGTGGTAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTGAAGGGGAATACCTA  >  NC_000913/385549‑386046
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                            
aGCGCTGGTGGTAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgtgt                                                                                                                                                                                                                                                            <  2:666792/251‑1 (MQ=255)
  cgcTGGTGGTAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTATCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTg                                                                                                                                                                                                                                                          <  1:420925/251‑1 (MQ=255)
   gcTGGTGGTAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAAGCCCCGGGAAAAAAAGGCAAGGGGTTgg                                                                                                                                                                                                                                                         >  1:94924/1‑250 (MQ=255)
        tggtAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGGTAACGCCCTGGGAAAAGACGGCCAGGGGTTGGCCGAt                                                                                                                                                                                                                                                   >  1:711596/1‑251 (MQ=255)
               gAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCCGGGAAAAGACGGCAAGGGGGTGACCGAGTCTgta                                                                                                                                                                                                                                             >  1:622485/1‑248 (MQ=255)
                aaaaCTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCGGCTGGCGGCCCTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACa                                                                                                                                                                                                                                           <  1:462018/251‑1 (MQ=255)
                        cAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCTTCCCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCggg                                                                                                                                                                                                                                    <  2:622485/250‑1 (MQ=255)
                            aaaCCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGACGGCAAGGTGGTGACCGGTTCTGAACAgggtggggcgtg                                                                                                                                                                                                                               >  1:164197/1‑241 (MQ=255)
                                cGGAAGATCGGATGGGCAAACGCATCGCCGTACCGTTGATCTCCACCACCCACCACAGCCTGCTGGCGGCACGGAAACACGGGGGCAGTAAACCCGGGAAAGTGGAGATGGTGAACCTGCAGCCGCCCGCGAGTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgc                                                                                                                                                                                                                           <  2:47892/251‑1 (MQ=255)
                                          gATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGACGGCAAGGGGGTGACCGATTCTGAACAGGGCGGGGAGGGGGGCGCGCCAACGc                                                                                                                                                                                                                 >  1:177920/1‑251 (MQ=255)
                                            ttGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACCGCCTGGGAAAAAAAGGCCAAGTGTTGGCCCCATCTGAACAGGGCGGGCAGGGGTGGGCGCCAAccccg                                                                                                                                                                                                               >  2:558219/1‑247 (MQ=255)
                                               gCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGAc                                                                                                                                                                                                            <  2:303326/251‑1 (MQ=255)
                                                             gTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAAGCCCTGGGAAAAGACGGCAAGGGGGTGGCCGATTCTGAAAAGGGCCGGGAGGGGGGGGCGCCCACGCTGGGCGGCTGGGTGGGgcg                                                                                                                                                                                              >  1:633872/1‑251 (MQ=255)
                                                                  gTTTATGTCCACCACCCACTACAGCCTGCTGGCGGCAATGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCAGCGATTATAGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGGTTATGCCCCGGCAAAAAACTGCAAGGTGTTGACAGACGCTGTCCAGCGCTGGCAGTGGGGCGAGCCAACGCTGGACGACTCTGTGGGGCGCaacg                                                                                                                                                                                        >  2:619286/1‑249 (MQ=255)
                                                                      aTCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGAtt                                                                                                                                                                                     <  1:337052/250‑1 (MQ=255)
                                                                         tcCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTTACGCCCCGGGAAAGGAGGGCAAGGGGTTTGCCGAGTTCGTACAGGGCGGGGGGGGGGGGGCGCCAACGCGGGCCGTCTgggggggggggaaaaatttttc                                                                                                                                                                                  >  2:789173/1‑232 (MQ=255)
                                                                          ccaccaCCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGcc                                                                                                                                                                                <  2:731966/251‑1 (MQ=255)
                                                                           caccacCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGACGGCAAGGGGTTGACCGATTCTGAACAGGGCGGGCAGTGGGGGGCGCCCACGCTGGGCGGCTGGGTGGGGGGGAAAGATTTTTCCg                                                                                                                                                                               >  1:340193/1‑251 (MQ=255)
                                                                             ccacccacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACCGCCCGGACAAAGAAGGGAAGGGGGTGACCGCTTCTTAAAAGGGGGGGCAGGGGGGGGCGCCCACGCCGGGCGGCTGGGTGGGGCGCaaaaagttttccgcg                                                                                                                                                                             >  2:662187/1‑239 (MQ=255)
                                                                              cacccacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCCGGGAAAAAACGGCAAAGAGTTGAACGATTCTGAACAGGGCGGGGAGGGGGGCGCGCCAACGCCGGGCGGCTGGGGGGGGGGCGAAGATTTTGTCgagg                                                                                                                                                                            >  2:119925/1‑250 (MQ=255)
                                                                                                         gAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTTACGCCCTGGAAAAAAACGGCCAGGGGTTGGCCGATTCTGAACAGGGCGGGCCGTGGGGCGGGCCAACGCTGGACGTCTGGGTGGGGCGCCAAGATTTTTCCCGGAAACATCCTGAGGGCGGGAAAACGTTCg                                                                                                                                                 >  1:653168/1‑251 (MQ=255)
                                                                                                         gAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTc                                                                                                                                                  <  2:177920/250‑1 (MQ=255)
                                                                                                          aaaCACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGc                                                                                                                                                <  1:690214/251‑1 (MQ=255)
                                                                                                                  gggCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAGAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAgcg                                                                                                                                        <  1:306023/251‑1 (MQ=255)
                                                                                                                    gCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGcc                                                                                                                                      <  1:154324/251‑1 (MQ=255)
                                                                                                                          aaCCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGAt                                                                                                                                <  1:119925/251‑1 (MQ=255)
                                                                                                                          aaCCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGAt                                                                                                                                <  1:709752/251‑1 (MQ=255)
                                                                                                                              cGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTc                                                                                                                                                              <  1:360724/217‑1 (MQ=255)
                                                                                                                                    aGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAAc                                                                                                                      <  1:439745/251‑1 (MQ=255)
                                                                                                                                           aTTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGAGAGATATTGATGGTGCTTCTGTCTGGGCACCGGCGGTTAACGCCCCGGGAAAAGACGGCAAGGGGGTGAGCGCTTATGACAAAGTCGGGCCGAGGGGCGCGCCAACGCTGGCGGGCTCGGGGGGGGGCAAAGAATTTGGCGAGAAAAATCCTTGGGGCGCGAAAACGCTCGCTAAAAACGCCCTCCATGCTCCGCaaacgaaca                                                                                                                >  1:387718/1‑243 (MQ=255)
                                                                                                                                              gtgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCCGGAAAAAGACGGCAAGGGGTTGACCGAATCTGTACAGGGTGGGCAGGGGGGGGGGCCAACGCTGGGCGGCTGGGGGGGGCGCCAAGGATTTTCCCAGAAACATCCCCGGGGCGTGGAAGAGTTCTCCAAAAGCGCCAGCGATGCTCacccaccggacaatgg                                                                                                            >  2:226143/1‑236 (MQ=255)
                                                                                                                                                   ccTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAAcc                                                                                                       <  1:455348/251‑1 (MQ=255)
                                                                                                                                                                 cgATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTg                                                                                         <  1:589314/251‑1 (MQ=255)
                                                                                                                                                                  gATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTgn                                                                                        <  2:340193/251‑2 (MQ=255)
                                                                                                                                                                  gATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCCGGGAAAAGACGGCAAGGGGTTGACCGATTCTGAACAGGGCGGGCAGGGGGGCGCGCCAACGCTGGACGTCTGGGGGGGGGGCAAAGATTTTGCCGAGAAACAACCTGAGGGCGGGAAAGCGTTCGCTAAAAACGCCATCCATGCTCAGCAAACGGACATTGCTAAACCAGACCGGTGGCTg                                                                                         >  1:461323/1‑250 (MQ=255)
                                                                                                                                                                  gATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCCGGAAAAAGACGGCAAGGGGGTGACCGATTCTGAAAAAGTCGGGCAGGGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGGGCAAAAATTTTGCCGAGAAAAATCCCGAGGGCGGGAAAACGGTCGCTAAAAACGCCATCGATTGTCAGCCAACGTAAAATTCTAACCCAAACGGGTGGCtt                                                                                         >  1:79292/1‑249 (MQ=255)
                                                                                                                                                                      aTCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGGGTTGACCGATTCTGAAAAGGGCGGGGCGTGGGGGGCGCCAACGCTGGACGTCTGGGGGGTGCGGAAAGATTTTTCCGCGAAACATCCTGAGGGCGTGAAAGCGGTCGCTAAAAACGCCCTCGATGCTCAGCCACCGTACATTTCTAACCCCGACGGGGGGCTGAAAc                                                                                    >  1:688346/1‑250 (MQ=255)
                                                                                                                                                                         gCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGACGGGAAGGTGTTGACCGATTCTGGAAAGGGCGGGGAGGGGGGCGCGCCCAAGCTGGGCGTCTGGGGGGGGCGCAAAGATTTTGCCGAGGAACACCCTGAGGGCGGGAAAGCGGTCGGTAAAAGAGCCATCGATGCTCCGCAACCGGACATTGCTAACCCAAACGTGTGGCTGAAACaacc                                                                                >  2:517430/1‑248 (MQ=255)
                                                                                                                                                                           tGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGACGGGAAGGGGGTGACCGAATCTGAAAAGGGCGGGCAGTGGGGGGCGCCCACGCTGGGCGTCTGGGTGGGGCGCAAAGATTTTTTCGAGAAACAACCTGAGGGCGGGGAAACGGTCGCGAAAAACGCCATCCATGCTCCGCAAACGGACATTTCTAAACCAAACGAGTGGGTGGAACAGACgg                                                                              >  2:728974/1‑251 (MQ=255)
                                                                                                                                                                             ccTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCCGGAAAAAAACGGCCAGGGGTTGAACGGTTCTGGACAGGTCGGGCGGTGGGGGGCGCCCACGCTGGGCGTCTGGGGGGGGCGCAAAAAATTTGCCGAGGAAAATCCTGAGGGCGTGGAAGAGTTCTCTAAAAGCGCCCTCGATGCTCAGCAACACTACATTTCTAACCCAGACCTGTGGCTGGAACAGCCggc                                                                             >  2:737416/1‑249 (MQ=255)
                                                                                                                                                                                   aGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAAAGCCCTGGAAAAAAAAGGCAAAGGGGTGACCGATTCTTAAAAGGGCGGGCCGTGGGGCGCGCCAACGCCGGACGTCTGGGTGGGGCGCAAAGAATTTGCCGAGAAAAATCCTGAGGGCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAAACACCGGGAAACATc                                                                      >  1:681042/1‑251 (MQ=255)
                                                                                                                                                                                     cGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGGAAAAGACGGCCAGGGGGTGGCCGATTCTGAACAGGGCGGGGAGGGGGGGGCGCCAAAGCTGGACGGCTGGGGGGGGGGCCAAGATTTTGGCGGGGAAAATCCTGAGGGCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGGAAACGGACAATGCTAACCCAGAAGGGGGGCTGAAAAAGCCGGAAAAAATcc                                                                     >  1:727586/1‑249 (MQ=255)
                                                                                                                                                                                     cGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGCGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACan                                                                                                               <  2:363799/209‑2 (MQ=255)
                                                                                                                                                                                            atatTGATGGTGCTTATGTCTGGGCACCGGCGGTGAACGCGCTGGAAAAAGGCGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCGAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTg                                                             <  2:21931/251‑1 (MQ=255)
                                                                                                                                                                                                    ggTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGttt                                                     <  1:342916/251‑1 (MQ=255)
                                                                                                                                                                                                     gTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTa                                                    <  1:737416/251‑1 (MQ=255)
                                                                                                                                                                                                             gTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAAc                                                               <  1:447770/232‑1 (MQ=255)
                                                                                                                                                                                                             gTCTGGGCACCGGCGGTTAACGCCATGGAAAAAAACGGGCAGGGGTTGACCGATTTTTAAAAGGGCGGGGACAGGGGCGCGCCAACGCTCGGCGCCTGGGGGGTGCGCAAAGAGTTTGCGGAGAAACACCCTGAGGTGGGGGAAACGCTCGCTAAAAACGCCATCGAAGCTCAGCAACCGTACATTGCTAACCCAGAGGAGGGGCTGAAACACACGGaaaaaacccgcaaaa                                                               >  2:447770/1‑221 (MQ=255)
                                                                                                                                                                                                                      ccGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGaaa                                                                                                                                                        >  1:153040/1‑135 (MQ=255)
                                                                                                                                                                                                                      ccGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGaaa                                                                                                                                                        <  2:153040/135‑1 (MQ=255)
                                                                                                                                                                                                                        ggcggTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGGGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAATCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAACCAGCCGGAAAACATCAGCAAAACGGCGCGTTTAAGCGGAGTGCCTGAAGGTg                                 >  2:775516/1‑251 (MQ=255)
                                                                                                                                                                                                                         gcggTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGa                                <  1:330238/251‑1 (MQ=255)
                                                                                                                                                                                                                         gcggTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGa                                <  1:488181/251‑1 (MQ=255)
                                                                                                                                                                                                                          cggTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGAc                               <  1:22070/251‑1 (MQ=255)
                                                                                                                                                                                                                           ggTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAgc                                                                                                                                         <  2:193808/145‑1 (MQ=255)
                                                                                                                                                                                                                           ggTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAgc                                                                                                                                         >  1:193808/1‑145 (MQ=255)
                                                                                                                                                                                                                            gTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGc                                                                                                                         >  2:600879/1‑160 (MQ=255)
                                                                                                                                                                                                                            gTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGc                                                                                                                         <  1:600879/160‑1 (MQ=255)
                                                                                                                                                                                                                            gTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGt                             >  1:219611/1‑251 (MQ=255)
                                                                                                                                                                                                                               aaCGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTcc                          >  1:317206/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ccTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTaaa                                                                                                                                            <  2:473574/133‑1 (MQ=255)
                                                                                                                                                                                                                                     ccTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTaaa                                                                                                                                            >  1:473574/1‑133 (MQ=255)
                                                                                                                                                                                                                                        ggAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTgg                  >  1:537275/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTg                   <  2:204093/250‑1 (MQ=255)
                                                                                                                                                                                                                                         gAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGt                 <  2:69510/251‑1 (MQ=255)
                                                                                                                                                                                                                                          naaaaGACGGCAAGGAGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGCACGTCTGGGTGGCGCGCAAAGATTTTGCCGAGAAACATCCTGAGGGCGTGAAAGCGTTGGCTAAAAGCGCCATCGATGCTCAGAAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAACACATCAGCACACTGGAGAGTTTAAGCGGCGGGCCTGAAGGAGACGTTCCGGGGCTGGTg                >  2:357842/2‑251 (MQ=255)
                                                                                                                                                                                                                                           aaaaGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTg                >  2:703700/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   gCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTg                   >  2:162820/1‑239 (MQ=255)
                                                                                                                                                                                                                                                   gCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTg                   <  1:162820/239‑1 (MQ=255)
                                                                                                                                                                                                                                                   gCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTGAAGGGGAAt       <  1:728974/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       ggTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTGAAGGGGAATACCt   <  1:226143/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       ggTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACGGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGCTGGTGAAGGGGAATACCTa  <  2:13308/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                            
AGCGCTGGTGGTAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGCCGGAAAACATCAGCAAACTGGCGCGTTTAAGCGGCGTGCCTGAAGGTGACGTTCCGGGGCTGGTGAAGGGGAATACCTA  >  NC_000913/385549‑386046

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