Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,137,543 G→A intergenic (‑300/‑359) wza ← / → yegH outer membrane polysaccharide export protein Wza/inner membrane protein YegH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,137,5430GA93.8% 43.9 / ‑4.5 16intergenic (‑300/‑359)wza/yegHouter membrane polysaccharide export protein Wza/inner membrane protein YegH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (7/8);  minor base C (1/0);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

GTCGCCATGCACATTTACCTTGCAGTTAATTGAATAAAAATTTAACTGGCATCAGTCCTAAAAAAATTGATTTCATCCGCAGGCTATTGACAGAATAATTCAGACTGGTCTTTCAGGCATCCAGACACGCTACCGCCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATA  >  NC_000913/2137297‑2137778
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
gTCGCCATGCACATTTACCTTGCAGTTAATTGAATAAAAATTTAACTGGCATCAGTCCTAAAAAAATTGATTTCATCCGCAGGCTATTGACAGCATAATTCAGACTGGTCTTTCAGGCATCCAGACACGCTACCGCCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAAGTg                                                                                                                                                                                                                                          >  2:215367/1‑250 (MQ=255)
 tCGCCATGCCCATTTACCTTGCAGTTAATTCCATCAAAATTTAACTGGCATCAGTCCTAAAACAATTGACTTCATCCGCAGGCTATTGAAAGACTAATTCAGCCTGGTCTTTCAGGCATCCAGACACGCTCCCGCACCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAACGTGa                                                                                                                                                                                                                                         >  2:167508/1‑250 (MQ=255)
                                   aaaaaTTTTACTGGTATCGGTCCTCAAAACATTGATTTAATCCGCAGGCTATTGACAGAATAATTAAGACTGGTATTTCAAGCATCCAGACACGCTACCGCCCCTGGCTTTTCAGCTACCAATACACTGCTTTAGTTTAATTTTTAACCCCCTATCAGCATGCAGTCGTTGATGCGGAAGGGTTATTACGGAAATTACCTTCCGAATATAAAGAGCCATTATGGTAATTGAATATTGGCTTTCCAATAATg                                                                                                                                                                                                      >  2:381591/1‑251 (MQ=255)
                                                       tCCTAAAAAAATTGATTTCATCCGCAGGCTATTGACAGAATAATTCAGACTGGTCTTTCAGGCATCCAGACACGCTACCGCCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAACGGGTTATTACGGAACTTAACTTCCGAATATAAAGTGACATTATGGTAATTGAATATTGGCTTTCCAATCATGCAGGAGGAAGTGTTACAGCt                                                                                                                                                                                  >  2:67753/1‑251 (MQ=255)
                                                                          aTCCGCAGGCTATTGACAGAATAATTCAGACTGGTCTTTCAGGCATCCAGACACGCTACCGCCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAAGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCa                                                                                                                                                                  >  2:392484/1‑248 (MQ=255)
                                                                                           aGAAAAATTCATGCTGTTCTTTCAGGCCTCCAGACACGCTAACGCCTCTGGCTTTTTAGCTAACAATACACTGATTTAGTTTAATTTTTCCCACACTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAAGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCATGCAAGATTACAATTCGGTAATTg                                                                                                                                              <  2:284661/251‑1 (MQ=255)
                                                                                                                       accAGGCACGCTACCGCCCCCGGCTTTTTAGCTTCCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCCGCATGCCGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAAGTGACATTATGGTAATTTAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTc                                                                                                                  <  2:151342/250‑1 (MQ=255)
                                                                                                                          aGACACGCTACCGCCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAAGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCAtt                                                                                                               <  1:253646/251‑1 (MQ=255)
                                                                                                                                           tGGCTTTATAGCTACCACTACACTGATTTGGTTTAATTTTTCACACCCTCTAAGCATGCAGTCGTTGATGAGAAAGGGTTATTAAGGAAATTAACTTCGGAATATAAAGAGACATTATGGTAATTGAATATTGGATTTACAATAATGAAAGAAGAAGTGTTACAGCTAACGGAATAGCAGGCAAAATAAAAATTCAGAAAATGGCTATTTTTAAGAATTATATTAAGAGTCATTCTATATGGATTTTAgg                                                                                               >  2:403600/1‑250 (MQ=255)
                                                                                                                                                     gCTACCAATACACTGATTTAGTTTAATTTTTCAAACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAAGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCttt                                                                                      >  1:128758/1‑249 (MQ=255)
                                                                                                                                                                              aTTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAAGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGtct                                                           <  1:136133/251‑1 (MQ=255)
                                                                                                                                                                                              cAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAAGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATAttt                                            <  1:299961/250‑1 (MQ=255)
                                                                                                                                                                                                          gTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAAGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGtt                               <  1:77465/251‑1 (MQ=255)
                                                                                                                                                                                                                             taCGGAAATTAACTTCCGAATATAAAGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTa             <  1:392484/250‑1 (MQ=255)
                                                                                                                                                                                                                             taCGGAAATTAACTTACGAATATAAAGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAATGTTACAGCTCACGGAATAGAAGGCAAGCTAACAATTCGGAAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCACTATGGTTTTTAGGAGTTTCTTTAGGTTGACACTATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTat            >  1:268724/1‑250 (MQ=255)
                                                                                                                                                                                                                                       aaCTTCCGAATATAAAGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATa  <  1:135506/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                           
GTCGCCATGCACATTTACCTTGCAGTTAATTGAATAAAAATTTAACTGGCATCAGTCCTAAAAAAATTGATTTCATCCGCAGGCTATTGACAGAATAATTCAGACTGGTCTTTCAGGCATCCAGACACGCTACCGCCCCTGGCTTTTTAGCTACCAATACACTGATTTAGTTTAATTTTTCACACCCTCTCAGCATGCAGTCGTTGATGAGAAAGGGTTATTACGGAAATTAACTTCCGAATATAAGGTGACATTATGGTAATTGAATATTGGCTTTCCAATAATGCAGGAGGAAGTGTTACAGCTAACGGAATAGCAGGCAAGATAACAATTCGGTAATTGGCTATTTTTAAGAATTATATTAAGTGTCATTCAATATGGTTTTTAGGAGTTTCTTTAGGTTGACAATATTTAATATAGTGTCTCCACATGCGATATTTCTTAAATAATGTTTTATTATTACCACTTTTAATTCAGGGATA  >  NC_000913/2137297‑2137778

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