Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,670,627	T→G	10.4%	intergenic (+50/‑45)	rplS → / → PP_1466	50S ribosomal protein L19/acyl‑CoA thioesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,670,627	0	T	G	10.4%	135.3 / 8.1	48	intergenic (+50/‑45)	rplS/PP_1466	50S ribosomal protein L19/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base T (21/22);  new base G (3/2);  total (24/24)
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.57e-01

GCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACCCT  >  NC_002947/1670482‑1670772                                                                                                                                                   |                                                                                                                                                  gCCGTGAAACGGCGTGGTGACGTGCGTAAAGCCAAGATGTACTACCTGAGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTc                                                                                                                                                  <  2:18565/147‑1 (MQ=255)  ccGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCggg                                                                                                                                               <  5:106936/149‑1 (MQ=255)  ccGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCggg                                                                                                                                               >  3:20411/1‑149 (MQ=255)                  tgACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGttt                                                                                                                               <  8:224135/149‑1 (MQ=255)                        gCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAAc                                                                                                                         >  7:68661/1‑149 (MQ=255)                             aaGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCg                                                                                                                     >  2:2114/1‑148 (MQ=255)                              aGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGa                                                                                                                    >  6:134407/1‑148 (MQ=255)                                   aGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAAtgt                                                                                                              <  3:101420/149‑1 (MQ=255)                                             ccTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCCGAAGGCTGCttttttct                                                                                                                                   >  5:257335/1‑116 (MQ=39)                                             ccTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCCGAAGGCTGCttttttct                                                                                                                                   <  6:257335/118‑3 (MQ=39)                                              cTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATga                                                                                                   <  4:20411/149‑1 (MQ=255)                                                 cgcgACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTTCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATgac                                                                                                  <  3:21589/147‑1 (MQ=255)                                                        tGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCgaga                                                                                         >  4:258276/1‑149 (MQ=255)                                                             ggCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGCCGACCCGAGACCAg                                                                                     >  4:67873/1‑148 (MQ=255)                                                             ggCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAg                                                                                     >  7:172552/1‑148 (MQ=255)                                                               cAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGa                                                                                   >  4:227360/1‑148 (MQ=255)                                                                aaaGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGaaa                                                                                 >  6:554/1‑149 (MQ=255)                                                                   gCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATcc                                                                              <  5:554/149‑1 (MQ=255)                                                                              cAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCgaga                                                                                         <  1:162244/127‑1 (MQ=255)                                                                              cAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCgaga                                                                                         >  2:162244/1‑127 (MQ=255)                                                                               aaGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGATAATGAGATTGCCATGACGAACCGAGACCCGGAAATCCAGCGCCGCAcc                                                                   >  8:266900/1‑148 (MQ=255)                                                                                    aaaaCTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTg                                                             >  2:41930/1‑149 (MQ=255)                                                                                           tCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGa                                                      >  4:288858/1‑149 (MQ=255)                                                                                             cTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGAc                                                     <  3:258276/148‑1 (MQ=255)                                                                                                aGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCcgc                                                 <  6:34905/149‑1 (MQ=255)                                                                                                    cAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGtt                                              <  4:55406/148‑1 (MQ=255)                                                                                                     aGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTg                                                       <  7:293775/138‑1 (MQ=255)                                                                                                     aGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTg                                                       >  8:293775/1‑138 (MQ=255)                                                                                                       tttGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAGAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGa                                                                                                            >  1:202154/1‑84 (MQ=255)                                                                                                       tttGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGa                                                                                                            <  2:202154/83‑1 (MQ=255)                                                                                                            ccGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCg                                     >  6:298363/1‑149 (MQ=255)                                                                                                             cGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCgg                                    <  8:68661/149‑1 (MQ=255)                                                                                                                   cctcggCCACCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTcc                              <  8:225041/145‑1 (MQ=255)                                                                                                                       ggCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCaa                          <  3:288858/149‑1 (MQ=255)                                                                                                                         ccgccTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATa                        <  5:50674/149‑1 (MQ=255)                                                                                                                            ccTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCga                                                                                           <  4:87616/79‑1 (MQ=255)                                                                                                                            ccTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCga                                                                                           >  3:87616/1‑79 (MQ=255)                                                                                                                            ccTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATacca                     <  5:134407/149‑1 (MQ=255)                                                                                                                             cTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccac                    <  8:51515/149‑1 (MQ=255)                                                                                                                                  gagaAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCaa                          <  6:196388/138‑1 (MQ=255)                                                                                                                                  gagaAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCaa                          >  5:196388/1‑138 (MQ=255)                                                                                                                                  gagaAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaacc               >  1:183279/1‑149 (MQ=255)                                                                                                                                    gaAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaacc               <  3:227360/147‑1 (MQ=255)                                                                                                                                     aaaa‑aaGCAGCCCGAAGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCg                                     >  4:301940/1‑124 (MQ=255)                                                                                                                                     aaaa‑aaGCAGCCCGAAGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCg                                     <  3:301940/124‑1 (MQ=255)                                                                                                                                      aaa‑aaGCAGCCCGAAGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaaccacca           >  6:29303/1‑149 (MQ=255)                                                                                                                                            gcagcCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCacaaca      >  3:249328/1‑149 (MQ=255)                                                                                                                                                ccTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACCCt  >  2:244369/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACCCT  >  NC_002947/1670482‑1670772

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