Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,313,276	C→G	9.6%	intergenic (‑595/+119)	PP_16SG ← / ← msrQ	16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,313,276	0	C	G	9.6%	152.1 / 6.8	52	intergenic (‑595/+119)	PP_16SG/msrQ	16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit
Reads supporting (aligned to +/- strand):  ref base C (21/25);  new base G (2/3);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

GAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCGCG  >  NC_002947/5313137‑5313423                                                                                                                                             |                                                                                                                                                    gAAGTGGGGCGCATTATAACGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCgcggcg                                                                                                                                            >  8:189045/1‑149 (MQ=255)   aGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CGCGGGCGGCGCt                                                                                                                                          <  6:76496/149‑1 (MQ=255)   aGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCt                                                                                                                                          <  7:189045/149‑1 (MQ=255)     tGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCg                                                                                                                                        <  4:243316/149‑1 (MQ=255)       gggCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTCCAGCGCCTATGAGATCGAGCGCCG‑CGCGGGCGGCGCTCGAt                                                                                                                                      <  2:232209/149‑1 (MQ=255)        ggCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGAtc                                                                                                                                     <  2:72875/149‑1 (MQ=255)          cgcATTATAAGGGGATGCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTAGCGGGGCTGCCAAGTGCTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGAtctc                                                                                                                                   <  3:175177/149‑1 (MQ=255)           gcATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTca                                                                                                                                  >  7:192179/1‑149 (MQ=255)                    ggggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCa                                                                                                                         <  2:56178/149‑1 (MQ=255)                     gggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGc                                                                                                                           >  1:257010/1‑146 (MQ=255)                     gggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGc                                                                                                                           >  7:173639/1‑146 (MQ=255)                     gggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGc                                                                                                                           <  6:153963/146‑1 (MQ=255)                     gggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGc                                                                                                                           >  5:153963/1‑146 (MQ=255)                        aTTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAgaga                                                                                                                     <  5:61475/149‑1 (MQ=255)                                   cGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg                                                                                                          <  8:92631/149‑1 (MQ=255)                                       aaCCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGccc                                                                                                      >  2:62441/1‑149 (MQ=255)                                        aCCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAgaga                                                                                                                     <  1:30327/133‑1 (MQ=255)                                        aCCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAgaga                                                                                                                     >  2:30327/1‑133 (MQ=255)                                            ttAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCACCTATGAGATCGAGCGCCG‑CACGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCtctc                                                                                                  >  7:19521/1‑148 (MQ=255)                                                tttCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAg                                                                                             >  4:162871/1‑149 (MQ=255)                                                    aaGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCAt                                                                                         <  4:213091/149‑1 (MQ=255)                                                     aGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg                                                                                                          >  6:270139/1‑131 (MQ=255)                                                     aGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg                                                                                                          <  5:270139/131‑1 (MQ=255)                                                      gAAAGTTAAATATAAGGGTAAATACCTTATCGGGGCTTCCAAGTTTTCGCCACAAGAGTTTCAGCACCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTg                                                                                       <  8:19521/149‑1 (MQ=25)                                                       aaaGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGc                                                                                      <  3:162871/149‑1 (MQ=255)                                                            gAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAtc                                                                                  <  2:257010/148‑1 (MQ=255)                                                               atatAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCaa                                                                              <  5:50021/149‑1 (MQ=255)                                                                    aGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg                                                                                                          <  5:2508/116‑1 (MQ=255)                                                                    aGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg                                                                                                          >  6:2508/1‑116 (MQ=255)                                                                                gTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGc                                                             >  5:179742/1‑149 (MQ=255)                                                                                gTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGc                                                             >  4:141116/1‑149 (MQ=255)                                                                                 tATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAtct                                                                                 <  6:92093/128‑1 (MQ=255)                                                                                 tATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTAtct                                                                                 >  5:92093/1‑128 (MQ=255)                                                                                     ggggCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGccccac                                                        <  8:192179/149‑1 (MQ=255)                                                                                      gggCTGCCAAGTTTTCGCCACAACAGTTTCAGCGCCTATGAGATCGACCGCCG‑CCCGCGCGGCGCTAGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCGAGTGCTATCTCACGCCAGGCCCATGGCAcccccccc                                                       >  6:62822/1‑149 (MQ=25)                                                                                           gCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATCAGATCGAGCGCCG‑CCCCCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCAct                                                  <  3:22032/149‑1 (MQ=255)                                                                                            ccAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGACCGGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCAct                                                  >  1:67019/1‑149 (MQ=255)                                                                                            ccAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCAct                                                  >  3:201872/1‑148 (MQ=255)                                                                                               aGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACtctc                                               >  1:97886/1‑148 (MQ=255)                                                                                               aGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATATCAAGACAGGCCCTTGGCAGCCCCACCCCACtctc                                               >  8:61228/1‑148 (MQ=255)                                                                                                      gCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGc                                        <  7:158945/148‑1 (MQ=255)                                                                                                             caGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGGTCGATCTCACAGAGGCCAGAGAAACAACGGCGGGGCCACTCCACAGCCATTGCTATCTCAAGACCCGCCCTTCGCAGCCCGACCCCACTATCAATCCGCCCTAcc                                  >  5:131195/1‑147 (MQ=37)                                                                                                                tttCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTa                                                                                    >  6:255472/1‑94 (MQ=255)                                                                                                                tttCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTa                                                                                    <  5:255472/94‑1 (MQ=255)                                                                                                                 ttCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCCGCCCCACCCCCCTCTCCAACCGCCCTACCCTGCTa                            >  7:61894/1‑149 (MQ=255)                                                                                                                       gcCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGc                      <  4:201872/149‑1 (MQ=255)                                                                                                                       gcCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGc                      <  6:94879/149‑1 (MQ=255)                                                                                                                                 tCGAGCGCCG‑CGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTggg             >  1:71643/1‑148 (MQ=255)                                                                                                                                 tCGAGCGCCG‑CGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTggg             <  2:71643/148‑1 (MQ=255)                                                                                                                                    aGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAgc         <  2:97886/149‑1 (MQ=255)                                                                                                                                     gcgcCG‑CGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCGCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAgcg        >  5:232992/1‑149 (MQ=255)                                                                                                                                      cgccc‑cCCGCGCGGCGCGCTATCACACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAgcgc       <  3:12820/149‑1 (MQ=255)                                                                                                                                         cg‑cgcggGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTggg             >  1:98120/1‑140 (MQ=255)                                                                                                                                         cg‑cgcggGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTggg             <  2:98120/140‑1 (MQ=255)                                                                                                                                          g‑cccgcGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCcgc   <  3:141116/149‑1 (MQ=255)                                                                                                                                            cccgcGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCcgcg  >  1:19730/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                    GAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCG‑CCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGCGCCCGCG  >  NC_002947/5313137‑5313423

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