Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,670,626	T→C	10.4%	intergenic (+49/‑46)	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,626	0	T	C	10.4%	126.6 / 4.7	48	intergenic (+49/‑46)	rplS/PP_1466	50S ribosomal protein L19/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base T (21/22);  new base C (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 = 1.00e+00

CCTGGCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACCCT  >  NC_002947/1670478‑1670772                                                                                                                                                      |                                                                                                                                                   ccTGGTCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCt                                                                                                                                                    >  3:55406/1‑149 (MQ=255)     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)                                                                                                                                                      |                                                                                                                                                   CCTGGCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACCCT  >  NC_002947/1670478‑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: ‑