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.95.9% 213.6 / 1.5 49coding (718/900 nt)rfbDdTDP‑4‑dehydrorhamnose reductase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base . (24/23);  total (25/24)
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.93e-01

TCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTCCTTTTCC  >  NC_000913/2111021‑2111508
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         
tcatcaTTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCttt                                                                                                                                                                                                                                               <  2:208734/250‑1 (MQ=255)
              aaTTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCaa                                                                                                                                                                                                                                 <  2:370725/250‑1 (MQ=255)
                ttCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGCGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCaaaa                                                                                                                                                                                                                               >  1:221216/1‑250 (MQ=255)
                       tGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAg                                                                                                                                                                                                                          <  1:218113/248‑1 (MQ=255)
                           ccaccaTCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAATTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCgcag                                                                                                                                                                                                                     <  1:358820/249‑1 (MQ=255)
                            caccaTCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAACTAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGAGTTGAGATTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCAt                                                                                                                                                                                                                  >  1:283965/1‑251 (MQ=255)
                              ccaTCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATa                                                                                                                                                                                                                 >  2:440208/1‑250 (MQ=255)
                               caTCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAAt                                                                                                                                                                                                               >  1:359277/1‑251 (MQ=255)
                                              gCAAAAACAATTAAATTGCTGTAGTCTTAAATAATTCATTGAGCAATCGTTTTACGCCAACCTGCCAGTCAGGCAGGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTG‑AAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAgtgtgtgg                                                                                                                                                                                                <  2:468534/251‑8 (MQ=255)
                                                aaaaaCTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGt                                                                                                                                                                                               >  1:409544/1‑250 (MQ=255)
                                                     cTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTCCCAcc                                                                                                                                                                                         >  2:120717/1‑250 (MQ=255)
                                                       aTTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTgg                                                                                                                                                                                       <  1:144545/251‑1 (MQ=255)
                                                           aaTTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTAc                                                                                                                                                                                   >  1:95000/1‑251 (MQ=255)
                                                               gCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCaa                                                                                                                                                                                >  1:35866/1‑250 (MQ=255)
                                                                          aaatcatacATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGc                                                                                                                                                                     <  1:238625/242‑1 (MQ=255)
                                                                               atTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCCAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGCGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCGTCTTCAACACCAAGCGCAGCATAATCGTACCAGGTTGTGGTAACACTGGCTACCAAATGGTACAAGACTGCg                                                                                                                                                                >  1:229380/1‑250 (MQ=255)
                                                                                  cATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACAt                                                                                                                                                            >  1:294862/1‑251 (MQ=255)
                                                                                            tagttgcaccccaagcttccagcaaGGCAAGACATGCGCAAAGTTCTGCGTAAATTTTTCTGTATTAATGCGAGAGTGATGTGGACGACGAGCTGGTGTAGGATATGCTTTTTTTGGTCCTGCGTTGAGCTTGTTGAGTGCAAGGT‑AAGGCCTGCTTGGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGTTACCCCTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTAt                                                                                                                                                   <  2:332539/226‑1 (MQ=255)
                                                                                              gTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTATGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGGGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTAACAgggtg                                                                                                                                                                                                   >  2:393867/1‑197 (MQ=255)
                                                                                              gTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTtg                                                                                                                                                                                                   <  1:393867/200‑1 (MQ=255)
                                                                                                                aGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCAt                                                                                                                              <  1:402577/251‑1 (MQ=255)
                                                                                                                             aaGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGATTTATGTGGACGACGTGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAACCCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATTTGCTGTACAATc                                                                                                                 <  1:117336/251‑1 (MQ=255)
                                                                                                                                      gTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTCTGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGGGCTGTACAATAAGCccgcat                                                                                                        >  2:393997/1‑250 (MQ=255)
                                                                                                                                      gTTCGGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCcagcag                                                                                                        <  1:282431/251‑1 (MQ=255)
                                                                                                                                       ttCTGCTTAAATTTTTCTGTATTAAGGCGAGAGTAATTTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTTGTTCTGCGTTTAGCTTGTTGAGAGCTAGGT‑AATGCCGGCTTTGCGCGCCCCTTCAAAAACCAGCGCAGTATAATCGTACCAGGTTGTGGTACTACTGGCTACCAAATGGTACAAGCCTGCTACCTCCGGTTTATTCAGTGCGACCCGAATGGCATTTGCTTTACAATCCGCcagcag                                                                                                        <  2:450737/250‑1 (MQ=255)
                                                                                                                                             tgAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACCTCCGGTTTATTCAGTGCGACACGAATGGCATGGGCTGTACAATCAGCCAGCAGTACAGc                                                                                                  >  1:13370/1‑250 (MQ=255)
                                                                                                                                                aaTTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTtg                                                                                                                                                                                                   >  2:84685/1‑150 (MQ=255)
                                                                                                                                                aaTTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTtg                                                                                                                                                                                                   <  1:84685/150‑1 (MQ=255)
                                                                                                                                                       cTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGgcg                                                                                       >  1:462568/1‑251 (MQ=255)
                                                                                                                                                                   gagaGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGAt                                                                           <  2:501981/251‑1 (MQ=255)
                                                                                                                                                                       gttatgtggtgtacgaGCGGGTGGAGGATAGGCTGTTGTTGGTACTGTGTTGAGCTTGTTGAGTGCAAGGG‑AATTCCTGCTTTGCGCGCCACTTCAAAATACAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGGTGGCGCACCAAACTGATCGtt                                                                       <  2:60213/239‑1 (MQ=255)
                                                                                                                                                                                  gacgaGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGAACCAAACTGAACGTTAATAACCGCt                                                             >  2:36068/1‑250 (MQ=255)
                                                                                                                                                                                    ggatcTGGTGTAGGATAGGCTGTTGTGGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAAGCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAt                                                          <  1:252395/247‑1 (MQ=255)
                                                                                                                                                                                      aGCTGGAGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCTCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATATGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGTTCGTTAATCACCGCTAAttt                                                        >  2:328234/1‑250 (MQ=255)
                                                                                                                                                                                        cTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAAttctt                                                      >  1:289228/1‑251 (MQ=255)
                                                                                                                                                                                          ggTGTAGGATAGGCTGTTGTTGGTACTGCGTGGAGCTTGTTGAGAGCAATGG‑ATTGCCTGCTTTGCGCGCCTCTTCAAAACCCAGCGCTGCATAAGCGTTCCAGGTGGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCTGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTTGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCa                                                    <  2:378873/251‑1 (MQ=255)
                                                                                                                                                                                           gtgtAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGt                                                                                                                                                                                                            <  1:155787/98‑1 (MQ=255)
                                                                                                                                                                                           gtgtAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGt                                                                                                                                                                                                            >  2:155787/1‑98 (MQ=255)
                                                                                                                                                                                                       cTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTa                                                                                                                      <  2:173121/172‑1 (MQ=255)
                                                                                                                                                                                                       cTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTa                                                                                                                      >  1:173121/1‑172 (MQ=255)
                                                                                                                                                                                                           tgttgGTACTGCGTTGAGCTTGTTGAGGGCAAGTG‑AATGCCTGCCTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGAc                                    <  1:328234/250‑1 (MQ=255)
                                                                                                                                                                                                            gttgGTACTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGt                                  >  1:250410/1‑251 (MQ=255)
                                                                                                                                                                                                                  aCTGCGTTGAGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATc                            <  2:35866/251‑1 (MQ=255)
                                                                                                                                                                                                                           aGCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCg                   <  1:292403/251‑1 (MQ=255)
                                                                                                                                                                                                                            gCTTGTTGAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGAAGTTCAGCACCTGTTGGCGAACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCg                   >  2:38892/1‑250 (MQ=255)
                                                                                                                                                                                                                              ttgttgAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCa                                                                                                    >  2:342847/1‑167 (MQ=255)
                                                                                                                                                                                                                              ttgttgAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCa                                                                                                    <  1:342847/167‑1 (MQ=255)
                                                                                                                                                                                                                              ttgttgAGTGCAAGGG‑AATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGa                  >  1:227627/1‑249 (MQ=255)
                                                                                                                                                                                                                                         aaGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAACGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAtttcctt      >  1:63118/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aagggAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTAc                                                                                                                     >  2:414225/3‑139 (MQ=255)
                                                                                                                                                                                                                                          aagggAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTAc                                                                                                                     <  1:414225/137‑1 (MQ=255)
                                                                                                                                                                                                                                             ggAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAtttccttttcc  <  1:440208/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                         
TCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTCCTTTTCC  >  NC_000913/2111021‑2111508

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