Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,312 G→A 70.1% 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,3120GA70.1% 86.2 / 58.7 80intergenic (‑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 (28/27);  total (49/31)
Fisher's exact test for biased strand distribution p-value = 6.95e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCC  >  NC_000913/2731220‑2731408
                                                                                            |                                                                                                
gtcttCATTTTTGGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCGCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCttt                                                                                              <  2:1153938/94‑1 (MQ=0)
       ctcctgagaactccggcagagacagcaaaattacatgctccactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAccgg                                                                                    >  1:412052‑M1/78‑100 (MQ=255)
       ctcctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAgcgg                                                                                    >  2:1536445‑M1/78‑100 (MQ=255)
        ttgtCGTCTTGGGAGGTTAAGAAGTCGTAGTTTCGAGTTCGCACACATATTGTCTGATATCTTGTTAAAGAGCAGGGCCGCTGCTCCTTTTCTcggcgg                                                                                    <  1:2875174/96‑1 (MQ=0)
         cctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAgcggcg                                                                                  >  2:2093503‑M1/76‑100 (MQ=255)
          ctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGgcgc                                                                                 >  1:1155706‑M1/75‑100 (MQ=255)
          ctgagaactccggcagagaaagcaaaaataaaggcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGgcgc                                                                                 >  2:1206072‑M1/75‑100 (MQ=255)
           tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGcgcg                                                                                >  2:399683‑M1/74‑100 (MQ=255)
           tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGcgcg                                                                                >  2:2825682‑M1/74‑100 (MQ=255)
           tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGcgcg                                                                                >  2:2160448‑M1/74‑100 (MQ=255)
           tcaaaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggtaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGcgcg                                                                                >  1:1764925‑M1/74‑100 (MQ=255)
            gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGcgcg                                                                                >  1:2565837‑M1/73‑99 (MQ=255)
            gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgg                                                                               >  2:125985‑M1/73‑100 (MQ=255)
            gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgg                                                                               >  2:1300516‑M1/73‑100 (MQ=255)
             agaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCggg                                                                              >  1:465212‑M1/72‑100 (MQ=255)
              gaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAgcggcg                                                                                  >  1:2066952‑M1/71‑95 (MQ=255)
              gaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgggg                                                                             >  2:2796782‑M1/71‑100 (MQ=255)
              gaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgggg                                                                             >  2:297563‑M1/71‑100 (MQ=255)
              gaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCgggg                                                                             >  2:2407603‑M1/71‑100 (MQ=255)
                        agagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCataat                                                                   >  1:524322‑M1/61‑100 (MQ=255)
                           agaaaagcaaagaaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACg                                                                >  2:1779852‑M1/58‑100 (MQ=255)
                            aaagaaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCataata                                                                  >  1:1022289‑M1/57‑97 (MQ=255)
                                    aataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGc                                                       >  2:907429‑M1/49‑100 (MQ=255)
                                                         ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt                                  <  1:517971/100‑1 (MQ=17)
                                                          tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                                 <  2:2520720/100‑1 (MQ=18)
                                                          tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                                 <  1:874716/100‑1 (MQ=18)
                                                           gTCTGATACATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                                <  1:1734580/100‑1 (MQ=255)
                                                           gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTTTTTTTTTTTCAGGGCGTGCATAATACGCCTTCTCTCTACAGAGTCAAGCATTTCTTTg                                 >  1:2222712/1‑99 (MQ=2)
                                                           gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                                <  1:2806362/100‑1 (MQ=25)
                                                           gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                                >  1:1884721/1‑100 (MQ=25)
                                                            tCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCt                               <  1:121128/100‑1 (MQ=38)
                                                            tCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCt                               >  2:1319461/1‑100 (MQ=255)
                                                            tCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCt                               >  2:2706216/1‑100 (MQ=255)
                                                             cTGATAAATTGTTCAAGAGCAGTGCCGCTTCACTTTTTGTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                              >  1:1796893/1‑100 (MQ=255)
                                                             cTGATAAATTGTTAATTAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTAAAGCATTTCTTTGCtt                              >  1:630189/1‑100 (MQ=14)
                                                             cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTTTTTTTTTTTTTTGGTGTGAATAATACGCCTTCCCGTTTCAGAGTTAAGCATTTCTTTGCtt                              >  1:960788/1‑100 (MQ=2)
                                                             cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTTTTTTCTTTGCTGAGTGTGCATAATACGCCTTCTTGCTACAATGGCAAGGATTTCTTTGCt                               >  1:630392/1‑99 (MQ=2)
                                                             cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                              >  1:279510/1‑100 (MQ=255)
                                                             cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                              >  2:518575/1‑100 (MQ=255)
                                                             cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAACACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                              <  2:1071866/100‑1 (MQ=255)
                                                              tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt                             >  1:895302/1‑100 (MQ=255)
                                                              tGATAAATTGATAAAGAGAAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGAGTGCATAATACGCCTTCACGCTACAGAGTCACGCATCTCTTTCCttt                             >  1:513964/1‑100 (MQ=11)
                                                               gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTTTTTTTTTTTTGGGGTGTGCATAATACGCCTTCCCTCTACATAGTCAAGCATTTCTTTGCtttt                            >  1:635361/1‑100 (MQ=9)
                                                               gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtttt                            >  2:646033/1‑100 (MQ=255)
                                                               gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtttt                            >  2:22343/1‑100 (MQ=255)
                                                               gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt                             >  1:1189053/1‑99 (MQ=255)
                                                                 tAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCCTAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtttt                            >  2:997442/1‑98 (MQ=255)
                                                                  aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                         <  2:279511/100‑1 (MQ=255)
                                                                  aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                         <  1:1962001/100‑1 (MQ=255)
                                                                  aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                         <  1:2109051/100‑1 (MQ=255)
                                                                    aTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTg                       >  2:1521622/1‑100 (MQ=255)
                                                                    aTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTg                       <  1:1676902/100‑1 (MQ=255)
                                                                     ttGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGt                      <  2:739002/100‑1 (MQ=255)
                                                                       gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                    <  2:1689193/100‑1 (MQ=255)
                                                                       gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                    >  2:2624993/1‑100 (MQ=255)
                                                                       gTTAAAGAGCAGCGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                    <  1:1156685/100‑1 (MQ=255)
                                                                        ttAATGAGGAGTGCCGCTTCTCTTTTTCTCAGCGGCGCGGGGTGTGCATGGTACGGCTGCTCGTTACAGAGTCAAGCATTGCTTTGCTTTTCTCTGTTga                   <  1:120515/100‑1 (MQ=17)
                                                                        ttAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                    >  2:1137776/1‑99 (MQ=255)
                                                                        ttAAAGAGCAGTGCCGCTTCACTTTTTTTTTTTTTCTCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtct                          >  1:1724223/1‑93 (MQ=14)
                                                                         tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag                  >  2:2080747/1‑100 (MQ=255)
                                                                          aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga                 >  1:1078825/1‑100 (MQ=255)
                                                                          aaaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga                 <  2:2606611/100‑1 (MQ=255)
                                                                           aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCTGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt                <  1:960849/100‑1 (MQ=255)
                                                                           aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGTCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga                 <  1:505991/99‑1 (MQ=255)
                                                                           aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt                >  1:622806/1‑100 (MQ=255)
                                                                             gagCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtc              <  2:2017946/100‑1 (MQ=255)
                                                                              agCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtct             >  1:835718/1‑100 (MQ=255)
                                                                              agCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtct             <  2:366876/100‑1 (MQ=255)
                                                                               gCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtctc            <  2:783905/100‑1 (MQ=255)
                                                                               gCAGTGCCGCTTCACTTTTTCTCAGAGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtctc            >  1:1209369/1‑100 (MQ=255)
                                                                                cAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCa           <  2:682635/100‑1 (MQ=255)
                                                                                 aGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg          <  1:2379424/100‑1 (MQ=255)
                                                                                   tGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGga        <  2:1832231/100‑1 (MQ=255)
                                                                                    GCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg          >  1:140062‑M1/1‑97 (MQ=255)
                                                                                    GCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgag       <  2:2643145‑M1/100‑1 (MQ=255)
                                                                                     CCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgaga      <  2:2377308‑M1/100‑1 (MQ=255)
                                                                                       GCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAc    <  1:1451911‑M1/100‑1 (MQ=255)
                                                                                       GCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAc    <  2:1511648‑M1/100‑1 (MQ=255)
                                                                                        CTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAcc   >  2:2387910‑M1/1‑100 (MQ=255)
                                                                                         TTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGCGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAccc  <  2:255311‑M1/100‑1 (MQ=255)
                                                                                            |                                                                                                
GTATTCATTTTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCC  >  NC_000913/2731220‑2731408

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