Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,873,316 T→C T8A (ACA→GCA)  ygbE ← DUF3561 domain‑containing inner membrane protein YgbE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,873,3160TC97.4% 118.8 / NA 39T8A (ACA→GCA) ygbEDUF3561 domain‑containing inner membrane protein YgbE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (17/21);  minor base A (1/0);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

GCGTAGCTTTCCGTCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGTAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTT  >  NC_000913/2873086‑2873559
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
gcgtagcATTCCGTCCATCAATGAATTCAGCGCAATCCCCAGCACTACCTCTACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATATAGAAGAACAGCGTGTTAGAGACATAAATCAACATCGGCATCGCCAGCGCAAATAGCCAGGAGATAAAACCGTCCTCGGCACCAGGCAGTGACCATGTTGGTTCTTCATCGTCATTTAGGCTGTCGTTATTTGTTAGTGCTATGTTATTGGTATGACGCa                                                                                                                                                                                                                                 <  2:33283/251‑1 (MQ=255)
  gtagcTTTCCGTCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCATGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAGCCGACCACGGCTCCAGTCAGTTAGCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATTTTATGGCTATTACGCata                                                                                                                                                                                                                               <  1:29285/251‑1 (MQ=255)
              ccATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCATCGCAAATAACCAGGAGATAAAACCGACCACGGCACCATGCAGTGACCATGTGGTTTCTTCATCCTCAGTAATGCTGTCGTTATGTGTTATTTCAATGTTATGGCTATTACGCATATTTGATCCTGtt                                                                                                                                                                                                                   <  1:182056/251‑1 (MQ=255)
                       aatgCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCATGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATATCCAGGAGATAAAACCGACCGCGGCACCAGGCAGTGTCCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATTTTATGGCGATTACGCATATTTGATCCTGTTACTTTGACg                                                                                                                                                                                                          <  1:455982/251‑1 (MQ=255)
                       aatgCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACg                                                                                                                                                                                                          >  1:284405/1‑251 (MQ=255)
                                     cccACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAATGAAGAACAGCGTGTTAGATCCATAAAGCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCTTCCACGGCACCATTCAGTTACCATGTGGTTTCTTCATCCTCAGTAATGCTGTCGTTGTTTGTTAGTGCAATGTTATGGCTATTACGTATAtttt                                                                                                                                                                                                                           <  1:372625/220‑2 (MQ=255)
                                     cccACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTg                                                                                                                                                                                                                           <  1:372628/220‑1 (MQ=255)
                                     cccACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTg                                                                                                                                                                                                                           >  2:372628/1‑220 (MQ=255)
                                       caccacTACAGCGACGGGCATCAGCGCCACAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCCAATAACCAGGCGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTATCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGCTCCTGTTACTTTGACGACCAGGGAATGGAAAc                                                                                                                                                                                          >  2:208032/1‑251 (MQ=255)
                                                   gACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGTTAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGtct                                                                                                                                                                              >  1:191468/1‑251 (MQ=255)
                                                    aCGGGCAGCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTTGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGtctc                                                                                                                                                                             <  1:251315/251‑1 (MQ=255)
                                                                     aaagaaagGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCAGAAATCACCAGCGGCATCTCCAGCGCAAATAGCCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCGTCATCCTCAGTAAGGCTGTCGTTATTTGTTATTGCAATGTTAGGGCGATGACGCATATTTGATCCTTGTGCTTTGACTAACCGGGGTTGGAAATCCGCTTGTGGCGCAGGAGCTGATAATATCg                                                                                                                                                            <  2:204908/251‑1 (MQ=255)
                                                                          aagGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTg                                                                                                                                                       <  2:361092/251‑1 (MQ=255)
                                                                           agGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGATTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGTAAACATGGTGGTGTCTCAGGATCTGATAATATCGTTcgtt                                                                                                                                                      >  2:376661/1‑248 (MQ=255)
                                                                             gCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGtct                                                                                                                                                    >  1:221415/1‑251 (MQ=255)
                                                                                 ggTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAACTCAACAACGGCATCGCCAGCGCAAATAGCCAGTAGAGAAAACCGACCACGTCACCAGGCATTGACCATGTGGTTTCTGTATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTGAGGGCGATTACGCATATTTGATCCTGTTACTTTGAAGAACCGGGGATGGAAACCCGGTGGTGTCGCATGATCTGAAAATAGCGTTCTTTCTCAAc                                                                                                                                                <  2:396645/251‑1 (MQ=255)
                                                                                      aGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGAGCAAATAACCAGGCGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGATTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGAGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAAc                                                                                                                                                >  1:423802/1‑246 (MQ=255)
                                                                                               agaaCAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTg                                                                                                                                  <  1:115727/251‑1 (MQ=255)
                                                                                                    aGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGt                                                                                                                              >  1:278977/1‑250 (MQ=255)
                                                                                                       gtgtTAGATACATAAATCAACATAGGCATCGCCCGCGCAAATAACCAGGAGATAAAACCAACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGCTAGTGAAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGACCCGGGCCTGGAAACCCGGTGGTGTCTCATGATCTGACAATATCGTTCTGTCTCAACAGATCTAATAATTGCTTTAcc                                                                                                                           >  2:248650/1‑250 (MQ=255)
                                                                                                       gtgtTAGAGCCATAAATCAACAGCGGCATCGCCAGCGCAAATAACCAGGAGATAAATCCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGGAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACATATCTAATAATTGCTGTACCa                                                                                                                          <  2:374984/251‑1 (MQ=255)
                                                                                                         gtTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCaa                                                                                                                         >  1:383074/1‑250 (MQ=255)
                                                                                                             gagCCATAAATCAACATCGGAACAGCCAGCGCAAATAACCAGTAGATAAAACCTACCACGGCACCAGGCAGTGACCATGTGGTTTCTGCAGCCTCAGTAAGGCTGTCGTTATTTGTTACTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACTAACCGTGCATGGAGACCCTGTGGTGTCTCAGGATCTGATAAGATCGTTCTGTCTCAATAGATCTAATAATTGCTGTTCCaaaattc                                                                                                                    <  2:239201/251‑5 (MQ=255)
                                                                                                                  ataaagccaCATCGGCATCGCCAGCGCAAATAACCAGAAGATAAAACCGACCTCGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTTTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATGTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCGGTCTCAACAGATCTAATAATTGCTGTCCCAAATTTGTTACt                                                                                                               <  1:287944/243‑1 (MQ=255)
                                                                                                                   aaaaTCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTa                                                                                                              <  1:62762/250‑1 (MQ=255)
                                                                                                                      aTCAACATCGGCATCGCCAGCGCACATAACCAGGCGCTAACACCAACCACGGCCCCAGGCAGTGACCATGTGGTTTCTTCATCCTCCGTAAGGCTGTCGTTATTTGTTAGTGCACTGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAAtt                                                                                                           >  1:141079/1‑251 (MQ=255)
                                                                                                                            aTCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCTCGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTGGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGAGCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCAc                                                                                                     <  2:114297/251‑1 (MQ=255)
                                                                                                                            aTCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCAc                                                                                                     >  1:189884/1‑251 (MQ=255)
                                                                                                                                gCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCAtt                                                                                                 >  1:195947/1‑251 (MQ=255)
                                                                                                                                     gCCAGCGCAAATAACCATGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTTGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTAGTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATTGAAACACGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGAt                                                                                            <  2:286251/251‑1 (MQ=255)
                                                                                                                                      ccAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATg                                                                                           >  1:329316/1‑251 (MQ=255)
                                                                                                                                           gcAAATAACCAGTAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATTTTATGGCTATTACGCATATTTTATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGAGCTAATAATTGCTGTACCAAATTTGTTACGTATTGTTCACCATTGAGCTGAAttt                                                                                      <  2:319748/251‑1 (MQ=255)
                                                                                                                                                          aTAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGgcg                                                                       <  1:1434/251‑1 (MQ=255)
                                                                                                                                                                   aCCACGGCACCAGGCAGTGACCCTGTGGTTTCTTCATCATCAGTAAGGATGTCGTTATTTGTTAGTGCCATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGCGTCTCAGGATATGATAATATCGTTCTGTCTCAACAGATCTACTAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTaa                                                               >  1:364205/1‑250 (MQ=255)
                                                                                                                                                                                      aCCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTaa                                            >  1:376080/1‑250 (MQ=255)
                                                                                                                                                                                        cATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGt                                         >  1:152495/1‑251 (MQ=255)
                                                                                                                                                                                                                  gCTGTCGTTATTTGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTc               >  1:315101/1‑251 (MQ=255)
                                                                                                                                                                                                                              tGTTAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGcc    <  2:329316/250‑1 (MQ=255)
                                                                                                                                                                                                                                ttAGTGCAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTTATAATATCGTTCTTTCTCAACAGATCTAATAATTTCTGTACCAAATTTGTTACTAATTTTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCTCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCtt  <  2:340204/250‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
GCGTAGCTTTCCGTCCATCAATGAATGCAGCGCAATCCCCACCACTACCGCGACGGGCATCAGCGCCAGAAAGAAAGGCCAGGTGTAGATAAAGAAGAACAGCGTGTTAGAGCCATAAATCAACATCGGCATCGCCAGCGCAAATAACCAGGAGATAAAACCGACCACGGCACCAGGCAGTGACCATGTGGTTTCTTCATCCTCAGTAAGGCTGTCGTTATTTGTTAGTGTAATGTTATGGCTATTACGCATATTTGATCCTGTTACTTTGACGAACCGGGCATGGAAACCCGGTGGTGTCTCAGGATCTGATAATATCGTTCTGTCTCAACAGATCTAATAATTGCTGTACCAAATTTGTTACTAATTGTTCACCATTGAGATGAATTTCTGCCGATTCAGGCGCTTCGTAAACGGAATCTATTCCCGTAAAGTTGCGCAGTTCACCGGCACGCGCTTTCTTATATAAGCCTT  >  NC_000913/2873086‑2873559

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