Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,312 G→A 64.0% intergenic (‑155/+288) rrsG ← / ← clpB 16S ribosomal RNA of rrnG operon/protein disaggregation chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,731,3120GA64.0% 51.0 / 59.5 64intergenic (‑155/+288)rrsG/clpB16S ribosomal RNA of rrnG operon/protein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base G (21/2);  new base A (18/23);  total (39/25)
Fisher's exact test for biased strand distribution p-value = 1.62e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAC  >  NC_000913/2731229‑2731406
                                                                                   |                                                                                              
cctgagatctccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAgcggcg                                                                                >  1:986902‑M1/76‑100 (MQ=255)
cctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAgcggcg                                                                                >  2:1322126‑M1/76‑100 (MQ=255)
cctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAgcggcg                                                                                >  1:1774628‑M1/76‑100 (MQ=255)
cctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAgcggcg                                                                                >  1:2382563‑M1/76‑100 (MQ=255)
cctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAgcggcg                                                                                >  1:2403811‑M1/76‑100 (MQ=255)
  tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGcgcg                                                                              >  2:2452401‑M1/74‑100 (MQ=255)
  tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGcgcg                                                                              >  2:110763‑M1/74‑100 (MQ=255)
  tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGcgcg                                                                              >  1:2143905‑M1/74‑100 (MQ=255)
  tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGcgcg                                                                              >  2:2020790‑M1/74‑100 (MQ=255)
  agagacctccggcagaaaaagcaaaaataaatgcttgactctgtagcgagaaggcgtattatgcacaccccgcGCCACTTCGCTTTTTCTCAGCGGcgcg                                                                              >  1:1789093‑M1/74‑100 (MQ=255)
   gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgg                                                                             >  2:673386‑M1/73‑100 (MQ=255)
   gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgg                                                                             >  2:730851‑M1/73‑100 (MQ=255)
   gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtatgatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgg                                                                             >  1:1736646‑M1/73‑100 (MQ=255)
   gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCgg                                                                             >  1:399210‑M1/73‑100 (MQ=255)
   gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCgg                                                                             >  2:22915‑M1/73‑100 (MQ=255)
    agaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgg                                                                             >  1:1468472‑M1/72‑99 (MQ=255)
    agaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCggg                                                                            >  1:50697‑M1/72‑100 (MQ=255)
     gaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgggg                                                                           >  1:954581‑M1/71‑100 (MQ=255)
     gaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCgggg                                                                           >  2:309696‑M1/71‑100 (MQ=255)
     gaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCgggg                                                                           >  1:941735‑M1/71‑100 (MQ=255)
       ttGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTTTTCGGCTCGGCgtg                                                                         >  1:449530/1‑100 (MQ=1)
                                cgcgTGCCCACACCGCTTGTCTGATCAACTGCTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTAc                                                  <  2:483026/95‑1 (MQ=1)
                                               aTTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTcat                                 <  2:1759290/100‑3 (MQ=2)
                                                ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt                                <  2:1373540/100‑1 (MQ=17)
                                                 tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                               <  1:1012704/100‑1 (MQ=18)
                                                 ggtcatgcgtattatgcacaccccgtGCCGCTTCACGGTTTCTCAGCGGCGTGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                               >  1:2388905‑M1/27‑100 (MQ=255)
                                                  gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                              <  1:384379/100‑1 (MQ=25)
                                                  gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                              >  2:2133235/1‑100 (MQ=25)
                                                   tCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCt                             <  2:2039152/100‑1 (MQ=255)
                                                   tCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGATACAGAGTCAAGCATTTCTTTGCt                             <  1:385642/100‑1 (MQ=255)
                                                    cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                            >  1:2545089/1‑100 (MQ=255)
                                                    cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                            >  1:708256/1‑100 (MQ=255)
                                                    cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                            >  1:2099100/1‑100 (MQ=255)
                                                    cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                            >  2:1350484/1‑100 (MQ=255)
                                                     tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGAGGCGAGGGGTGTGCATAATACGCCTTCTCGCTACAGAGTCAAGCCTTTCTTTGCttt                           >  1:1630198/1‑100 (MQ=14)
                                                        tAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtct                        >  2:2519061/1‑100 (MQ=255)
                                                          aaTTGTTAAAGAGCAGTGCCACTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTctct                      <  2:680372/100‑1 (MQ=255)
                                                           aTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTg                     >  1:637812/1‑100 (MQ=255)
                                                           aTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTg                     >  2:2375262/1‑100 (MQ=255)
                                                            ttGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGt                    <  2:756721/100‑1 (MQ=255)
                                                             tGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGtt                   <  2:715012/100‑1 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATTCGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  >  1:564117/1‑100 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  <  2:270618/100‑1 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  <  1:2021714/100‑1 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  >  2:2225052/1‑100 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  <  2:637815/100‑1 (MQ=255)
                                                              atgcacaccccgcGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  >  2:740368‑M1/14‑100 (MQ=255)
                                                               ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCACGGGGTGTGCATAATACGCTTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  >  1:101309/1‑99 (MQ=255)
                                                                tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag                >  2:495745/1‑100 (MQ=255)
                                                                  aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt              <  2:225987/100‑1 (MQ=255)
                                                                  aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt              <  1:1434720/100‑1 (MQ=255)
                                                                    gagCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtc            >  2:1339572/1‑100 (MQ=255)
                                                                     agCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtct           <  1:1194935/100‑1 (MQ=255)
                                                                        aGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg        <  1:89197/100‑1 (MQ=255)
                                                                         gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg       <  2:2145913/100‑1 (MQ=255)
                                                                         gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg       >  2:2112032/1‑100 (MQ=255)
                                                                         gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg       <  1:1291183/100‑1 (MQ=255)
                                                                         gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg       >  1:1056983/1‑100 (MQ=255)
                                                                          tGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGga      <  2:574366/100‑1 (MQ=255)
                                                                           GCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgag     <  1:2188894‑M1/100‑1 (MQ=255)
                                                                           GCCGCTTCACTTTTTCTCAGCGGCGCGGGGAGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgag     <  1:393772‑M1/100‑1 (MQ=255)
                                                                            CCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgaga    <  1:2498747‑M1/100‑1 (MQ=255)
                                                                              GCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAc  <  2:923658‑M1/100‑1 (MQ=255)
                                                                              GCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAc  <  2:1079234‑M1/100‑1 (MQ=255)
                                                                                   |                                                                                              
TTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAC  >  NC_000913/2731229‑2731406

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