Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,370,206 C→A 17.8% E192* (GAG→TAG)  sspA ← stringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,370,2060CA17.8% 63.7 / 10.0 34E192* (GAG→TAG) sspAstringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog
Reads supporting (aligned to +/- strand):  ref base C (11/17);  new base A (3/3);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GTTATCCAGCAACCACTCATAGAATGCACGCAGCAGATAGGGACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCACAGCAGCGGAGCAAGATAGCAATCGACCAGGCTGAACTCATCGCTCAGGAAGTACGGCTTCTGACCGA  >  CP009273/3370064‑3370352
                                                                                                                                              |                                                                                                                                                  
gTTATCCAGCAACCACTCATAGAGTGCACGCAGCAGATAGGGACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTAAAAGAc                                                                                                                                              <  2:254025/149‑1 (MQ=255)
gTTATCCAGCAACCACTCATAGAATGCACGCAGCAGATAGGGACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTc                                                                                                                                                    >  1:172838/1‑143 (MQ=255)
     ccAGCAACCACTCATAGAATGCACGCAGCAGATAGGGACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGa                                                                                                                                               <  1:171135/143‑1 (MQ=255)
      cAGCAACCACTCATAGAATGCACGCAGCAGATAGGGCCGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAACGAAAGAGTCACGCCCAAAGAc                                                                                                                                              >  1:64202/1‑143 (MQ=255)
      cAGCAACCACTCATAGAATGCACGCAGCAGATAGGGACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGt                                                                                                                                        >  2:189213/1‑149 (MQ=255)
         cAACCACTCATAGAATGCACGCAGCAGATAGGGACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCCGGTAAAGAAGCAAGGGAAGAGTCCCGCTCAAAGACGCGGGTCat                                                                                                                                     >  2:118947/1‑149 (MQ=255)
         cAACCACTCATAGAATGCACGCAGCAGATAGGGACGACGTGATGTTAGCTGTGACAAATCCATACGGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGAcgcg                                                                                                                                           >  1:231803/1‑143 (MQ=255)
                  aTAGAATGCACGCAGCAGATAGGGACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCttt                                                                                                                            >  2:348424/1‑149 (MQ=255)
                         gCACGCAGCAGATAGGGACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTAAAAGACGCGGGTCATATAGCCTTTc                                                                                                                           <  1:15094/143‑1 (MQ=255)
                                  agATAGGGACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTc                                                                                                                  <  1:61306/143‑1 (MQ=255)
                                    aTAGGGACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGc                                                                                                                <  1:495636/143‑1 (MQ=255)
                                        ggACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTAAAAGACGCGGGTCATATAGCCTTTCAGCCCTTTCGCAccc                                                                                                            >  1:341176/1‑143 (MQ=255)
                                                  gtgtTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGa                                                                                            <  2:511641/149‑1 (MQ=255)
                                                        gCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCATCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTAATAAAGCCTTTCAGCTCTTTCGCACCCCGGCCGCTGAACTCGa                                                                                            >  1:199883/1‑143 (MQ=255)
                                                         cTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAGGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCCGAAATCGag                                                                                           >  1:116178/1‑142 (MQ=255)
                                                            tgACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCATCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCCGCTCTTTCGCACCACGGCCGCTGAACGCCATGcc                                                                                        >  1:364566/1‑143 (MQ=255)
                                                             gACAAATCCATACAGAATAACGCCGGCCCAGACGCATTGCACGTTCTGCTGCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGc                                                                                 <  2:263400/149‑1 (MQ=255)
                                                                   tCCATACAGATTAAGTCCGGCCCAGACGCATTTCACATTATGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCTCAGTTGCGGCAGa                                                                           <  2:451245/149‑1 (MQ=255)
                                                                   tCCATACAGATTAACTCCGGCCCAGACGCATTTCCCGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGa                                                                           <  2:92842/149‑1 (MQ=255)
                                                                   tCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGa                                                                           <  2:92834/149‑1 (MQ=255)
                                                                    ccgtACCGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTAAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCg                                                                                <  1:548185/140‑1 (MQ=255)
                                                                     cATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCgg                                                                               <  1:270144/143‑1 (MQ=255)
                                                                      aTACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTAAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGc                                                                              >  1:326620/1‑143 (MQ=255)
                                                                       tACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGGCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCa                                                                             >  1:218955/1‑143 (MQ=255)
                                                                                  tCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCAcagc                                                                  >  1:512422/1‑143 (MQ=255)
                                                                                  tCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCACAGCAGCGGa                                                            <  2:475822/149‑1 (MQ=255)
                                                                                  tCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTAAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCACAACAGCGGa                                                            >  2:260650/1‑149 (MQ=255)
                                                                                                 tttCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCACAGCAGCGGAGCAAGATAg                                                   <  1:153407/143‑1 (MQ=255)
                                                                                                       gTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCACAGCAGCGGAGCAAGATAGCAATCGACCAgg                                       <  2:262335/149‑1 (MQ=255)
                                                                                                              ttCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCACAGCAGCGGAGCAAGATAGCAATCGACCAGGc                                      <  1:58255/143‑1 (MQ=255)
                                                                                                                   ttAGAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCACAGCAGCGGAGCAAGATAGCAATCGACCAGGCTGAACTCATCg                           <  2:413932/149‑1 (MQ=255)
                                                                                                                   ttAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCACAGCAGCGGAGCAAGATAGCAATCGACCAGGCTGAACTCATCg                           <  2:414725/149‑1 (MQ=255)
                                                                                                                     aaaGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCCTTTCTCGCCGCAGCAGCGGAGCAAGATAGCAATCGACCAGGCTGAACTCaacccg                         >  2:18694/1‑144 (MQ=255)
                                                                                                                                      gTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCACAGAAGCGGATTAAGATAGCAATCGACCAGGCTGAACTCATCGCTCAGGAAGTACGGCTTCt        <  2:338158/149‑1 (MQ=255)
                                                                                                                                            cTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCACAGCAGCGGAGCAAGATAGCAATCGACCAGGCTGAACTCATCGCTCAGGAAGTACGGCTTCTGACCGa  <  2:229591/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTTATCCAGCAACCACTCATAGAATGCACGCAGCAGATAGGGACGACGTGGTGTTAGCTGTGACAAATCCATACAGATTAACTCCGGCCCAGACGCATTTCACGTTCTGCTTCAGTTAAAGAAGCAAGGAAAGAGTCACGCTCAAAGACGCGGGTCATATAGCCTTTCAGCTCTTTCGCACCCGGGCCGCTGAACTCGATGCCCAGTTGCGGCAGACGCCACAGCAGCGGAGCAAGATAGCAATCGACCAGGCTGAACTCATCGCTCAGGAAGTACGGCTTCTGACCGA  >  CP009273/3370064‑3370352

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: