Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,816,294 C→T 25.0% R575H (CGC→CAC)  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,816,2940CT25.0% 50.5 / 13.6 28R575H (CGC→CAC) spoTbifunctional (p)ppGpp synthetase II and guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base C (8/13);  new base T (3/4);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGGTGGATCACCAGACCTTTACCGGGGCTGACGTGGGCGATAATCGGGTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATGGGCAGATGTCCGTGGCTTTG  >  W3110S.gb/3816229‑3816361
                                                                 |                                                                   
cATGGTGGATCACCAGACCTTTACCGGGGCTGACGTGGGCGATAATCGGGTCGCCAGGAATAGGGCGgcag                                                                >  1:648215/1‑71 (MQ=255)
cATGGTGGATCACCAGACCTTTACCGGGGCTGACGTGGGCGATAATCGGGTCGCCAGGAATAGGGCGgca                                                                 >  1:875792/1‑70 (MQ=255)
    gtggATCACCAGACCTTTACCGGGGCTGACGTGGGCGATAATCGGGTCGCCAGGAATAGGGCGGCAGCAtt                                                            <  1:2215921/71‑1 (MQ=255)
               gACCTTTACCGGGGCTGACGTGGGCGATAATCGGGTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTg                                                 <  1:1713039/71‑1 (MQ=255)
                     tACCGGGGCTGACGTGGGGGATAATCGGGTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGAt                                               >  1:804427/1‑67 (MQ=255)
                     tACCGGGGCTGACGTGGGCGATAATCGGGTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGa                                                >  1:509236/1‑66 (MQ=255)
                     tACCGGGGCTGACGTGGGCGATAATCGGGTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAg                                            >  1:239675/1‑70 (MQ=255)
                           ggCTGACGTGGGCGATAATCGGGTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCa                                     <  1:1842825/71‑1 (MQ=255)
                                aCGTGGGCGATAATCGGGTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGc                                <  1:1437642/71‑1 (MQ=255)
                                aCGTGGGCGATAATCGGGTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGc                                <  1:1826893/71‑1 (MQ=255)
                                       cGATAATCGGGTCGCCAGGAATAGGGTGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTa                         <  1:1680468/71‑1 (MQ=255)
                                       cGATAATCGGGTCGCCAGGAATAGGGTGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTa                         <  1:1683540/71‑1 (MQ=255)
                                          tAATCGGGTCGCCAGGAATAGGGTGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTa                         <  1:2053641/68‑1 (MQ=255)
                                          tAATCGGGTCGCCAGGAATAGGGTGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTa                         <  1:444904/68‑1 (MQ=255)
                                            aTCGGGTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATgg                     >  1:2100975/1‑70 (MQ=255)
                                             tCGGGTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATggg                    >  1:25300/1‑70 (MQ=255)
                                              cGGGTCGCCAGGAATAGGGTGGCAGCATTTCGCAAAGGTGATCa                                             >  1:56912/1‑44 (MQ=255)
                                              cGGGTCGCCAGGAATAGGGTGGCAGCATTTCGCAAAGGTGATCa                                             >  1:1802586/1‑44 (MQ=255)
                                              cGGGTCGCCAGGAATAGGGTGGCAGCATTTCGCAAAGGTGATCa                                             >  1:2405534/1‑44 (MQ=255)
                                               gggTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATGGGCa                  <  1:54969/70‑1 (MQ=255)
                                                  tCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATGGGCAGATg              >  1:36915/1‑71 (MQ=255)
                                                    gCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATGGGCAGATGt             <  1:2153565/70‑1 (MQ=255)
                                                    gCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATGGGCAGATGt             <  1:1241682/70‑1 (MQ=255)
                                                    gCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATGGGCAGATGTc            <  1:2181577/71‑1 (MQ=255)
                                                    gCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATGGGCAGATGTc            <  1:604945/71‑1 (MQ=255)
                                                         gAATAGGGCGGCAGCATTTCGCAAACGTGATCAGCACGCCATCGGCACCTTTAATGGGCAg                 <  1:920355/61‑1 (MQ=255)
                                                              gggCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATGGGCAGATGTCCGTGGCttt   <  1:2035248/70‑1 (MQ=255)
                                                              gggCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATGGGCAGATGTCCGTGGCTTTg  <  1:1463794/71‑1 (MQ=255)
                                                                 |                                                                   
CATGGTGGATCACCAGACCTTTACCGGGGCTGACGTGGGCGATAATCGGGTCGCCAGGAATAGGGCGGCAGCATTTCGCAAAGGTGATCAGCACGCCATCGGCACCTTTAATGGGCAGATGTCCGTGGCTTTG  >  W3110S.gb/3816229‑3816361

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: ‑