Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,312 G→A 100% 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,3120GA84.6% 112.0 / 16.6 52intergenic (‑155/+288)rrsG/clpB16S ribosomal RNA of rrnG operon/protein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base G (8/0);  new base A (19/25);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 4.39e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

TTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAC  >  NC_000913/2731229‑2731406
                                                                                   |                                                                                              
cctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgc‑‑‑‑‑TTCGCTTTTTCTCAGCGGCGCgggg                                                                           >  1:231610‑M1/76‑100 (MQ=255)
cctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAgcggcg                                                                                >  2:2109012‑M1/76‑100 (MQ=255)
  tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGcgcg                                                                              >  2:951157‑M1/74‑100 (MQ=255)
    agaactccggcagagaaagcaaaaataaatgcttgtctctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCggg                                                                            >  1:1735855‑M1/72‑100 (MQ=255)
     gaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgggg                                                                           >  1:1119106‑M1/71‑100 (MQ=255)
        ctccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCGGGgtgt                                                                        >  1:1522183‑M1/68‑100 (MQ=255)
         tccggcagagaaagcaaaaataaatgcttgactctgtagtgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCGGGGtgtg                                                                       >  1:1096705‑M1/67‑100 (MQ=255)
         tccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGtgtg                                                                       >  2:5430‑M1/67‑100 (MQ=255)
                 aaaaaTCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATAc                                                               <  1:1846230/97‑1 (MQ=0)
                      aagcaaagaaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCtt                                                          >  2:686098‑M1/54‑100 (MQ=255)
                                              gATTGTCTGATAAATTGTTAAAGAGCAGTGTCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCt                                  <  1:2067870/100‑1 (MQ=9)
                                                ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt                                <  1:2506259/100‑1 (MQ=17)
                                                ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt                                <  2:577465/100‑1 (MQ=17)
                                                  gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                              <  1:1216727/100‑1 (MQ=25)
                                                    cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                            >  2:1055460/1‑100 (MQ=38)
                                                    cTGATAAATTGTTAAAGAGCAGTGCCGCTACACTTTTTCTCAGCGGCGCGGGGAGTGCATAATACGCCTTCCCGCTACAGAGACAAGCATTTCTTTGCtt                            >  2:2445888/1‑100 (MQ=14)
                                                     tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt                           <  1:1152928/100‑1 (MQ=255)
                                                     tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt                           <  2:2427847/100‑1 (MQ=255)
                                                     tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                            <  2:2226917/99‑1 (MQ=38)
                                                      gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtttt                          >  2:2520102/1‑100 (MQ=255)
                                                      gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt                           >  2:230747/1‑99 (MQ=255)
                                                       aTAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtc                         <  1:1851773/100‑1 (MQ=255)
                                                         aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCTGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                       <  2:2521902/100‑1 (MQ=255)
                                                         aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                       <  2:1899249/100‑1 (MQ=255)
                                                         aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                       >  1:2067601/1‑100 (MQ=255)
                                                         aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGAATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                       >  1:2409165/1‑100 (MQ=255)
                                                           aTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTg                     <  2:1921165/100‑1 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  <  2:122625/100‑1 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  <  2:496387/100‑1 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  <  2:2627545/100‑1 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  >  1:1688598/1‑100 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  <  1:571844/100‑1 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                  >  1:551653/1‑100 (MQ=255)
                                                              gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGACAAGCATTTCTTTGCTTTTCTCTGTTg                  >  1:1850822/1‑100 (MQ=255)
                                                               ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTga                 <  1:2333922/100‑1 (MQ=255)
                                                                tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag                <  1:1753713/100‑1 (MQ=255)
                                                                tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTga                 >  1:476650/1‑99 (MQ=255)
                                                                 aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga               >  1:1366899/1‑100 (MQ=255)
                                                                 aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga               >  1:1366966/1‑100 (MQ=255)
                                                                 aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCACGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga               >  2:1626886/1‑100 (MQ=255)
                                                                      gCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtctc          <  2:158965/100‑1 (MQ=255)
                                                                       cttccgatCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCa         >  2:265766‑M1/9‑100 (MQ=255)
                                                                       cAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCa         <  2:2305361/100‑1 (MQ=255)
                                                                         gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg       <  2:2075657/100‑1 (MQ=255)
                                                                         gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg       <  2:2069518/100‑1 (MQ=255)
                                                                         gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg        >  1:181104/1‑99 (MQ=255)
                                                                          tGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGga      <  2:1544230/100‑1 (MQ=255)
                                                                           GCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgag     >  2:181102‑M1/1‑100 (MQ=255)
                                                                           GCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgag     >  2:1347309‑M1/1‑100 (MQ=255)
                                                                            CCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACACAGTCAAACATTTCTTTGCTTTTCTATGTTGAGATTCTCAGgaga    <  1:245224‑M1/100‑1 (MQ=255)
                                                                              GCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAc  >  1:1358713‑M1/1‑100 (MQ=255)
                                                                              GCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAc  <  1:1056562‑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: