Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,815,914 T→C 76.7% 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,9140TC76.7% 146.2 / 40.2 73N702D (AAT→GAT) spoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base C (26/30);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCGCCTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGATG  >  W3110S.gb/3815852‑3815980
                                                              |                                                                  
cTGCCGCCTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtctcggtttc                                                            <  1:1572512/71‑1 (MQ=255)
 tGCCGCCTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtctcggtttc                                                            <  1:1734950/70‑1 (MQ=255)
   ccgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtctcggtttcgg                                                          >  1:2767988/1‑70 (MQ=255)
   ccgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtctcggtttcgg                                                          >  1:2884700/1‑70 (MQ=255)
   ccgccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCggg                                                         >  1:458459/1‑71 (MQ=255)
     gccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTg                                                       <  1:2574415/71‑1 (MQ=255)
     gccTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTg                                                       <  1:1370031/71‑1 (MQ=255)
        tGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACt                                                    >  1:3317009/1‑71 (MQ=255)
          gCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACt                                                    <  1:2663889/69‑1 (MQ=255)
          gCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACt                                                    <  1:889340/69‑1 (MQ=255)
          gCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACt                                                    <  1:3306332/69‑1 (MQ=255)
           cGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTa                                                 <  1:1779859/71‑1 (MQ=255)
           cGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTa                                                 <  1:2272581/71‑1 (MQ=255)
           cGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTa                                                 <  1:3629738/71‑1 (MQ=255)
                aTTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAAtctcggttt                                                             >  1:2089225/1‑54 (MQ=255)
                    cGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCt                                        >  1:3616709/1‑71 (MQ=255)
                    cGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCt                                        >  1:3544331/1‑71 (MQ=255)
                    cGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCt                                        >  1:3548431/1‑71 (MQ=255)
                      cAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTgg                                      <  1:1284014/71‑1 (MQ=255)
                         aTGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                   <  1:3350680/71‑1 (MQ=255)
                         aTGCGTGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                   <  1:37185/71‑1 (MQ=255)
                          tGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                   >  1:303687/1‑70 (MQ=255)
                          tGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                   >  1:2724300/1‑70 (MQ=255)
                          tGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                   >  1:220910/1‑70 (MQ=255)
                          tGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCAt                                   >  1:3460702/1‑70 (MQ=255)
                              tGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCTGGTGACTTTa                                                 <  1:1181870/52‑1 (MQ=255)
                              tGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTa                                                 <  1:1363704/52‑1 (MQ=255)
                              tGCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTa                                                 <  1:567273/52‑1 (MQ=255)
                               gCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTc                                         >  1:3045131/1‑59 (MQ=255)
                               gCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTc                                         >  1:989367/1‑59 (MQ=255)
                               gCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTc                                         >  1:915509/1‑59 (MQ=255)
                               gCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTc                                         >  1:142330/1‑59 (MQ=255)
                               gCATAACGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTc                                         >  1:2205368/1‑59 (MQ=255)
                                    aCGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGAttt                        >  1:2921837/1‑71 (MQ=255)
                                     cGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGAtttt                       >  1:2360850/1‑71 (MQ=255)
                                     cGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGAtttt                       >  1:2762620/1‑71 (MQ=255)
                                     cGTGTTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGAtttt                       >  1:2954821/1‑71 (MQ=255)
                                      gtgtTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTg                      >  1:2916855/1‑71 (MQ=255)
                                      gtgtTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTg                      >  1:2010084/1‑71 (MQ=255)
                                      gtgtTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTg                      >  1:3172576/1‑71 (MQ=255)
                                      gtgtTGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTg                      >  1:204509/1‑71 (MQ=255)
                                         ttGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTgcg                    >  1:669584/1‑70 (MQ=255)
                                         ttGGGTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTgcg                    <  1:1425887/70‑1 (MQ=255)
                                            ggTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATg                <  1:235959/71‑1 (MQ=255)
                                            ggTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATg                <  1:1021555/71‑1 (MQ=255)
                                            ggTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATg                <  1:3306625/71‑1 (MQ=255)
                                            ggTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATg                <  1:2324595/71‑1 (MQ=255)
                                             gTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATg                <  1:2039690/70‑1 (MQ=255)
                                             gTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATg                <  1:2989605/70‑1 (MQ=255)
                                             gTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGa               <  1:1642073/71‑1 (MQ=255)
                                             gTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGa               <  1:565790/71‑1 (MQ=255)
                                             gTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGa               <  1:2956671/71‑1 (MQ=255)
                                             gTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGa               <  1:3227332/71‑1 (MQ=255)
                                             gTTCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGa               <  1:3042441/71‑1 (MQ=255)
                                              ttCATAAAACATTAATTTCGGTTTCGGGTGACttt                                                  <  1:3098254/35‑1 (MQ=255)
                                              ttCATAAAACATTAATTTCGGTTTCGGGTGACttt                                                  <  1:1883635/35‑1 (MQ=255)
                                               tCATAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGata             >  1:1879884/1‑71 (MQ=255)
                                                cATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGata             <  1:1789641/70‑1 (MQ=255)
                                                 aTAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGGCTGGCATCACGCGGATTTTGCGCATGata             <  1:3076742/69‑1 (MQ=255)
                                                 aTAAAACATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGata             <  1:820132/69‑1 (MQ=255)
                                                  tAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGata             <  1:3006503/68‑1 (MQ=255)
                                                  tAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGata             <  1:1455925/68‑1 (MQ=255)
                                                  tAAAACATTAATCTCGGTTTCGGGTGACTTTAATc                                              <  1:3018853/35‑1 (MQ=255)
                                                  tAAAACATTAATCTCGGTTTCGGGTGACTTTAATc                                              <  1:3689060/35‑1 (MQ=255)
                                                  tAAAACATTAATCTCGGTTTCGGGTGACTTTAATc                                              <  1:2118174/35‑1 (MQ=255)
                                                   aaaaCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCg         >  1:2237456/1‑71 (MQ=255)
                                                    aaaCATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGc        >  1:314901/1‑71 (MQ=255)
                                                    aaaCATTAATCTCGGTTTCGGGTGACTTTAATCAc                                            >  1:3352393/1‑35 (MQ=255)
                                                    aaaCATTAATCTCGGTTTCGGGTGACTTTAATCAc                                            >  1:1309708/1‑35 (MQ=255)
                                                    aaaCATTAATCTCGGTTTCGGGTGACTTTAATCAc                                            >  1:2307731/1‑35 (MQ=255)
                                                    aaaCATTAATCTCGGTTTCGGGTGACTTTAATCAc                                            >  1:910156/1‑35 (MQ=255)
                                                    aaaCATTAATCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATGTTCGc        >  1:1303563/1‑71 (MQ=255)
                                                      aCATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCa      <  1:1797279/71‑1 (MQ=255)
                                                      aCATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCa      <  1:2275700/71‑1 (MQ=255)
                                                         ttAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGAt   >  1:3607133/1‑71 (MQ=255)
                                                         ttAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGAt   >  1:1433647/1‑71 (MQ=255)
                                                         ttAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGAt   >  1:60260/1‑71 (MQ=255)
                                                           aaTCTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGATg  >  1:1061853/1‑70 (MQ=255)
                                                              |                                                                  
CTGCCGCCTGGCGAGCATTTCGCAGATGCGTGCATAACGTGTTGGGTTCATAAAACATTAATTTCGGTTTCGGGTGACTTTAATCACGTCTGGCATCACGCGGATTTTGCGCATGATATTCGCCAGATG  >  W3110S.gb/3815852‑3815980

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