Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,371,135 G→T 25.0% 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,1350GT25.0% 35.3 / 6.2 24intergenic (‑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 (11/7);  new base T (2/4);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACG  >  CP009273/3370997‑3371272
                                                                                                                                          |                                                                                                                                          
ttAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTATTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGtttttt                                                                                                                                     >  2:302983/1‑145 (MQ=255)
gtAGTCTGTAATAGTTGCCGTCCTCCCTTATTGGCTACGGTTCTTTTGTTTCTCTTCTCTGTCAGGCTGCGCCAAATAGCAACCTTTCTGTTTTTACTCTTATTCATATTTTTGCTATAAAAAAACCCGCCGAAGCGGGtttttt                                                                                                                                     <  1:58563/144‑1 (MQ=255)
    tCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAACAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGaa                                                                                                                                 >  2:204088/1‑145 (MQ=255)
      tGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGaaaa                                                                                                                               <  2:390789/145‑1 (MQ=255)
                  cTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTGTTACGCTGATTCAGGTGTTAGCTATAAAAAAACCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGcc                                                                                                                   >  1:345388/1‑145 (MQ=255)
                                      gtgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTACTGCCGGAGCAGAAGCCAATTAACg                                                                                               <  2:253608/145‑1 (MQ=255)
                                      ggggttttggttttACCTATCTGTCAGGTAGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCACATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACg                                                                                               <  1:71142/143‑1 (MQ=255)
                                       tgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGCGGCATAAGCCAAGTACCGt                                                                                              >  2:152980/1‑145 (MQ=255)
                                                            gccAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacga                                                                         <  2:265636/143‑1 (MQ=255)
                                                                 gTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGc                                                                    >  2:318793/1‑145 (MQ=255)
                                                                    gcagcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCAGCGGAGGGGGGTTTTTTCGAAAGTTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGGGCttt                                                                 >  1:234169/1‑145 (MQ=255)
                                                                       gcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGGGCTTTACg                                                              >  2:244656/1‑145 (MQ=255)
                                                                         aaaaaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCCGAAGCCAAATAACGCTTGGCGAACTACGgaagacgacatacgatacgcc                                                            >  2:240401/1‑133 (MQ=255)
                                                                           aaaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGa                                                          >  1:361040/1‑145 (MQ=255)
                                                                            aaGCAACTTTCCAGTTTGTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGAc                                                         <  1:163718/145‑1 (MQ=255)
                                                                              gCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCCGACCg                                                       >  2:298543/1‑145 (MQ=255)
                                                                                   tttCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGGGGACGACGACGTGCTTTACGCAGACCGACctt                                                  >  1:43776/1‑145 (MQ=255)
                                                                                     tCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGTTCTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCt                                                <  1:304457/145‑1 (MQ=255)
                                                                                        aGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGAAGAACTGCGGACGACGACGTGCTTTACGCAGACAGACTGTCTTAc                                             >  2:257123/1‑145 (MQ=255)
                                                                                         gTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACg                                            <  1:130931/145‑1 (MQ=255)
                                                                                               tacgctagtccaaTTTTAGGAATAAAAAAACCCGGCGAAGCGGTTTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCaa                                       <  2:263274/133‑1 (MQ=255)
                                                                                                                    ataAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAAc                 >  1:362262/1‑145 (MQ=255)
                                                                                                                      aaaaaaaCCCGCCGAAGCGGTTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGa               <  1:305366/145‑1 (MQ=255)
                                                                                                                                   gAAGCGGGGTTTTGCGAAAATGGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACg  <  2:83969/145‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
TTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACG  >  CP009273/3370997‑3371272

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