Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,371,135 G→T 29.3% intergenic (‑356/+39) sspA ← / ← rpsI stringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog/30S ribosomal subunit protein S9

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,371,1350GT29.3% 17.2 / 6.2 17intergenic (‑356/+39)sspA/rpsIstringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog/30S ribosomal subunit protein S9
Reads supporting (aligned to +/- strand):  ref base G (4/8);  new base T (3/2);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTT  >  CP009273/3370994‑3371277
                                                                                                                                             |                                                                                                                                              
gcggttgTCTGAAATAGTTGACTAAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCGGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTc                                                                                                                                         <  1:117992/143‑1 (MQ=255)
    tAGGCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGaaaa                                                                                                                                    <  2:200896/149‑1 (MQ=255)
       tCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAAtt                                                                                                                                  >  2:14647/1‑149 (MQ=255)
                   gACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATGGTTTTCTGCCgg                                                                                                                      >  1:133294/1‑149 (MQ=255)
                       tAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGCCGGAGCa                                                                                                                  <  1:102182/149‑1 (MQ=255)
                            ccTTATTGGCGATGTGGTTTTTGTTTTACTTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGc                                                                                                             <  2:158079/149‑1 (MQ=255)
                                aTTGGCGATGTGGTGTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAGGCAACTTACCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGTTTTTTTCGATAATTGATTTCTGCCGGAGCAGAGGCCAAt                                                                                                         >  2:294231/1‑149 (MQ=255)
                                                  gTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGc                                                                                       <  1:14647/149‑1 (MQ=255)
                                                               gtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTg                                                                          <  2:261934/149‑1 (MQ=255)
                                                                tcAGGTGGCCGCAAAAAGCAACTTTCCAGTTTTTACGTTGATTCAGATTTTAGCTGTAAAAATACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGTCGGAGCCGAAGCCAAGTACCGCTTGGAGAACTGTGGGCGACAAGGTGc                                                                         >  2:38982/1‑149 (MQ=255)
                                                                 cAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGAATTTAGCTATAAAAAAACCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGGCGACGACGGGCt                                                                        >  2:14108/1‑149 (MQ=255)
                                                                                        tCCAGTTTTTACGCTGGTTCAGATTTTAGCTATAAAATAACCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGCACTGCGGGCGACGACGTGCTTTACGCAGACCGACTTTCTTACg                                                 <  2:246253/149‑1 (MQ=255)
                                                                                                              aTTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTGTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGa                           <  2:9462/149‑1 (MQ=255)
                                                                                                                      tataAAAAAACCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCGATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTAAGTTCAACCTGACGCGCGGAGCGGGTAACGag                   >  1:385327/1‑148 (MQ=255)
                                                                                                                          aaaaaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCa               >  1:376714/1‑149 (MQ=255)
                                                                                                                                  cgcggAAGCGGGTTTTTTCGAAAACTGTTGTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCGGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACg       <  2:14560/149‑1 (MQ=255)
                                                                                                                                       aaGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGtt  <  1:853/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTT  >  CP009273/3370994‑3371277

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