Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,815,914 T→C 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,9140TC91.8% 193.2 / 6.0 61N702D (AAT→GAT) spoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (5/0);  new base C (50/6);  total (55/6)
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.00e-01

CTGCCGCCTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGATG  >  W3110S.gb/3815852‑3815980
                                                              |                                                                  
cTGCCGCCTGTCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtctcggtttc                                                            <  1:4061230/71‑1 (MQ=255)
           cGAGCATTTCGCAGATGCGTTCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTa                                                 <  1:3868376/71‑1 (MQ=255)
           cGAGCATTTCGCAGATGCGTGCATAACTTGTTTGGTTCATAAAACATTAAtctcggtttc                                                            <  1:3820479/60‑1 (MQ=255)
            gAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtctcgg                                                                >  1:372829/1‑55 (MQ=255)
            gAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTaa                                                >  1:3754703/1‑71 (MQ=255)
            gAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTaa                                                >  1:3123955/1‑71 (MQ=255)
            gAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTaa                                                >  1:3128660/1‑71 (MQ=255)
            gAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTaa                                                >  1:3003394/1‑71 (MQ=255)
            gAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTaa                                                >  1:3973651/1‑71 (MQ=255)
            gAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTaa                                                >  1:3351946/1‑71 (MQ=255)
            gAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTaa                                                >  1:1897866/1‑71 (MQ=255)
            gAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTaa                                                >  1:4124085/1‑71 (MQ=255)
            gAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGAGACTTTaa                                                >  1:3628491/1‑71 (MQ=255)
            gAGCATTTCGCAGATGCGTGCATAACGTCTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACttt                                                  >  1:511658/1‑69 (MQ=255)
            gAGCATTTCGCAGATGCCTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTaa                                                >  1:2323926/1‑71 (MQ=255)
            gAGCATATCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTaa                                                >  1:532602/1‑71 (MQ=255)
                   tCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAAt                                               >  1:3681182/1‑65 (MQ=255)
                   tCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAAt                                               >  1:2163113/1‑65 (MQ=255)
                   tCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAAt                                               >  1:2264653/1‑65 (MQ=255)
                   tCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAAt                                               >  1:2443780/1‑65 (MQ=255)
                   tCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTc                                         >  1:3404150/1‑71 (MQ=255)
                   tCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTc                                         >  1:1106106/1‑71 (MQ=255)
                   tCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTc                                         >  1:3215260/1‑71 (MQ=255)
                   tCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTc                                         >  1:529011/1‑71 (MQ=255)
                           gCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCa                                 >  1:1110455/1‑71 (MQ=255)
                                cATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCAcgcg                             >  1:2575164/1‑70 (MQ=255)
                                cATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCAcgcg                             >  1:2981872/1‑70 (MQ=255)
                                    aCGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGAttt                        >  1:1600211/1‑71 (MQ=255)
                                           gggTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCAt                 >  1:2504341/1‑71 (MQ=255)
                                           gggTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCAt                 >  1:2448585/1‑71 (MQ=255)
                                           gggTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCAt                 >  1:2192926/1‑71 (MQ=255)
                                           gggTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCAt                 >  1:4098697/1‑71 (MQ=255)
                                           gggTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCAt                 >  1:1182438/1‑71 (MQ=255)
                                           gggTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCAt                 >  1:100063/1‑71 (MQ=255)
                                              ttCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGa               >  1:3038239/1‑70 (MQ=255)
                                              ttCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGa               >  1:2525800/1‑70 (MQ=255)
                                                cATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGatat            >  1:3905850/1‑71 (MQ=255)
                                                cATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGatat            >  1:4001100/1‑71 (MQ=255)
                                                cATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGatat            >  1:966005/1‑71 (MQ=255)
                                                cATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGatat            >  1:477846/1‑71 (MQ=255)
                                                cATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGatat            >  1:3746692/1‑71 (MQ=255)
                                                cATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGatat            >  1:350449/1‑71 (MQ=255)
                                                cATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGatat            >  1:1804120/1‑71 (MQ=255)
                                                cATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGata             <  1:3013180/70‑1 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:3733108/1‑69 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:88317/1‑69 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:673269/1‑69 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:542652/1‑69 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:1270430/1‑69 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:1369759/1‑69 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:497077/1‑69 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:442843/1‑69 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:2317355/1‑69 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:2602562/1‑69 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:2782988/1‑69 (MQ=255)
                                                      aCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGAGCATGATATTCGc        >  1:2070257/1‑69 (MQ=255)
                                                         ttAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGAt   >  1:426153/1‑71 (MQ=255)
                                                         ttAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGAt   >  1:645721/1‑71 (MQ=255)
                                                         ttAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGAt   >  1:1016832/1‑71 (MQ=255)
                                                          tAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGATg  >  1:1627744/1‑71 (MQ=255)
                                                           aaTCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGATg  <  1:1727826/70‑1 (MQ=255)
                                                           aaTCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGATg  <  1:510592/70‑1 (MQ=255)
                                                           aaTCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGATg  <  1:3658496/70‑1 (MQ=255)
                                                           aaTCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGATg  <  1:3789020/70‑1 (MQ=255)
                                                              |                                                                  
CTGCCGCCTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGATG  >  W3110S.gb/3815852‑3815980

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