Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,793,377 Δ1 bp coding (182/750 nt) btuD ← vitamin B12 ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,793,3770T.98.2% 255.1 / ‑1.5 54coding (182/750 nt)btuDvitamin B12 ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (28/25);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTT  >  NC_000913/1793098‑1793667
                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                   
acttGGCCTTCGGGATTGGCTTGTGGTGTGATTTGCAACACCGCCGCAGCAAGACGTACGCTTTGCCATTCACCGCCTGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAACCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTAGCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTc                                                                                                                                                                                                                                                                                <  2:241094/298‑1 (MQ=255)
aaTTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTc                                                                                                                                                                                                                                                                                <  1:174580/300‑1 (MQ=255)
  tttgCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGANAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTccc                                                                                                                                                                                                                                                                              <  1:207011/297‑1 (MQ=255)
  ttGGCCTGCGGGATTGGCTTTTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGANAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTccc                                                                                                                                                                                                                                                                              <  1:207122/300‑1 (MQ=255)
  ttGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGANAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTccc                                                                                                                                                                                                                                                                              <  1:207112/300‑1 (MQ=255)
         gCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTNCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGc                                                                                                                                                                                                                                                                       <  2:741444/300‑1 (MQ=255)
         gCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGAGTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTTAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGc                                                                                                                                                                                                                                                                       <  2:362795/300‑1 (MQ=255)
            ggATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCt                                                                                                                                                                                                                                                                      >  1:161737/1‑298 (MQ=255)
                 ggCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCAtt                                                                                                                                                                                                                                                               >  1:13216/1‑300 (MQ=255)
                             gATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGANTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTACATTACCGCTGGTCATTCCGGCCATTcgc                                                                                                                                                                                                                                                   >  1:211927/1‑300 (MQ=255)
                             gATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTACCTTACCGCTGGTCATTCCGGCCATTcgc                                                                                                                                                                                                                                                   >  2:435024/1‑300 (MQ=255)
                             gATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTACAGTGGTGCCCAGCGAACTGAATGCTTACCTTACCGCTGGTCATTCCGGCCATTcgc                                                                                                                                                                                                                                                   >  2:433932/1‑300 (MQ=255)
                                   cAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGt                                                                                                                                                                                                                                             >  1:904745/1‑300 (MQ=255)
                                             gcagcaAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActc                                                                                                                                                                                                                                   <  1:343619/299‑1 (MQ=255)
                                               agcaAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGTTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTActct                                                                                                                                                                                                                                  <  2:84446/298‑1 (MQ=255)
                                                  aaGACGTACGCGTTGCCATCCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGNCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAccc                                                                                                                                                                                                                             <  2:670544/300‑1 (MQ=255)
                                                  aaGACGTACGCGTTGCAATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCNGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAccc                                                                                                                                                                                                                             <  1:201337/300‑1 (MQ=255)
                                                     aCGTACGCGTTGCCCTTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCATGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCACGAACTGAATGCTTCCCTTACCGCTGGTCCTTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCcgc                                                                                                                                                                                                                           >  2:99402/1‑299 (MQ=255)
                                                     aCGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCcgc                                                                                                                                                                                                                           >  1:138234/1‑299 (MQ=255)
                                                      cGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCgcgc                                                                                                                                                                                                                         >  1:130792/1‑300 (MQ=255)
                                                           gcgTTGCCATTCACCGCCGGAAAGGTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCNTTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCCTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCTGCGCCATTc                                                                                                                                                                                                                    <  1:189667/300‑1 (MQ=255)
                                                              ttGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGc                                                                                                                                                                                                                 >  1:796577/1‑300 (MQ=255)
                                                                gCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCTAGTTTGTCAGCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGccc                                                                                                                                                                                                               <  2:455991/300‑1 (MQ=255)
                                                                           ccgGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCa                                                                                                                                                                                                    >  1:724465/1‑300 (MQ=255)
                                                                              gAAAGTTGATTGGTGGCTACGTCGCGAGTTGGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCTGCCATTCGCGCCAGTAAGGGTACTCTTACCCgcgc                                                                                                                                                                                                                         <  1:319632/279‑1 (MQ=255)
                                                                              gAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCATTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAgta                                                                                                                                                                                                 <  2:531466/300‑3 (MQ=255)
                                                                               aaaGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGAt                                                                                                                                                                                                <  2:844218/300‑1 (MQ=255)
                                                                                         ggTGCTACGTCCGAGTTTGTCATCGAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCTGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCTGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCTGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGcc                                                                                                                                                                                      <  2:476408/300‑1 (MQ=255)
                                                                                                                       gCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGc                                                                                                                                                          >  1:49615/1‑300 (MQ=255)
                                                                                                                          aGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCg                                                                                                                                                         <  2:883110/298‑1 (MQ=255)
                                                                                                                                  tGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTc                                                                                                                                                <  1:61635/299‑1 (MQ=255)
                                                                                                                                           atTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCt                                                                                                                                      >  1:89480/1‑300 (MQ=255)
                                                                                                                                                                tttATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGAc                                                                                                                 <  1:154055/300‑1 (MQ=255)
                                                                                                                                                                             tGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGcc                                                                                                      >  2:623452/1‑298 (MQ=255)
                                                                                                                                                                             ggCAGTGTCAGGTAGTGCCGGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGATGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCAGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCACCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCGGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCAAAATAGAAATTAtttt                                                                                                         >  2:623444/2‑295 (MQ=255)
                                                                                                                                                                                 gtgtCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCGTGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGcc                                                                                                >  1:446959/1‑300 (MQ=255)
                                                                                                                                                                                 gtgtCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGcc                                                                                                >  1:8256/1‑300 (MQ=255)
                                                                                                                                                                                 gtgtCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGc                                                                                                 >  2:2458/1‑299 (MQ=255)
                                                                                                                                                                                 gtgtCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCACAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATCATTTTGCCAACGc                                                                                                 >  2:495335/1‑299 (MQ=255)
                                                                                                                                                                                 gtgtCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCAGCGAACTGAATGCTTCCATTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGcc                                                                                                >  2:830780/1‑300 (MQ=255)
                                                                                                                                                                                      aGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGttt                                                                                           <  1:446991/300‑1 (MQ=255)
                                                                                                                                                                                           gTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGc                                                                                      >  1:41483/1‑300 (MQ=255)
                                                                                                                                                                                            tGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGNCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCt                                                                                     >  2:657528/1‑300 (MQ=255)
                                                                                                                                                                                                         ttGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTGTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCATGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATcac                                                                            <  2:458524/296‑1 (MQ=255)
                                                                                                                                                                                                             aaaCGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGAt                                                                    <  1:637815/300‑1 (MQ=255)
                                                                                                                                                                                                                                        aTAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAAc                                         <  1:365513/300‑1 (MQ=255)
                                                                                                                                                                                                                                            gcgcGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCt                                      >  2:120676/1‑299 (MQ=255)
                                                                                                                                                                                                                                                 aTGCAGCGCGAGTTTTNTTGCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGANCTTCCGGCGTCATATCCGNGGAAAAACGCTGGATGa                                >  1:786445/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                 gttgCGGACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCTGGGAAAAACGCTGGATGACTTTTCCGTCCCTGcc                <  1:794508/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ggACCATGCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCa           >  1:244131/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                             gCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCATGCCAACCAGGAAttt    >  1:883141/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                             gCTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTAACCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCAAAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTACCTGCCAACCAGGAAttt    >  1:471345/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                              cTTCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGNTGCAGGATCTCCCCAGCCCGAACCTCGCCAGANAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCCCCATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGATGGATGACTTTTCCGTCCCTGCCAAACAGGAAtttt   >  2:707857/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                               ttCCAGTGGTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTACGTCCCTGCCAACCAGGAAttttt  >  2:628916/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                   
AATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGAGCCAGCGCCCCTGCGACATCATTCAGTAGTTCGGTACGCGTTTTATCGTGCTGATGCAGTGTCAGGTAGTGCCAGACCGGCGTTGCAAACGGCGGCGTCTGCTGTTGTGAAAGATAGGCGCGATGCAGCGCGAGTTTTGTTGCGGACCATGCTTCCAGTGGTTGCCCCGCGAACTGAATGCTTCCCTTACCGCTGGTCATTCCGGCCATTCGCGCCAGTAAGGTACTCTTACCCGCGCCATTCGGCCCCACCAGGTGCAGGATCTCCCCAGCCCGAACCTCGCCAGAAAGCGGCCCCAGGCGGGTAGATTCCGCAACATCTTGTAACTGCATCACAATAGACATTATTTTGCCAACGCCAGTTTAATGCTTTCCATCACAATGGGATCTTCCGGCGTCATATCCGGGGAAAAACGCTGGATGACTTTTCCGTCCCTGCCAACCAGGAATTTTT  >  NC_000913/1793098‑1793667

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