Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,101,857 G→A intergenic (‑113/+38) wbbL ← / ← insH1 pseudogene, lipopolysaccharide biosynthesis protein/IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,101,8570GA88.0% 53.4 / 1.8 25intergenic (‑113/+38)wbbL/insH1pseudogene, lipopolysaccharide biosynthesis protein/IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (10/12);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

TTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTT  >  NC_000913/2101722‑2101883
                                                                                                                                       |                          
tttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAAcc                         <  1:272552/139‑1 (MQ=255)
 ttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCg                        >  2:327677/1‑139 (MQ=255)
 ttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCg                        >  1:1100185/1‑139 (MQ=255)
 ttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCg                        >  2:175157/1‑139 (MQ=255)
 ttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCg                        >  1:239854/1‑139 (MQ=255)
 ttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCg                        >  2:1372861/1‑139 (MQ=255)
 ttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCg                        >  1:484720/1‑139 (MQ=255)
 ttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGAGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCg                        >  1:1042095/1‑139 (MQ=255)
  tctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGt                       >  1:461626/1‑139 (MQ=255)
         gAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCg                <  1:493353/139‑1 (MQ=255)
         gAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCg                <  2:1164278/139‑1 (MQ=255)
         gAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCg                <  1:251742/139‑1 (MQ=255)
              tttGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCGCCAtt           >  1:1295903/1‑139 (MQ=255)
               ttGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCGCCAttt          <  1:498441/139‑1 (MQ=37)
                tGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCGCCATTTa         <  2:458296/139‑1 (MQ=37)
                  ttttATTAAGGAAGGTGTGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCGCCATTTAAg       <  2:1460679/139‑1 (MQ=21)
                  ttttATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCGCCATTTAAg       <  2:461626/139‑1 (MQ=21)
                    ttattaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCGCCATTTAAGGc     >  1:685726/1‑139 (MQ=18)
                     tattaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCGCCATTTAAGGCg    <  1:902694/139‑1 (MQ=17)
                     tattaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCGCCATTTAAGGCg    <  2:239854/139‑1 (MQ=17)
                     atttaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg    >  2:889775/3‑139 (MQ=11)
                     atttaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCg    >  2:1043692/3‑139 (MQ=11)
                      attaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCGCCATTTAAGGCGt   <  2:1042095/139‑1 (MQ=14)
                       ttaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCTTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGtt  <  1:1043692/139‑1 (MQ=11)
                       ttaAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAAACCGTTTCTTCGCCATTTAAGGCGtt  <  2:685726/139‑1 (MQ=11)
                                                                                                                                       |                          
TTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTT  >  NC_000913/2101722‑2101883

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