Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,989,572	C→A	15.4%	G79C (GGT→TGT)	PP_4398 ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,989,572	0	C	A	15.4%	63.2 / 2.7	26	G79C (GGT→TGT)	PP_4398	MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (14/8);  new base A (2/2);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TCGATGACCATGTATTGGCACATCATCCCCAGGCCGACGATCACGCGCAGCCCGACCCAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGCACCCACCCACAGGCTGTCGACATGGTCGGCAGCCAGGCGCAGGG  >  NC_002947/4989460‑4989718                                                                                                                   |                                                                                                                                                   tcGATGACCATGTATTGGCACATCATCCCCAGGCCGACGATCACGCGCAGCCCGACCCAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGt                                                                                                               >  2:377657/1‑150 (MQ=255)         cATGTATTGGCACATCATCCCCAGGCCGACGATCACGCGCAGCCCGACCCAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGc                                                                                                                        >  1:292720/1‑133 (MQ=255)         cATGTATTGGCACATCATCCCCAGGCCGACGATCACGCGCAGCCCGACCCAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGc                                                                                                                        <  2:292720/133‑1 (MQ=255)          acgtcttGGCACATCATCCCCAGGCCGACGATCACGAGCAGCAACA‑TCTCCTCGGCAAGCA‑CTGGCATCAGGCCATGCCCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGc                                                                                                      <  1:59638/145‑1 (MQ=255)                  gCACATCATCCCCAGGCCGGCGATCACGCGCAGCCCGACCCAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGcc                                                                                              >  1:501478/1‑150 (MQ=255)                   cacaTCATCCCCAGGCCGACGATCACGCGCAGCCCGACCCAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCa                                                                                             >  1:225798/1‑150 (MQ=255)                      atcatcCCCAGGCCGACGATCACGCGCAGCCCGACCCAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTACGTGCCCTACCCGTCCGATCAGCCGGTGGCCAAcc                                                                                          >  2:271256/1‑150 (MQ=255)                                       gaTCACGCGCAGCCCGACCCAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAg                                                                         >  2:411394/1‑150 (MQ=255)                                           aCGCGCAGCCCGACCCAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCa                                                                     >  2:36826/1‑150 (MQ=255)                                                         ccAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGtt                                                             >  1:397169/1‑144 (MQ=255)                                                         ccAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGtt                                                             <  2:397169/144‑1 (MQ=255)                                                           aGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACAGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCa                                                     >  1:452841/1‑150 (MQ=255)                                                                ggCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACAGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTgg                                                           >  2:431929/1‑139 (MQ=255)                                                                ggCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACAGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTgg                                                           <  1:431929/139‑1 (MQ=255)                                                                       agcTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGcaccca                                         >  1:227122/1‑150 (MQ=255)                                                                                     cATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTAc                                                                                                                   >  2:50267/1‑64 (MQ=255)                                                                                     cATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTAc                                                                                                                   <  1:50267/64‑1 (MQ=255)                                                                                     cATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGCACCCACCCACAGGCTGt                             >  2:357110/1‑150 (MQ=255)                                                                                           cGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGCACCCACCCACAGGCTGTCGACAt                       >  1:481536/1‑150 (MQ=255)                                                                                               cACCCCAGCGCCGACGACTCCGGCAACGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGCACCCACCCACAGGCTGTCGACATGGTc                   <  1:7578/150‑1 (MQ=255)                                                                                                aCCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGCACCCACCCACAGGCTGTCGACATGGTCg                  <  2:227122/150‑1 (MQ=255)                                                                                                          cgacgaTACAGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGCACCCACCCACAGGCTGTCGACATGGTCGGCAGCCAGgc        <  2:452841/150‑1 (MQ=255)                                                                                                           gacgaTACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGCACCCACCCACAGGCTGTCGACATGGTCGGCAGCCAGgcg       >  2:423106/1‑150 (MQ=255)                                                                                                             cgaTACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACATTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGCACCCACCCACAGGCTGTCGACATGGTCGGCAGCCAGGCGCa     <  1:377657/150‑1 (MQ=255)                                                                                                                tACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGCACCCACCCACAGGCTGTCGACATGGTCGGCAGCCAGGCGCAggg  >  2:199654/1‑150 (MQ=255)                                                                                                                tACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGCACCCACCCACAGGCTGTCGACATGGTCGGCAGCCAGGCGCAggg  <  1:423106/150‑1 (MQ=255)                                                                                                                   |                                                                                                                                                   TCGATGACCATGTATTGGCACATCATCCCCAGGCCGACGATCACGCGCAGCCCGACCCAGGCCGGCAGCCAGCTGG‑‑TCAGGCCATGACCGAGCACCGCAGCGCCGACGATACCGGCACAGGTGGCATAGGCGCGTATGTGCCCTACCCGTCCGATCAGCCGGTGGCCAACCTTGCCGCCAACCGCCAGACCAAAGTAGTTGGCCGCCATCAACGCACCCACCCACAGGCTGTCGACATGGTCGGCAGCCAGGCGCAGGG  >  NC_002947/4989460‑4989718

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