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
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
catttcatCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGc                                                                                                                                                                                                                                           <  1:483467/251‑1 (MQ=255)
catttcatCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGc                                                                                                                                                                                                                                           >  1:256930/1‑251 (MQ=255)
        catTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAACAACTATTAAATTGCTGTAGTCGTAACTAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGCGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGTTTTgcgcg                                                                                                                                                                                                                                   >  2:304535/1‑251 (MQ=255)
                      cATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAAcc                                                                                                                                                                                                                     <  1:78249/251‑1 (MQ=255)
                           tgctcacacaTCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTTGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCa                                                                                                                                                                                                                    <  2:406668/239‑1 (MQ=255)
                           tGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAgc                                                                                                                                                                                                                  <  1:415920/249‑1 (MQ=255)
                             ccccaccaACACGAAATAGATGCAAAAACTATTAAATTGCTGTAGTCTGAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCgcag                                                                                                                                                                                                              <  2:641147/249‑1 (MQ=255)
                                    aTCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATc                                                                                                                                                                                                       <  1:304535/251‑1 (MQ=255)
                                        cGAACAAGATGCAAAAACTATTAAATGGCTGTAGTCGTATATATGTCTGTGAGCATTCGTGTGTGTCCAACCTGCCAGTCATGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTAGTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCGTTTGTTGGTACTGCGGTGAGCTTGGTGAGTGCTTGGGAATGCCTGCTTTGCGCGCCTCGTCAAAAACCAGCGCAGAATAATAGTaa                                                                                                                                                                                                   <  1:404140/251‑2 (MQ=255)
                                              aGATGAAAAAACTATTAAATTGCTGTCGTCGTAAATAATTCATTGAGCATTCGTTTCACACCAACCTGCCCGTCAGGCAAGACAAGCGCAAAGTTCTGCTGCAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCAGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGCGAATACCTGCTTTGCGCGCCTATTCAAAAACCAGCGCCGCATACACATACCAGGtt                                                                                                                                                                                             >  2:296487/1‑251 (MQ=255)
                                                aTGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAgg                                                                                                                                                                                               >  1:129919/1‑247 (MQ=255)
                                                aTGCAAAAACTATTAAATTGCTGGAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAgg                                                                                                                                                                                               <  2:129919/247‑1 (MQ=255)
                                                      aaaCGATTAAATTTCTGTAGTCGTAAATAATTCATGGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTGCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCGGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGATCTGGTTGAGTGCAAGGGAATGCCTGCGTTGCGCGCCTCTTCAAAAACCATCGCAGCATAATCGTACCAGGTTGTGGTAcc                                                                                                                                                                                     <  2:326619/251‑1 (MQ=255)
                                                          tATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTGCTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTAAAAAACCAGCGCAGCATAATCGTACAAGGTtgtgtgtccactg                                                                                                                                                                                 >  2:134920/1‑242 (MQ=255)
                                                                     tGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGAACCAGGTTGTGGTACCACTGGCTACCAAATg                                                                                                                                                                      >  2:298339/1‑251 (MQ=255)
                                                                         gTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTAc                                                                                                                                                                  <  1:554025/251‑1 (MQ=255)
                                                                            gTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGGTGTTGGTACTGCGTTGTGCTTGTTGAGTGCAAGGGAATGCCCGCTTTGAGCGCCTCTTCAAAAACCAGCGTAGACTAATAGTACCAGGTTGGAGTACCACTGGCTACCAAATGGTACa                                                                                                                                                                 >  2:459504/1‑249 (MQ=255)
                                                                               aataatTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAACTTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATAGTACCAGGTTGTGGTACCACTGGCTACAAAATGGACCAAGCCt                                                                                                                                                            >  2:70168/1‑251 (MQ=255)
                                                                                      cATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTACAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACAt                                                                                                                                                     >  1:307016/1‑251 (MQ=255)
                                                                                         taggcatgcgtttccagctaaCCTGCCAGTCAGTCAAGACAAACGCAAAGTTCTGCTGTAATTTTTCTGTATTAAGGTGAGAGTGATGTGGACGACGAGCTGGTGGAGGATAGGCTGTTGTTGGGCCTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGGGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATcc                                                                                                                                                   <  2:228966/231‑1 (MQ=255)
                                                                                            gCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGGGAGCAAGGGAAAGCCTGCTTTGCGCGACTTTTCAAAAAACAGCGCAGAATAATCGTACCAGGATGTGGTACCACGGGCCACCAAATGGTACAAACCTGCGACATCCGGtt                                                                                                                                               >  2:556503/1‑251 (MQ=255)
                                                                                                tCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGAACAAGCCTGAGACATCCGGTTTAtt                                                                                                                                           >  1:225375/1‑251 (MQ=255)
                                                                                                   tttCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTc                                                                                                                                          >  1:321469/1‑249 (MQ=255)
                                                                                                                                cAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGGGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAAt                                                                                                           >  1:499748/1‑251 (MQ=255)
                                                                                                                                                     aTTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAgg                                                                                                                                                                                               <  2:133851/146‑1 (MQ=255)
                                                                                                                                                     aTTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAgg                                                                                                                                                                                               >  1:133851/1‑146 (MQ=255)
                                                                                                                                                                ttAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTATTGGCGCAcc                                                                            >  2:201621/1‑250 (MQ=255)
                                                                                                                                                                  aaGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCaa                                                                          >  1:289598/1‑250 (MQ=255)
                                                                                                                                                                    ggCGAGAGTTAGGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATTCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGTTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGGTGGCGCTCCAAACt                                                                       <  2:292265/251‑1 (MQ=255)
                                                                                                                                                                       gagagaTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGGTGAGTGCAAGGGAATGCCTGCTTTgc                                                                                                                                                                                                                                      >  2:197359/1‑89 (MQ=255)
                                                                                                                                                                       gagaGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTgc                                                                                                                                                                                                                                      <  1:197359/89‑1 (MQ=255)
                                                                                                                                                                       gagaGTTATGTGGACGACGAGCTGGAGTAGGAGTGGCTGTTGTTGGTGCTGCGTTGAGCTTGGTGAGTGATAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAACCCAGCGAAGCATAATCGTTCCAGGTTGTGGTACCACTGGCTACCAAAGGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGGACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGAt                                                                    <  2:130683/251‑1 (MQ=255)
                                                                                                                                                                        ggagTTATGTGGACGACGGGCTGGTGTAGGATAGGGTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAAGGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCTTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATc                                                                   <  2:321469/250‑1 (MQ=255)
                                                                                                                                                                        agagTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATc                                                                   >  1:603559/1‑251 (MQ=255)
                                                                                                                                                                              aTGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTtaa                                                              <  1:298339/250‑1 (MQ=255)
                                                                                                                                                                               cccgctgcgtttccgatcttAGGATTGGCGGTTGTTGGTTCTGCGTTGAGCTTGTTGAGTCCAAGGGAATGCCTGTTTTGCGCGCCTCTTCAAAACCCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACTCGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGAGCGTtaata                                                            <  2:408091/232‑1 (MQ=255)
                                                                                                                                                                                         gaGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAtt                                                  >  1:581785/1‑251 (MQ=255)
                                                                                                                                                                                             tGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGGGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAttcttc                                              >  1:345841/1‑251 (MQ=255)
                                                                                                                                                                                             gggTGTAGGATTGGCTGTTGTTTGTACTGCGTTGAGATTGTGGAGTGCAAGTGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTCCCCCTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCAGGTTCTGTACAATCATCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAttcttc                                              <  2:225375/250‑1 (MQ=255)
                                                                                                                                                                                               gtgtAGTATAGTCTGTTTTTGGTACGTCGTTGAGCTTGTTGAGTGAAAGGGAAGGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCATCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAAGTGTACAAGCCGGCGACATCCGGTTTATGCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCATCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCAc                                            <  2:269672/251‑1 (MQ=255)
                                                                                                                                                                                               gtgtAGGTTAGGCTGTTTTTGGTACTGCGTTGAGTTTTTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCATCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCAc                                            <  2:87124/251‑1 (MQ=255)
                                                                                                                                                                                               gtgtAGGATAGGCTGTTGTTGGTACTGCGTTGAGATTGTTGAGTGAAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCAc                                            >  1:212810/1‑251 (MQ=255)
                                                                                                                                                                                                  tAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTtaata                                                            >  1:408091/1‑232 (MQ=255)
                                                                                                                                                                                                        aggtttttgtttgtACTGCGTTGAGCTTGTTGCGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCCCCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTg                                   <  2:193846/245‑1 (MQ=255)
                                                                                                                                                                                                               tgttgGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTg                                                                                      >  1:221526/1‑193 (MQ=255)
                                                                                                                                                                                                                         gCGTTGAGTGTGTTGAGTGAAAGGGAATGCCTGGTTTGGTAGCCTCTGCAAAAACCAGCGCGGCATATTCGTACCAGGGTGGGGTAACACTGGCTACCAAATGGGACAAGCCTGCGACAACCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTAAGCACCTGTTTGCGCACCAAACTGATCGTTAATAACCGCTAATTCGTCACGCTCTCTTGCCAGACGTAACATCGtt                  <  2:19408/251‑1 (MQ=255)
                                                                                                                                                                                                                          cGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGttt                 >  1:670163/1‑251 (MQ=255)
                                                                                                                                                                                                                               aGCTTGTTGAGTGCAAGGGAAGGCCTTCTTTGCGCGCCTCTGCAAAAACCAGCGCAGCAAAATCGTACCAGGTTGGGGTACCACGTGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCGCTTTTGCCAGACGTAACAGCGTTTTGGCg            <  2:280806/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aGTGCAGGGTAATGCCTGCTTTGCGCGCCTCTTCACAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAttt   <  1:300166/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gTGCAAGGGAATCCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCTTACCAGGTTGTGGTACCACTGGCTACCAACTGGTAAACGCCTACGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAGTCATCCAGCACTACAGCACCTGTTTGCCAACCACACTGATCGTTAATAACCACTCATTCTTCACGCTCTTTTCCCCGACGTAACCTCGttatgtagcatttatttc  >  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: