Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,815,914 T→C 63.2% N702D (AAT→GAT)  spoT ← bifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,815,9140TC63.2% 12.5 / 16.2 19N702D (AAT→GAT) spoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base C (4/8);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

CCGCCTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAG  >  W3110S.gb/3815855‑3815977
                                                           |                                                               
ccgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtttcggtttcgg                                                       <  1:180985/70‑1 (MQ=255)
ccgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtttcggtttcgg                                                       <  1:903401/70‑1 (MQ=255)
ccgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtttcggtttcgg                                                       <  1:1615445/70‑1 (MQ=255)
ccgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtttcggtttcg                                                        <  1:1126703/69‑1 (MQ=255)
ccgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtttcggtttcg                                                        <  1:17358/69‑1 (MQ=255)
ccgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtttcggtttcg                                                        <  1:1773877/69‑1 (MQ=255)
ccgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtttcggtttcg                                                        <  1:281992/69‑1 (MQ=255)
 cgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGt                                                     <  1:162804/71‑1 (MQ=255)
 cgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGt                                                     <  1:21144/71‑1 (MQ=255)
                tCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTc                                      >  1:1050895/1‑71 (MQ=255)
                tCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTc                                      >  1:1765280/1‑71 (MQ=255)
                      aTGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                <  1:736613/71‑1 (MQ=255)
                      aTGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                <  1:1038810/71‑1 (MQ=255)
                      aTGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                <  1:1880041/71‑1 (MQ=255)
                      aTGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                <  1:862092/71‑1 (MQ=255)
                      aTGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                <  1:113591/71‑1 (MQ=255)
                      aTGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                <  1:1609309/71‑1 (MQ=255)
                       tGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                <  1:885854/70‑1 (MQ=255)
                                 aCGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGAttt                     >  1:1803435/1‑71 (MQ=255)
                                 aCGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGAttt                     >  1:3554/1‑71 (MQ=255)
                                 aCGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGAttt                     >  1:561589/1‑71 (MQ=255)
                                 aCGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGAttt                     >  1:844045/1‑71 (MQ=255)
                                 aCGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGAttt                     >  1:1659928/1‑71 (MQ=255)
                                          gTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGa            <  1:299288/71‑1 (MQ=255)
                                          gTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGa            <  1:556954/71‑1 (MQ=255)
                                                    cATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAg  <  1:1684920/71‑1 (MQ=255)
                                                           |                                                               
CCGCCTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAG  >  W3110S.gb/3815855‑3815977

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