Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,123,948 T→A 16.2% intergenic (‑143/+78) PP_mr13 ← / ← gcvT‑I cspA regulatory region/aminomethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,123,9480TA16.2% 89.7 / 12.3 37intergenic (‑143/+78)PP_mr13/gcvT‑IcspA regulatory region/aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base T (17/14);  new base A (2/4);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

CCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCACGCAC  >  NC_002947/1123801‑1124097
                                                                                                                                                   |                                                                                                                                                     
ccAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCAAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGcctggcctgg                                                                                                                                                     >  1:157621/1‑150 (MQ=255)
      cACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGTCGATCACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCGCCAt                                                                                                                                               <  1:185519/150‑1 (MQ=255)
      cACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCAGGCCAGGCGCCAt                                                                                                                                               <  1:163359/150‑1 (MQ=255)
         gAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGtt                                                                                                                                            >  2:11779/1‑150 (MQ=255)
                   gACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTAt                                                                                                                                  >  2:735693/1‑150 (MQ=255)
                    aCCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTAtt                                                                                                                                 <  2:252432/150‑1 (MQ=255)
                      ccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTgg                                                                                                                               <  1:748299/150‑1 (MQ=255)
                      ccATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTgg                                                                                                                               >  1:600614/1‑150 (MQ=255)
                         tGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCa                                                                                                                            >  1:644064/1‑150 (MQ=255)
                              acacAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTc                                                                                                                       <  1:11779/150‑1 (MQ=255)
                                    tgcagtgACTCTACGCAAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAAc                                                                                                                 <  2:157621/150‑1 (MQ=255)
                                         tgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGt                                                                                                                                     <  2:313862/125‑1 (MQ=255)
                                         tgACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGt                                                                                                                                     >  1:313862/1‑125 (MQ=255)
                                             tctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATg                                                                                                          >  1:123449/1‑148 (MQ=255)
                                             tctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATg                                                                                                          <  2:123449/148‑1 (MQ=255)
                                             tctACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTc                                                                                                        >  1:59068/1‑150 (MQ=255)
                                                  gCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTAtgc                                                                                                   <  2:437911/150‑1 (MQ=255)
                                                            gAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGc                                                                                                                             >  2:110986/1‑114 (MQ=255)
                                                            gAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGc                                                                                                                             <  1:110986/114‑1 (MQ=255)
                                                                gAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTc                                                                                     >  1:493407/1‑150 (MQ=255)
                                                                gAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCAACTTGTCGAAAATGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTc                                                                                     >  2:108384/1‑150 (MQ=255)
                                                                                gTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATAGGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGc                                                                     >  2:421776/1‑150 (MQ=255)
                                                                                               ggCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtg                                                      <  2:289131/150‑1 (MQ=255)
                                                                                                      tGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCAc                                               <  1:735693/150‑1 (MQ=255)
                                                                                                         cGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGaa                                            >  1:121736/1‑150 (MQ=255)
                                                                                                          gATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGaaa                                           <  2:59068/150‑1 (MQ=255)
                                                                                                               ccgaccgaTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCAt                                      <  1:428244/150‑1 (MQ=255)
                                                                                                                      aTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtg                                                      <  1:437283/127‑1 (MQ=255)
                                                                                                                      aTTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACgtg                                                      >  2:437283/1‑127 (MQ=255)
                                                                                                                               aaaCGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACgg                                                                <  2:648174/108‑1 (MQ=255)
                                                                                                                               aaaCGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACgg                                                                >  1:648174/1‑108 (MQ=255)
                                                                                                                                   aacATGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGAc                             >  2:151936/3‑139 (MQ=255)
                                                                                                                                   aacATGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGAc                             <  1:151936/137‑1 (MQ=255)
                                                                                                                                    acATGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACgg                                                                >  2:674536/2‑103 (MQ=255)
                                                                                                                                    acATGGCGCCAGGCCAGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACgg                                                                <  1:674536/102‑1 (MQ=255)
                                                                                                                                      aTGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACct               >  2:697424/1‑150 (MQ=255)
                                                                                                                                      aTGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACct               >  2:502618/1‑150 (MQ=255)
                                                                                                                                      aTGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACct               >  1:718962/1‑150 (MQ=255)
                                                                                                                                             ctggcctggcGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTcc        <  2:600614/150‑1 (MQ=255)
                                                                                                                                              tggcctggcGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCa       <  1:697424/150‑1 (MQ=255)
                                                                                                                                              tggcctggcGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCa       <  1:502618/150‑1 (MQ=255)
                                                                                                                                                  ctggcGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACNAATATGNCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCACGca   >  1:436079/1‑150 (MQ=255)
                                                                                                                                                   tggcGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCACGcac  >  1:103972/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCAGGTCACGAACTGTTGCGACCCATGCAAACACAGTGCAGTGACTCTACGCCAACTCTCGAGCGAAAAACAAGCCCTTGGTCGTATGGAAAACTGGCGGTTTGCCGATAACCGACCGATTTGTCGGAAACGGCATGGCGCCTGGCCTGGCGCCATGTTCAAAAGTTATTGGGCAGGTTCGAAAACGAATATGTCTATGCTGCCCAAAACCTTCACAAAGCCGTTTCAGCCACGGTAGTAACGTGGTGTCACGAAAGGCATTTTGCTGACTTTCAAGGCCACCTTCTTTCCACGCAC  >  NC_002947/1123801‑1124097

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