Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,670,629	G→A	11.1%	intergenic (+52/‑43)	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,629	0	G	A	11.1%	114.7 / 5.2	45	intergenic (+52/‑43)	rplS/PP_1466	50S ribosomal protein L19/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base G (20/20);  new base A (3/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-02

CCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACCCT  >  NC_002947/1670483‑1670772                                                                                                                                                    |                                                                                                                                                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)                                                                                                                                                    |                                                                                                                                                CCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAA‑AAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACCCT  >  NC_002947/1670483‑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: ‑