Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 22,048 G→A L214L (CTG→CTA ribF → bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091322,0480GA98.1% 149.7 / NA 52L214L (CTG→CTAribFbifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/29);  minor base C (0/1);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGC  >  NC_000913/21807‑22289
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
aaGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTc                                                                                                                                                                                                                                            >  2:326307/1‑249 (MQ=255)
     gATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTc                                                                                                                                                                                                                                     <  1:118103/251‑1 (MQ=255)
               tattaCAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGAGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGGGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGaa                                                                                                                                                                                                                            >  2:263304/1‑250 (MQ=255)
                   aCAGAAAGCTGTCATGGAATACGGCTTCGATATCACCAATACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGAGGAATTGGGGCGACCTATAGGTTTCCCGACGGCTAATGTACCGCTACGCCGGCCGGTTTCCCCGGTGAAAgtgg                                                                                                                                                                                                                       >  2:304538/1‑248 (MQ=255)
                          gCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCAGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGAGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGGTTATGc                                                                                                                                                                                                                >  1:241268/1‑251 (MQ=255)
                             ggCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCgg                                                                                                                                                                                                              <  1:225162/250‑1 (MQ=255)
                                 tGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCCGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGa                                                                                                                                                                                                          >  2:455389/1‑250 (MQ=255)
                                   gAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGc                                                                                                                                                                                                                >  2:437307/1‑242 (MQ=255)
                                   gAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGc                                                                                                                                                                                                                <  1:437307/242‑1 (MQ=255)
                                       aCGGCTTCGATATCACCAGTACGCAAACTTTTTGGGAAGGTGGCGTGCGTAGCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCGGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCt                                                                                                                                                                                                    <  2:367974/250‑1 (MQ=255)
                                            ttCGATATCCCCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACACTCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGAGTAGTCCACGGTGATGCATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCCGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCt                                                                                                                                                                                              >  2:208432/1‑251 (MQ=255)
                                                          tACGCAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGCGCGTCAGGCGCTGGCGGATGACAATCTGGATCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCCGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGACAAGCCGt                                                                                                                                                                                >  2:272518/1‑251 (MQ=255)
                                                              cAAACTTTTTGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCGGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGTGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTCCCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTAcc                                                                                                                                                                            <  2:386737/251‑1 (MQ=255)
                                                                    ttttGCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCACGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCgg                                                                                                                                                                      >  2:1494/1‑250 (MQ=255)
                                                                        gCGAAGGTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGc                                                                                                                                                                   <  1:326307/250‑1 (MQ=255)
                                                                              gTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATc                                                                                                                                                            <  1:453134/251‑1 (MQ=255)
                                                                              gTGGCGTGCGTATCAGCAGCACCGCTGTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATc                                                                                                                                                            <  1:455389/251‑1 (MQ=255)
                                                                                gtcgtgtgtaTCAGCAGCTCCGCTGTTCGTCAGTCGCTGGCGGATGACCATCTGTCTCTCGTAGAAAGTTCACTGTTGCACCTTTTTTCTATCGCCGGGCGTGTAGTCCACGGTGATGAAGTATGGCGCGCTATATGTTGCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCTGTAGATGTGCTTGGCCTCGGTGAAAAGCCTTTACCCGGCGtcgccatcatcg                                                                                                                                                           <  2:365284/241‑12 (MQ=255)
                                                                                                   ccgctgTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTgg                                                                                                                                                                                                  >  1:478098/1‑192 (MQ=255)
                                                                                                   ccgctgTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTgg                                                                                                                                                                                                  <  2:478098/192‑1 (MQ=255)
                                                                                                        gTGCGTCAGGCGCTGGCGGATGACAATCTGGCTCTGGCAGAAAGTTTATTGGGGCACCCTTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCTCTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACTCCCAACGGTTGCCGGTAt                                                                                                                                  <  1:322533/251‑1 (MQ=255)
                                                                                                                        gggATAACAATCTGGCTCTGGCAGAAAGTTTACTGGGGCACCCGTTTGCCATCGCCGGGCGTGTAGTCCACGGTGAAGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTAAGCCGGCAGGTTGCCCCGGTGAAAGGGGTTTATGCGGTAGAAGGGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGaa                                                                                                                  <  2:89546/250‑1 (MQ=255)
                                                                                                                           gtGACAATCTGGCGCTGGCATAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTg                                                                                                               <  2:228898/250‑1 (MQ=255)
                                                                                                                              aCAATCTGGCTCTGGCAGAAAGTTTACTGTGTCACCCGTGTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTAGAGGTTTCCCGACGGCTAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGTTCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCttgccgtgcat                                                                                                            <  2:111248/251‑11 (MQ=255)
                                                                                                                                    tggctctggcAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGttt                                                                                                      >  2:254097/1‑250 (MQ=255)
                                                                                                                                    tggctctggcAGAAAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCACTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGtt                                                                                                       >  1:173636/1‑250 (MQ=255)
                                                                                                                                     ggggctggcAAAAAGTTTACGGGGGACCCCGTTTGCCATCTCCGTGCGGGTAGACCACGGTGATGAATGAGGGCGGACTATATGTTTCCCTACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTGTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCTCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCCTTTGTTa                                                                                                      <  2:164805/246‑1 (MQ=255)
                                                                                                                                              agaaagTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCCCTATAGTTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGGGATGGGCCTCGGTGAAAAGCCGTTACCCGTCGTGGCGAACATCGGAACACGCCCAACGGTTGCCGGGATTCGCCAGCCGCTGGAAGTGCATTTGTTAGATGTTGCa                                                                                             >  2:241407/5‑250 (MQ=255)
                                                                                                                                                   gTTTACTGGCGCACCCGTTTGCCATCTCCGGCCGTGTAGCCCACGGTGATGAACTAGGACGCACTATCGGCTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCAGGAAGTGCATTTGTTAGATTTTGCAATGGAc                                                                                       >  2:72497/1‑251 (MQ=255)
                                                                                                                                                       aCTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCttt                                                                                   <  1:44912/251‑1 (MQ=255)
                                                                                                                                                       aCTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCGCGGTGAAAAGCCGTTACCCGGCGTGGCAAACAACGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCtt                                                                                    <  1:241407/250‑1 (MQ=255)
                                                                                                                                                          ggggCCCCCGTTTGCCATCTCCGGGCGGTTAGTCCTCGGTGATGAATTAGGGCGCACTTTAGGGTTACCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTa                                                                                  <  2:322975/249‑1 (MQ=255)
                                                                                                                                                          ggggCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTAc                                                                                 <  1:349060/250‑1 (MQ=255)
                                                                                                                                                             gCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTc                                                                             >  1:412038/1‑251 (MQ=255)
                                                                                                                                                                  cGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGATTCTACCACTACGCCGTCAGGTTTCCCCGGTGCTAGGGGTTTATGAGGTCGAAGTGCTGGGCCTAGGTGAAAAACCTTTACCCGGCGTGGCTAACATCGGAACACGCCCACCGGTTCCCGGTATTCGCCAGCGGCTGGCAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGTCat                                                                        >  1:461422/1‑251 (MQ=255)
                                                                                                                                                                       gCCATCTCCGGGCGTGGAGTCCAAGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATAc                                                                    <  2:146883/250‑1 (MQ=255)
                                                                                                                                                                                  gCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTGACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGtgctgc                                                        <  2:442606/251‑1 (MQ=255)
                                                                                                                                                                                      gtAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTa                                                     >  2:395863/1‑250 (MQ=255)
                                                                                                                                                                                        aGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTaaaa                                                  <  1:108027/251‑1 (MQ=255)
                                                                                                                                                                                           ccACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTACATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTaaaaaa                                                <  1:160813/250‑1 (MQ=255)
                                                                                                                                                                                                 tgatgaATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGTGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAGCATCGGAACACGCCCAACGGTGGCCGGTATTCGCCAGCAGCt                                                                                                                      <  2:51433/174‑1 (MQ=255)
                                                                                                                                                                                                 tgatgaATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCt                                                                                                                      >  1:51433/1‑174 (MQ=255)
                                                                                                                                                                                                     gaATTAGGGCGCGATATAGGGTTCCCGGCGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAATTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCCGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGGGCTGCTTAAAAAAAGACGCAATGa                                     <  2:80576/251‑1 (MQ=255)
                                                                                                                                                                                                       aTTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGagc                                   <  1:100121/251‑1 (MQ=255)
                                                                                                                                                                                                          agggctccttataGGTTTCCCGACTGCTGATGTACCTCTACGCCGTCATGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTTGGTGAAAAGCCGTTTCCCTGCGTGGCAAACATCGGAACACGCCCTACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATTGACCTTTACGGTCGCCTTATTCAAGTAGTGCTGCGTAAAAAAATACGCAATGagcag                                 <  1:297172/241‑1 (MQ=255)
                                                                                                                                                                                                              cgcACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCCTATACAAGTAGTGCTGCGTAACAAAATACGCAATGAGCAGCGAtt                            >  1:453045/1‑251 (MQ=255)
                                                                                                                                                                                                                cACTATAGGTTTCCCGACGGCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTg                          <  1:395036/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCGAATGTACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGCCCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTg         >  1:186014/1‑250 (MQ=255)
                                                                                                                                                                                                                                   gCGAATGCACCGCTCCGCCGTCAGGTTTCCCCGGTGAAATGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCTGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGaa       <  2:412038/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gtaaggTACCGCTACGCGGTCCGGTTGCCCCGGTGAAAGGGGTTTATGCGGTAGAAATGGTGGGCGTCGGTGGAAAGCCGTTACCCGGCGCGGCAATCAGCGGTACACGCCCAACGGTGGCCGGTATTCGACCGCAACTGGAAGCTCAGTTGTTATATGTTGCAATGGACCTTTGCGGACGCCATATACAACTAGAGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGaaa      <  2:368859/246‑1 (MQ=255)
                                                                                                                                                                                                                                     gAATGTACCGCAAAGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGCAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGCAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGaa       >  2:497769/1‑249 (MQ=255)
                                                                                                                                                                                                                                          tACCGCTACGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTAACCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGGAGTGCTGCGTAAAAAAATACGCAATGAGCAGAGATTTGCGTCGCTGGACGAACTGAAAgcgc  >  1:488096/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
AAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGATATCACCAGTACGCAAACTTTTTGCGAAGGTGGCGTGCGCATCAGCAGCACCGCCGTGCGTCAGGCCCTTGCGGATGACAATCTGGCTCTGGCAGAGAGTTTACTGGGGCACCCGTTTGCCATCTCCGGGCGTGTAGTCCACGGTGATGAATTAGGGCGCACTATAGGTTTCCCGACGGCGAATGTACCGCTGCGCCGTCAGGTTTCCCCGGTGAAAGGGGTTTATGCGGTAGAAGTGCTGGGCCTCGGTGAAAAGCCGTTACCCGGCGTGGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGC  >  NC_000913/21807‑22289

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