Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,111,259 Δ1 bp coding (718/900 nt) rfbD ← dTDP‑4‑dehydrorhamnose reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,111,2590G.89.8% 175.0 / 13.0 50coding (718/900 nt)rfbDdTDP‑4‑dehydrorhamnose reductase
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base . (24/20);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CATTTCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTC  >  NC_000913/2111017‑2111501
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
catttcatCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGc                                                                                                                                                                                                                                           <  1:483467/251‑1 (MQ=255)
catttcatCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGc                                                                                                                                                                                                                                           >  1:256930/1‑251 (MQ=255)
        catTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAACAACTATTAAATTGCTGTAGTCGTAACTAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGCGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGTTTTgcgcg                                                                                                                                                                                                                                   >  2:304535/1‑251 (MQ=255)
                      cATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAAcc                                                                                                                                                                                                                     <  1:78249/251‑1 (MQ=255)
                           tgctcacacaTCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTTGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCa                                                                                                                                                                                                                    <  2:406668/239‑1 (MQ=255)
                           tGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAgc                                                                                                                                                                                                                  <  1:415920/249‑1 (MQ=255)
                             ccccaccaACACGAAATAGATGCAAAAACTATTAAATTGCTGTAGTCTGAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCgcag                                                                                                                                                                                                              <  2:641147/249‑1 (MQ=255)
                                    aTCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATc                                                                                                                                                                                                       <  1:304535/251‑1 (MQ=255)
                                        cGAACAAGATGCAAAAACTATTAAATGGCTGTAGTCGTATATATGTCTGTGAGCATTCGTGTGTGTCCAACCTGCCAGTCATGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTAGTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCGTTTGTTGGTACTGCGGTGAGCTTGGTGAGTGCT‑TGGGAATGCCTGCTTTGCGCGCCTCGTCAAAAACCAGCGCAGAATAATAGTaa                                                                                                                                                                                                   <  1:404140/251‑2 (MQ=255)
                                              aGATGAAAAAACTATTAAATTGCTGTCGTCGTAAATAATTCATTGAGCATTCGTTTCACACCAACCTGCCCGTCAGGCAAGACAAGCGCAAAGTTCTGCTGCAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCAGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGCG‑AATACCTGCTTTGCGCGCCTATTCAAAAACCAGCGCCGCATACACATACCAGGtt                                                                                                                                                                                             >  2:296487/1‑251 (MQ=255)
                                                aTGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAgg                                                                                                                                                                                               >  1:129919/1‑247 (MQ=255)
                                                aTGCAAAAACTATTAAATTGCTGGAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAgg                                                                                                                                                                                               <  2:129919/247‑1 (MQ=255)
                                                      aaaCGATTAAATTTCTGTAGTCGTAAATAATTCATGGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTGCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCGGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGATCTGGTTGAGTGCAAGGG‑AATGCCTGCGTTGCGCGCCTCTTCAAAAACCATCGCAGCATAATCGTACCAGGTTGTGGTAcc                                                                                                                                                                                     <  2:326619/251‑1 (MQ=255)
                                                          tATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTGCTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTAAAAAACCAGCGCAGCATAATCGTACAAGGTtgtgtgtccactg                                                                                                                                                                                 >  2:134920/1‑242 (MQ=255)
                                                                     tGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGAACCAGGTTGTGGTACCACTGGCTACCAAATg                                                                                                                                                                      >  2:298339/1‑251 (MQ=255)
                                                                         gTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTAc                                                                                                                                                                  <  1:554025/251‑1 (MQ=255)
                                                                            gTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGGTGTTGGTACTGCGTTGTGCTTGTTGAGTGCAAGGG‑AATGCCCGCTTTGAGCGCCTCTTCAAAAACCAGCGTAGACTAATAGTACCAGGTTGGAGTACCACTGGCTACCAAATGGTACa                                                                                                                                                                 >  2:459504/1‑249 (MQ=255)
                                                                               aataatTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAACTTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATAGTACCAGGTTGTGGTACCACTGGCTACAAAATGGACCAAGCCt                                                                                                                                                            >  2:70168/1‑251 (MQ=255)
                                                                                      cATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTACAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACAt                                                                                                                                                     >  1:307016/1‑251 (MQ=255)
                                                                                         taggcatgcgtttccagctaaCCTGCCAGTCAGTCAAGACAAACGCAAAGTTCTGCTGTAATTTTTCTGTATTAAGGTGAGAGTGATGTGGACGACGAGCTGGTGGAGGATAGGCTGTTGTTGGGCCTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGGGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATcc                                                                                                                                                   <  2:228966/231‑1 (MQ=255)
                                                                                            gCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGGGAGCAAGGG‑AAAGCCTGCTTTGCGCGACTTTTCAAAAAACAGCGCAGAATAATCGTACCAGGATGTGGTACCACGGGCCACCAAATGGTACAAACCTGCGACATCCGGtt                                                                                                                                               >  2:556503/1‑251 (MQ=255)
                                                                                                tCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGAACAAGCCTGAGACATCCGGTTTAtt                                                                                                                                           >  1:225375/1‑251 (MQ=255)
                                                                                                   tttCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTc                                                                                                                                          >  1:321469/1‑249 (MQ=255)
                                                                                                                                cAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGGGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAAt                                                                                                           >  1:499748/1‑251 (MQ=255)
                                                                                                                                                     aTTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAgg                                                                                                                                                                                               <  2:133851/146‑1 (MQ=255)
                                                                                                                                                     aTTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAgg                                                                                                                                                                                               >  1:133851/1‑146 (MQ=255)
                                                                                                                                                                ttAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTATTGGCGCAcc                                                                            >  2:201621/1‑250 (MQ=255)
                                                                                                                                                                  aaGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCaa                                                                          >  1:289598/1‑250 (MQ=255)
                                                                                                                                                                    ggCGAGAGTTAGGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATTCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGTTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGGTGGCGCTCCAAACt                                                                       <  2:292265/251‑1 (MQ=255)
                                                                                                                                                                       gagagaTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGGTGAGTGCAAGGG‑AATGCCTGCTTTgc                                                                                                                                                                                                                                      >  2:197359/1‑89 (MQ=255)
                                                                                                                                                                       gagaGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTgc                                                                                                                                                                                                                                      <  1:197359/89‑1 (MQ=255)
                                                                                                                                                                       gagaGTTATGTGGACGACGAGCTGGAGTAGGAGTGGCTGTTGTTGGTGCTGCGTTGAGCTTGGTGAGTG‑ATAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAACCCAGCGAAGCATAATCGTTCCAGGTTGTGGTACCACTGGCTACCAAAGGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGGACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGAt                                                                    <  2:130683/251‑1 (MQ=255)
                                                                                                                                                                        ggagTTATGTGGACGACGGGCTGGTGTAGGATAGGGTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AAGGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCTTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATc                                                                   <  2:321469/250‑1 (MQ=255)
                                                                                                                                                                        agagTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATc                                                                   >  1:603559/1‑251 (MQ=255)
                                                                                                                                                                              aTGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTtaa                                                              <  1:298339/250‑1 (MQ=255)
                                                                                                                                                                               cccgctgcgtttccgatcttAGGATTGGCGGTTGTTGGTTCTGCGTTGAGCTTGTTGAGTCCAAGGG‑AATGCCTGTTTTGCGCGCCTCTTCAAAACCCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACTCGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGAGCGTtaata                                                            <  2:408091/232‑1 (MQ=255)
                                                                                                                                                                                         gaGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAtt                                                  >  1:581785/1‑251 (MQ=255)
                                                                                                                                                                                             tGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGGGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAttcttc                                              >  1:345841/1‑251 (MQ=255)
                                                                                                                                                                                             gggTGTAGGATTGGCTGTTGTTTGTACTGCGTTGAGATTGTGGAGTGCAAGTG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTCCCCCTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCAGGTTCTGTACAATCATCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAttcttc                                              <  2:225375/250‑1 (MQ=255)
                                                                                                                                                                                               gtgtAGTATAGTCTGTTTTTGGTACGTCGTTGAGCTTGTTGAGTG‑AAAGGGAAGGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCATCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAAGTGTACAAGCCGGCGACATCCGGTTTATGCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCATCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCAc                                            <  2:269672/251‑1 (MQ=255)
                                                                                                                                                                                               gtgtAGGTTAGGCTGTTTTTGGTACTGCGTTGAGTTTTTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCATCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCAc                                            <  2:87124/251‑1 (MQ=255)
                                                                                                                                                                                               gtgtAGGATAGGCTGTTGTTGGTACTGCGTTGAGATTGTTGAGTG‑AAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCAc                                            >  1:212810/1‑251 (MQ=255)
                                                                                                                                                                                                  tAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTtaata                                                            >  1:408091/1‑232 (MQ=255)
                                                                                                                                                                                                        aggtttttgtttgtACTGCGTTGAGCTTGTTGCGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCCCCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTg                                   <  2:193846/245‑1 (MQ=255)
                                                                                                                                                                                                               tgttgGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTg                                                                                      >  1:221526/1‑193 (MQ=255)
                                                                                                                                                                                                                         gCGTTGAGTGTGTTGAGTG‑AAAGGGAATGCCTGGTTTGGTAGCCTCTGCAAAAACCAGCGCGGCATATTCGTACCAGGGTGGGGTAACACTGGCTACCAAATGGGACAAGCCTGCGACAACCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTAAGCACCTGTTTGCGCACCAAACTGATCGTTAATAACCGCTAATTCGTCACGCTCTCTTGCCAGACGTAACATCGtt                  <  2:19408/251‑1 (MQ=255)
                                                                                                                                                                                                                          cGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGttt                 >  1:670163/1‑251 (MQ=255)
                                                                                                                                                                                                                               aGCTTGTTGAGTGCAAGGG‑AAGGCCTTCTTTGCGCGCCTCTGCAAAAACCAGCGCAGCAAAATCGTACCAGGTTGGGGTACCACGTGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCGCTTTTGCCAGACGTAACAGCGTTTTGGCg            <  2:280806/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aGTGCA‑GGGTAATGCCTGCTTTGCGCGCCTCTTCACAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAttt   <  1:300166/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gTGCAAGGG‑AATCCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCTTACCAGGTTGTGGTACCACTGGCTACCAACTGGTAAACGCCTACGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAGTCATCCAGCACTACAGCACCTGTTTGCCAACCACACTGATCGTTAATAACCACTCATTCTTCACGCTCTTTTCCCCGACGTAACCTCGttatgtagcatttatttc  >  2:426386/1‑235 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
CATTTCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTC  >  NC_000913/2111017‑2111501

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