Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,731,330 G→A 9.1% intergenic (‑173/+270) rrsG ← / ← clpB 16S ribosomal RNA/ClpB chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,731,3300GA 9.1% 177.5 / 8.7 55intergenic (‑173/+270)rrsG/clpB16S ribosomal RNA/ClpB chaperone
Reads supporting (aligned to +/- strand):  ref base G (21/29);  new base A (2/3);  total (23/32)
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.96e-01

CGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCCGGCGGCGTGT  >  NC_000913/2731253‑2731427
                                                                             |                                                                                                 
cggaTCTTCGAGTGCCCACACAGATTGTCGAATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGGGGGGCGGGGTGTGCATGATACGCCTTc                                                                              <  2:131243/96‑1 (MQ=1)
                       aTTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCt                                                       >  2:1345283/1‑99 (MQ=32)
                       aTTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCt                                                       <  2:1446296/99‑1 (MQ=32)
                       aTTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCt                                                       <  2:791118/99‑1 (MQ=32)
                       aTTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCt                                                       <  1:129704/99‑1 (MQ=32)
                        ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCtt                                                      <  2:1239314/99‑1 (MQ=33)
                        ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCtt                                                      <  1:769541/99‑1 (MQ=33)
                        ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCtt                                                      >  2:996675/1‑99 (MQ=33)
                         tGTCTGGTAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt                                                     >  1:121028/1‑99 (MQ=21)
                         tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt                                                     <  1:1886202/99‑1 (MQ=33)
                          gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                                                    <  1:1260288/99‑1 (MQ=35)
                          gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCAGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg                                                    >  1:2162330/1‑99 (MQ=21)
                           tCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc                                                   <  2:1272956/99‑1 (MQ=37)
                            cTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCt                                                  <  1:852510/99‑1 (MQ=38)
                            cTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCt                                                  <  1:559924/99‑1 (MQ=38)
                             tGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt                                                 <  2:895034/99‑1 (MQ=39)
                               aTAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtttt                                               <  2:391872/99‑1 (MQ=255)
                                tAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtc                                              >  2:1064251/1‑99 (MQ=255)
                                tAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtc                                              <  2:1772368/99‑1 (MQ=255)
                                  aaTTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                                            >  1:939562/1‑99 (MQ=255)
                                  aaTTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCAGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc                                            >  2:205260/1‑99 (MQ=255)
                                   aTTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTctct                                           >  1:1043318/1‑99 (MQ=255)
                                     tGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGt                                         <  2:639319/99‑1 (MQ=255)
                                       ttAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                                       <  1:76296/99‑1 (MQ=255)
                                       ttAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCAGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg                                       <  1:1187154/99‑1 (MQ=255)
                                          aaGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga                                    >  2:697504/1‑99 (MQ=255)
                                          aaGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgaga                                    >  1:231921/1‑99 (MQ=255)
                                           agagCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt                                   >  1:1682523/1‑99 (MQ=255)
                                            gagCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAtt                                  >  2:818689/1‑99 (MQ=255)
                                            gagCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCAGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAtt                                  <  1:1084536/99‑1 (MQ=255)
                                                 gTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg                             <  2:1223927/99‑1 (MQ=255)
                                                 gTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg                             >  1:689549/1‑99 (MQ=255)
                                                 gTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg                             >  1:1128137/1‑99 (MQ=255)
                                                 gTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg                             >  2:1698796/1‑99 (MQ=255)
                                                  tGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg                            >  1:1726826/1‑99 (MQ=255)
                                                   gCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGga                           <  2:1524370/99‑1 (MQ=255)
                                                    ccGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgag                          <  2:1072338/99‑1 (MQ=255)
                                                      gcttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGaa                        <  1:1231771/99‑1 (MQ=255)
                                                       cttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAc                       <  2:911584/99‑1 (MQ=255)
                                                        ttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAcc                      <  1:1129861/99‑1 (MQ=255)
                                                        ttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAcc                      <  2:667048/99‑1 (MQ=255)
                                                        ttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAcc                      <  1:542580/99‑1 (MQ=255)
                                                        ttcgcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAcc                      <  1:848575/99‑1 (MQ=255)
                                                           gcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCccg                   >  1:342824/1‑99 (MQ=255)
                                                           gcttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCccg                   >  2:1339089/1‑99 (MQ=255)
                                                             tttttCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCccgcc                 >  1:719633/1‑99 (MQ=255)
                                                                ttCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGAc              >  1:499769/1‑99 (MQ=255)
                                                                ttCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGAc              <  2:970273/99‑1 (MQ=255)
                                                                ttCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGAc              >  2:2171607/1‑99 (MQ=255)
                                                                ttCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGAc              <  1:24277/99‑1 (MQ=255)
                                                                      gcggcgCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCcggcg        >  1:27253/1‑99 (MQ=255)
                                                                      gcggcgCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCcggcg        >  1:91690/1‑99 (MQ=255)
                                                                       cggcgCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCcggcgg       <  2:303930/99‑1 (MQ=255)
                                                                        ggcgCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCCggcggc      <  2:2106077/99‑1 (MQ=255)
                                                                         gcgcAGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCCGgcggcg     <  2:1226422/99‑1 (MQ=255)
                                                                          cgcgGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCCGGCGGCgt    <  2:2187769/99‑1 (MQ=255)
                                                                           gcgGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCCGGCGGCgtg   >  2:1516108/1‑99 (MQ=255)
                                                                            cgGGGTGTGCATAATACGCGTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAGCCCCGCCGACCCGGCGGCgtgt  <  2:2164820/99‑1 (MQ=255)
                                                                             |                                                                                                 
CGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAACCCCGCCGACCCGGCGGCGTGT  >  NC_000913/2731253‑2731427

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