Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,624,187	G→C	20.5%	intergenic (‑92/+81)	PP_3194 ← / ← PP_3195	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,624,187	0	G	C	20.5%	62.4 / 13.3	30	intergenic (‑92/+81)	PP_3194/PP_3195	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/10);  new base C (4/2);  total (18/12)
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.51e-01

CAGGCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTG  >  NC_002947/3624045‑3624334                                                                                                                                               |                                                                                                                                                    cAGACCAGTGCCGTGGCACAAGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGCTGGTGCAGGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTTATGGCCAGGTGTCGATAGCCTCGAt                                                                                                                                              >  2:670211/1‑150 (MQ=255)                                    ggtgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGgcg                                                                                                           <  2:733959/150‑1 (MQ=255)                                              aCTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGgcg                                                                                                 >  2:19269/1‑150 (MQ=255)                                              aCTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTg                                                                                                                                         >  1:678597/1‑110 (MQ=255)                                              aCTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTg                                                                                                                                         <  2:678597/110‑1 (MQ=255)                                                 cATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGcgccg                                                                                              >  2:266855/1‑150 (MQ=255)                                                      gTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCaa                                                                                         <  1:306077/150‑1 (MQ=255)                                                      gTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCaa                                                                                         <  2:411113/150‑1 (MQ=255)                                                           ctctctCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAg                                                                                    >  2:341685/1‑150 (MQ=255)                                                              tctcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTgg                                                                                                                                        <  2:783828/95‑1 (MQ=255)                                                              tctcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTgg                                                                                                                                        >  1:783828/1‑95 (MQ=255)                                                                             ggtgCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCa                                                                  <  1:19269/150‑1 (MQ=255)                                                                              gtgCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAg                                                                 <  1:266855/150‑1 (MQ=255)                                                                                             gtgGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGTCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGcc                                                                       >  1:20147/1‑129 (MQ=255)                                                                                             gtgGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGcc                                                                       <  2:20147/129‑1 (MQ=255)                                                                                                                 gcagATGCTAATGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTc                              >  2:314749/1‑150 (MQ=255)                                                                                                                     aTGCTATTGATGGCCAGGTATCGATGGACTCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAgatcgga                          >  1:20966/1‑144 (MQ=255)                                                                                                                        cTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATg                       >  1:669211/1‑150 (MQ=255)                                                                                                                           ttGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAg                    >  2:597473/1‑150 (MQ=255)                                                                                                                              aTGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTa                     <  1:10531/146‑1 (MQ=255)                                                                                                                              aTGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTa                     >  2:10531/1‑146 (MQ=255)                                                                                                                              aTGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGc                 <  1:597473/150‑1 (MQ=255)                                                                                                                                 gCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGACGTAGTGCCCg              <  1:706920/150‑1 (MQ=255)                                                                                                                                    aGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGTCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCa           >  1:361965/1‑150 (MQ=255)                                                                                                                                      gTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCAtt         >  2:623290/1‑150 (MQ=255)                                                                                                                                      gTATCGATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCAtt         >  1:528629/1‑150 (MQ=255)                                                                                                                                       tATCGATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTg        >  1:169877/1‑150 (MQ=255)                                                                                                                                           gATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGTCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGc    <  1:341685/150‑1 (MQ=255)                                                                                                                                           gATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGc    >  1:691516/1‑150 (MQ=255)                                                                                                                                             tGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGAGTTTCGTCGATCCAGTCGGCGGTGTAGTGCCCGGCATTGATGCTg  >  1:775524/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                    CAGGCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGCTG  >  NC_002947/3624045‑3624334

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