Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,111,727 2 bp→TT coding (865‑866/2544 nt) opgH → osmoregulated periplasmic glucans biosynthesis protein H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,111,7270CT97.2% 116.3 / ‑3.9 37P289M (CCG→TTG) opgHosmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (23/12);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01
*NC_0009131,111,7280CT97.2% 123.8 / ‑3.7 37P289M (CCG→TTG) opgHosmoregulated periplasmic glucans biosynthesis protein H
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (23/12);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

GACGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCG  >  NC_000913/1111492‑1111970
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                                  
gACGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCg                                                                                                                                                                                                                                      >  1:221545/1‑251 (MQ=255)
         aaTGGGCTTGCTGCAACTGCTTATTGGTAGCGATAAATACAGGATATCGGCGTCAACAGGTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCaaa                                                                                                                                                                                                                             <  2:9766/251‑1 (MQ=255)
             ggCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACTAAGACGTGAACCGTGTTTTTGCTAGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGc                                                                                                                                                                                                                         <  1:452265/251‑1 (MQ=255)
                  ccTGCAACTGCTTATTGGGCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGCCGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCtgg                                                                                                                                                                                                                     <  1:233542/250‑1 (MQ=255)
                     gCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCtggatgg                                                                                                                                                                                                                 >  1:28726/1‑251 (MQ=255)
                             ttATTGGTCGCTAGAAATACAGTATAGCTGCGTCAGCAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTTAACCGTGGTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACGTTGATGTCTACATTCTTAGTGACAGTTATAACTTTGAGATCTGCGTCGCAGATCAAAAAGCCTGGATGTAGCTTAt                                                                                                                                                                                                          <  2:525131/250‑1 (MQ=255)
                                            aaTACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAACGCCTGGATGGAGCTTATCGCTGAAGTCGGTggg                                                                                                                                                                                          >  2:29184/1‑250 (MQ=255)
                                               aCAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGCGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAACGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGAGATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTAACGCGGAAGTCGGTGGCgg                                                                                                                                                                                        >  2:371267/1‑249 (MQ=255)
                                                  gTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTCAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGt                                                                                                                                                                                    >  1:4964/1‑251 (MQ=255)
                                                      atCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTc                                                                                                                                                                                   >  1:159403/1‑248 (MQ=255)
                                                       tCTGCGTCAACAGTTGGCGATGAACCATTAAACCAGGAGCATCGCACGGCGTTGCTCATGCCTCTCTGTACCGACGACGTGAACCGTGATTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTATACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGGGGCGAAggac                                                                                                                                                                                   >  2:326148/1‑245 (MQ=255)
                                                       tCTGCGTCAACAGTTGGAGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGCCGACGTAAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTc                                                                                                                                                                                   >  2:254455/1‑247 (MQ=255)
                                                               aaCAGTTGGAGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATTGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATc                                                                                                                                                                       <  1:289429/251‑1 (MQ=255)
                                                                        cGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCgccgccgcc                                                                                                                                                              <  1:89352/251‑1 (MQ=255)
                                                                                       cccGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGc                                                                                                                                               >  2:188089/1‑251 (MQ=255)
                                                                                        ccGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCg                                                                                                                                              >  2:387554/1‑251 (MQ=255)
                                                                                                  gCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTTGCCGCGTGAAGCGTAAAAGCGGt                                                                                                                                    >  1:86975/1‑251 (MQ=255)
                                                                                                            gATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATg                                                                                                                          <  1:524481/251‑1 (MQ=255)
                                                                                                                aTGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACtt                                                                                                                      >  2:27427/1‑251 (MQ=255)
                                                                                                                         tGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGACTGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCg                                                                                                             >  2:325501/1‑251 (MQ=255)
                                                                                                                             aCGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACGTTGATGTCTACATTCTTAGTGACAGTTATAACGTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGGCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTgg                                                                                                         <  2:159403/251‑1 (MQ=255)
                                                                                                                                            gtgtTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTg                                                                                                                                                   >  1:42290/1‑194 (MQ=255)
                                                                                                                                                    gCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACAtggt                                                                                  >  2:127023/1‑251 (MQ=255)
                                                                                                                                                    gCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACAtggt                                                                                  <  1:27427/251‑1 (MQ=255)
                                                                                                                                                     cTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACAtggtg                                                                                 >  2:245956/1‑251 (MQ=255)
                                                                                                                                                          cTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGCATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGAATATCTGCGTCGCAGAGCa                                                                                                                                                                                                                               >  1:254107/1‑104 (MQ=255)
                                                                                                                                                                   aCGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTg                                                                    >  2:377146/1‑250 (MQ=255)
                                                                                                                                                                                    aGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTg                                                  >  1:487691/1‑251 (MQ=255)
                                                                                                                                                                                           gggAATTCCAAACCCTGTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGGAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGGCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGGCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGAttgttt                                           <  1:127865/251‑1 (MQ=255)
                                                                                                                                                                                                 gCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCgg                                     <  1:245956/251‑1 (MQ=255)
                                                                                                                                                                                                      acacTTTGATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTgg                                >  1:127183/1‑251 (MQ=255)
                                                                                                                                                                                                             gATGTCTACATTCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGcc                          <  1:321570/250‑1 (MQ=255)
                                                                                                                                                                                                                       ttCTTAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAACAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGcc               >  2:320036/1‑251 (MQ=255)
                                                                                                                                                                                                                           tAGTGACAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCaa             <  1:371267/249‑1 (MQ=255)
                                                                                                                                                                                                                                aCAGTTATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGCAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGGGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAAc      >  2:413253/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gTTATAACTGGGATAGCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCTAACGcc   <  2:221545/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ttATAACTTGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGACGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCg  >  1:309935/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ||                                                                                                                                                                                                                                                  
GACGGCGTTAATGGGCTTCCTGCAACTGCTTATTGGTCGCGATAAATACAGTATATCTGCGTCAACAGTTGGCGATGAACCATTAAACCCGGAGCATCGCACGGCGTTGATCATGCCTATCTGTAACGAAGACGTGAACCGTGTTTTTGCTGGCCTGCGTGCAACGTGGGAATCAGTAAAAGCCACCGGGAATGCCAAACACTTTGATGTCTACATTCTTAGTGACAGTTATAACCCGGATATCTGCGTCGCAGAGCAAAAAGCCTGGATGGAGCTTATCGCTGAAGTCGGTGGCGAAGGTCAGATTTTCTATCGCCGCCGCCGTCGCCGCGTGAAGCGTAAAAGCGGTAATATCGATGACTTCTGCCGTCGCTGGGGCAGCCAGTACAGCTACATGGTGGTGCTGGATGCTGACTCGGTAATGACCGGTGATTGTTTGTGCGGGCTGGTGCGCCTGATGGAAGCCAACCCGAACGCCG  >  NC_000913/1111492‑1111970

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: ‑