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.97.4% 168.1 / ‑4.3 38coding (718/900 nt)rfbDdTDP‑4‑dehydrorhamnose reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (19/18);  minor base C (0/1);  total (19/19)
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

TCATTTCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTCCTTTT  >  NC_000913/2111016‑2111506
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
tcatttcatCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTg                                                                                                                                                                                                                                                 >  2:115075/1‑251 (MQ=255)
  atttcatCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCt                                                                                                                                                                                                                                               >  1:37520/1‑251 (MQ=255)
      catcatTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCTCGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATGAAGGCGAGAGTTAGGTGGACGACTAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGATCTTGTTGAGTGCAAGGGAATGCCTGCTTTg                                                                                                                                                                                                                                            <  1:319337/250‑1 (MQ=255)
      caacatTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTg                                                                                                                                                                                                                                            <  2:135000/247‑1 (MQ=255)
         catTCCTTTTAATTTATCTTGCTCCACCAGCTCGAACAAGATGCAAAATCTATTAAATTTCTGTAGTCGTAAATAATTCATTGAGCATGCGGTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTGCTGTTGAAATTTTTCTGTATTAAGGCGAGAGTTATTTGGACGACGAGCTGGTGTAGGATATGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTgcgcg                                                                                                                                                                                                                                        <  1:324477/251‑1 (MQ=255)
         catTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTCGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGCGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGACTGCCTGCTTTgcgcg                                                                                                                                                                                                                                        >  2:200659/1‑251 (MQ=255)
              cTTTTAATTCCTCTTGCTCCACCATCACGAACAAGATGCAAACACTATTAAATTGCTGTAGTCGTCAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGCCAAGCGCCATGTTATGCTGACATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGCTAGGCTGTTGTTGGTAATCGGTTGCGCTTGTTGCGTGCAAGGGCACACCTGCTTTGCGCGCcttt                                                                                                                                                                                                                                   >  2:154639/1‑249 (MQ=255)
                 ttAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTc                                                                                                                                                                                                                                 <  1:321935/250‑1 (MQ=255)
                  tAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGATCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGTCGACGAGCTGTTGTAGGATATGCTGTTGTTGGTACTGCGTTAAGCTTTTTTGAGGCAAGGGAATGCCGTCTTTGCGCCCCTCTTcc                                                                                                                                                                                                                                >  2:103128/1‑249 (MQ=255)
                      tCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAAc                                                                                                                                                                                                                           >  2:151464/1‑251 (MQ=255)
                                ccaccaTAACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAACTAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACACGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGACGAGCGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGGGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCgaagc                                                                                                                                                                                                                  >  2:194862/1‑250 (MQ=255)
                                         cGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTCTTAAGGCGAGAGTTATGTGGCCGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAACCCAGCgcagcgtatcgga                                                                                                                                                                                                          >  2:126452/1‑241 (MQ=255)
                                         cGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCgcagc                                                                                                                                                                                                                  <  1:126452/241‑1 (MQ=255)
                                          gAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTAcc                                                                                                                                                                                                       >  1:217023/1‑251 (MQ=255)
                                                             ttAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGc                                                                                                                                                                                    >  1:119278/1‑251 (MQ=255)
                                                                                   aatTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCg                                                                                                                                                              >  1:248366/1‑251 (MQ=255)
                                                                                          tGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATcc                                                                                                                                                        <  1:344604/250‑1 (MQ=255)
                                                                                            aGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGt                                                                                                                                                     >  2:286679/1‑251 (MQ=255)
                                                                                                      tCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGt                                                                                                                                            <  1:286679/250‑1 (MQ=255)
                                                                                                       cACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTGCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGc                                                                                                                                          <  2:37520/251‑1 (MQ=255)
                                                                                                            aaaCCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTTCGTAAATTTTTCTGTAGTAAGGCTAGAGTGATGTGGAAGACTAGCTGGTGAAGGAGATGCTGTGGTTGGTGCTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTGCAAAAACCAGCGCACCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGacac                                                                                                                                     <  2:140523/250‑1 (MQ=255)
                                                                                                                  gCCAGTCAGGCAAGACAAGCGCAAAGTTCGGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATg                                                                                                                                <  1:114622/250‑1 (MQ=255)
                                                                                                                            cAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTCTTAAGGCGAGAGTTATGTGGCCGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGACTGGCATGTGCTGt                                                                                                                     >  2:349195/1‑251 (MQ=255)
                                                                                                                            cAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCCCTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGt                                                                                                                     >  2:349191/1‑251 (MQ=255)
                                                                                                                            cAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGt                                                                                                                     >  1:99686/1‑251 (MQ=255)
                                                                                                                              aGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTAc                                                                                                                   <  1:151464/251‑1 (MQ=255)
                                                                                                                                          aGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATcagccagc                                                                                                        >  1:14019/1‑250 (MQ=255)
                                                                                                                                          aGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCcagca                                                                                                       <  2:14019/251‑1 (MQ=255)
                                                                                                                                                                   aaggggggagTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTTCCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCaa                                                                               <  2:217023/243‑1 (MQ=255)
                                                                                                                                                                    aGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAAc                                                                             >  1:78171/1‑251 (MQ=255)
                                                                                                                                                                       cGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTTTTGGTACTGCGTTGAGCTGGTTGAGTGCAAGGGAATGCCTGCTTGGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGa                                                                          <  2:119278/251‑1 (MQ=255)
                                                                                                                                                                                tgtgGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACg                                                                                                                                    >  1:179452/1‑184 (MQ=255)
                                                                                                                                                                                tgtgGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACg                                                                                                                                    <  2:179452/184‑1 (MQ=255)
                                                                                                                                                                                    gacgacGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCg                                                             <  1:92166/251‑1 (MQ=255)
                                                                                                                                                                                                     ggagagtctgttggttgatcTGCGTTGACCTTGTTGCCGTCAAGGGCATGCCTGCGTTGTGCGACTCTTCTATCACCAGCGCAGCTTAATCGCGCCAGGTTGTGGTATCACTGGCTACAAAATGGTACAAGCCTGCGACATTCGGTTTATTCAGTGCGACACGAATGGCATGCGCGGTACAATCAGCCTGCAGTTCAGCACCTGTTAGCGCAACAAACTGATCGTTAATAACCGCGAATTCTTCACGctct                                            <  2:8560/232‑1 (MQ=255)
                                                                                                                                                                                                         aGGCTGTTGTTGGTACTGCGGTGAGCTTGTTGAGTGCAAGGGAATGCCTGCTTTTCGCGCCTCTTCAAAAACCAGCGCAGCATAAGCGTTCCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACGCGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTg                                        <  1:318844/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tgttgAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAACAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGCTTTGGCGAAGt             >  2:82019/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gttgAGTGCAAGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGtt            <  1:200659/251‑1 (MQ=255)
                                                                                                                                                                                                                                                caagggAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAtttcctt    <  1:116718/248‑1 (MQ=255)
                                                                                                                                                                                                                                                  agggAATGCCTGCTTTGCTCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAtttcctgtt  <  1:349195/250‑4 (MQ=255)
                                                                                                                                                                                                                                                  agggAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTAtttcctttt  <  1:349191/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
TCATTTCATCATTCCTTTTAATTCATCTTGCTCCACCATCACGAACAAGATGCAAAAACTATTAAATTGCTGTAGTCGTAAATAATTCATTGAGCATTCGTTTCACGCCAACCTGCCAGTCAGGCAAGACAAGCGCAAAGTTCTGCTGAAATTTTTCTGTATTAAGGCGAGAGTTATGTGGACGACGAGCTGGTGTAGGATAGGCTGTTGTTGGTACTGCGTTGAGCTTGTTGAGTGCAAGGGGAATGCCTGCTTTGCGCGCCTCTTCAAAAACCAGCGCAGCATAATCGTACCAGGTTGTGGTACCACTGGCTACCAAATGGTACAAGCCTGCGACATCCGGTTTATTCAGTGCGACACGAATGGCATGTGCTGTACAATCAGCCAGCAGTTCAGCACCTGTTGGCGCACCAAACTGATCGTTAATAACCGCTAATTCTTCACGCTCTTTTGCCAGACGTAACATCGTTTTGGCGAAGTTATTTCCTTTT  >  NC_000913/2111016‑2111506

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