Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,989,071 C→T 10.5% G17S (GGT→AGT)  cspR ← TrmH family RNA methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,989,0710CT10.5% 110.8 / 5.0 38G17S (GGT→AGT) cspRTrmH family RNA methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (17/17);  new base T (2/2);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

ATACGTAATATTCACAAATTCCCAATAATCTAAACCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCCTCTTGATCATAATTAGCATTGATAAAATCTCT  >  CP000730/1988940‑1989210
                                                                                                                                   |                                                                                                                                           
aTACGTAATATTCACAAATTCCCAATAATCTAAACCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGc                                                                                                                                    >  1:56042/1‑141 (MQ=255)
 tACGTAATATTCACAAATTCCCAATAATCTAAACCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCt                                                                                                                                   >  2:345682/1‑141 (MQ=255)
   cGTAATATTCACAAATTCCCAATAATCTAAACCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTgg                                                                                                                                 <  2:313460/141‑1 (MQ=255)
      aaTATTCACAAATTCCCAATAATCTAAACCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACTAGTATTAGCTGGGAt                                                                                                                              >  1:360397/1‑141 (MQ=255)
      aaTATTCACAAATTCCCAATAATCTAAACCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACTAGTATTAGCTGGGAt                                                                                                                              >  1:360407/1‑141 (MQ=255)
         atTCACAAATTCCCAATAATCTAAACCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTc                                                                                                                           <  1:131758/141‑1 (MQ=255)
                 aTTCCCAATAATCTAAACCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGat                                                                                                                   <  1:369910/141‑1 (MQ=255)
                    cccAATAATCTAAACCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGAtata                                                                                                                <  2:326469/141‑1 (MQ=255)
                      cAATGATCTAACCCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATaaa                                                                                                              <  2:69137/141‑1 (MQ=255)
                       aataatCTAAACCGGCTCTTTTTAACATTTTATCAGCAGTTCGAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGCTTTCTGGTTGATATaaaa                                                                                                             >  1:324234/1‑141 (MQ=255)
                       aataatCTAAACCGGCTCTTTTTAACATTTTATCAACAGTTCGAAATCCAAGAGGTTTAATTAAATGCAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATaaaa                                                                                                             >  1:324256/1‑141 (MQ=255)
                             cTAGACCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGatat                                                                                                       <  2:2815/141‑1 (MQ=255)
                                    ggCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGt                                                                                                <  2:74454/141‑1 (MQ=255)
                                    ggCTCTTTGTAACATTTTATCATCAGTTATAAATCCAAGAGGATTAATTAAATGTAAATGTGTGTTTGTAACTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATGTCTGGTTGATATAAAACTATATGATTTGt                                                                                                >  2:18903/1‑141 (MQ=255)
                                       tctTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCat                                                                                             <  1:187061/141‑1 (MQ=255)
                                            ttAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTAc                                                                                        >  1:218254/1‑141 (MQ=255)
                                               aCATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTAt                                                                                     >  2:349436/1‑141 (MQ=255)
                                                  ttttATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATtct                                                                                  >  2:99953/1‑141 (MQ=255)
                                                        catcaGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACTAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCtt                                                                            <  2:360407/141‑1 (MQ=255)
                                                        catcaGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGCGTTTGTACCTGCACACGTACGTGCAATGTTACTAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCtt                                                                            <  2:360397/141‑1 (MQ=255)
                                                         atcaGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTACACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTtg                                                                           >  2:348692/1‑141 (MQ=255)
                                                                                      aatgtaaatgtGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGa                                              >  2:168944/1‑141 (MQ=255)
                                                                                       atgtaaatgtGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGaa                                             >  1:8672/1‑141 (MQ=255)
                                                                                            aatgtGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCt                                        >  1:368122/1‑141 (MQ=255)
                                                                                            aatgtGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCt                                        >  2:271625/1‑141 (MQ=255)
                                                                                               gtgtgtTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCa                                     <  1:190071/141‑1 (MQ=255)
                                                                                                  tgtTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATcc                                  <  2:95937/141‑1 (MQ=255)
                                                                                                  tgtTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATcc                                  >  2:193770/1‑141 (MQ=255)
                                                                                                   gtTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCct                                 <  2:73282/141‑1 (MQ=255)
                                                                                                   gtTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCct                                 <  1:96473/141‑1 (MQ=255)
                                                                                                      tGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCCTCtt                              <  1:158555/141‑1 (MQ=255)
                                                                                                       gTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATCACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCCTCTTg                             >  2:70107/1‑141 (MQ=255)
                                                                                                         aCCTGCCCACGTACGTGCAATGTTACCAGTTTTAGCGGTGATTACTGGTTGATATAAAAAGATATGTTTTGTCAGATTAGTATTCTCTCCTGGTGTGGATTGCTTTTGTCATTTTATTATGAACTTGTGCATCctccggag                           >  2:147467/1‑136 (MQ=255)
                                                                                                           cTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCCTCTTGATCa                         <  1:349436/141‑1 (MQ=255)
                                                                                                           cTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCCTCTTGATCa                         <  1:130506/141‑1 (MQ=255)
                                                                                                            tGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCCTCTTGATCAt                        >  1:96527/1‑141 (MQ=255)
                                                                                                                 acgtacgtGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCCTCTTGATCATAATTa                   <  2:239995/141‑1 (MQ=255)
                                                                                                                  cgtacgtGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCCTCTTGATCATAATTAg                  <  1:29160/141‑1 (MQ=255)
                                                                                                                                  ccAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCCTCTTGATCATAATTAGCATTGATAAAATctct  <  1:97967/141‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                           
ATACGTAATATTCACAAATTCCCAATAATCTAAACCGGCTCTTTTTAACATTTTATCATCAGTTCTAAATCCAAGAGGTTTAATTAAATGTAAATGTGTGTTTGTACCTGCACACGTACGTGCAATGTTACCAGTATTAGCTGGGATTTCTGGTTGATATAAAACGATATGATTTGTCATATTACTATTCTCTCCTTGTGTCTAATCCTTTTATCATTTCATTCTGAACTTCTGCATCCTCTTGATCATAATTAGCATTGATAAAATCTCT  >  CP000730/1988940‑1989210

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