Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,078,789	G→T	14.8%	L165M (CTG→ATG)	PP_5330 ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,078,789	0	G	T	14.8%	68.0 / 4.2	27	L165M (CTG→ATG)	PP_5330	MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (11/12);  new base T (2/2);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

AGGTACACCGAGAACAGGCCGAACACACCGCCAAGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTCGCCGCCGATTGCCGCGCCTTGGATCACGCGCATCAGCAACAGCAGGATCGGTGC  >  NC_002947/6078641‑6078915                                                                                                                                                     |                                                                                                                               aGGTACACCGAGAACAGGCCGAAAAGTGGCCCAAGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGccggcggtcaggccggcggtcagg                                                                                                                               <  1:247374/150‑1 (MQ=255)      cacCGAGAACAGGCCGAACACACCGCCAAGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGcc                                                                                                                          <  2:4076/150‑1 (MQ=255)      cacCGACAACAGGCCGAACACACCGACCCGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGcc                                                                                                                          <  1:348332/150‑1 (MQ=255)       acCGAGAACAGGCCGAACACACCGCCAAGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCg                                                                                                                         >  1:2679/1‑150 (MQ=255)         cGAGAACAGGCCGAACACACCGCCAAGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGCACCCACTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCa                                                                                                                       >  1:59461/1‑150 (MQ=255)                  gCCGAACACACCGCCAAGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAg                                                                                                                      >  2:452792/1‑142 (MQ=255)                  gCCGAACACACCGCCAAGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAg                                                                                                                      <  1:452792/142‑1 (MQ=255)                  gCCGAACACACCGCCAAGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGcc                                                                                                              >  1:323292/1‑150 (MQ=255)                   ccGAACACACCGCCAAGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCg                                                                                                             >  2:268322/1‑150 (MQ=255)                           acCGCCAAGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCg                                                                                                     >  2:61502/1‑150 (MQ=255)                                   ggaagAAGGGGATTAGCCAGGCGTAGGCGGCCACGTCGGCCGCGCTGTACCCGCAGGTGATAAGGGTCGCTATTAGCGAGCCCAGTACAATGCCGGCGGTGAGGCCGGCGGTCAGGATGCCGCAGGCATAGCCGGTGTTTCGAGCCGGCa                                                                                             <  2:11502/146‑1 (MQ=255)                                        aaCGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACgtgt                                                                                        >  2:196460/1‑150 (MQ=255)                                                   ccAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCATGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCa                                                                                     <  1:318155/142‑1 (MQ=255)                                                   ccAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCATGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCa                                                                                     >  2:318155/1‑142 (MQ=255)                                                       gCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGGGGTCAGGCCGGCGGTCAGGGTGACGCAGGCCTAGCCGGTGTTGCCAGTCGGCACGTGTTCAGAGACTAACAcc                                                                         <  2:199183/150‑1 (MQ=255)                                                            gTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATgc                                                                    >  2:142882/1‑150 (MQ=255)                                                                 ccACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCATGGTGCCGCAg                                                                                                                      <  2:400964/95‑1 (MQ=255)                                                                 ccACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCATGGTGCCGCAg                                                                                                                      >  1:400964/1‑95 (MQ=255)                                                                             cgcTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTcgccgc                                                   <  1:196460/150‑1 (MQ=255)                                                                                  tACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTCGCCGCCGAtt                                              >  2:186198/1‑150 (MQ=255)                                                                                  tACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTCGCCGCCGAtt                                              >  1:195528/1‑150 (MQ=255)                                                                                                   ggTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTCGCCGCCGATTGCCGCGCCTTGGATCAc                             <  2:2679/150‑1 (MQ=255)                                                                                                    gTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTCGCCGCCGATTGCcgc                                         >  1:43344/1‑137 (MQ=255)                                                                                                    gTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTCGCCGCCGATTGCcgc                                         <  2:43344/137‑1 (MQ=255)                                                                                                     tCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTCGCCGCCGATTGCCGCGCCTTGGATCAcgc                           <  2:195528/150‑1 (MQ=255)                                                                                                     tCGCTACCAGCGAGCCCAGCAGAATGCCGGAGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTCCAGAGCCGACCCACCATGGGTCCGTTACTTCGCCGCCGATTGCCGCCCGTTGGagttctg                           <  2:327240/150‑7 (MQ=255)                                                                                                                   cccAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTCGCCGCCGATTGCCGCGCCTTGGATCACGCGCATCAGCAAcagc             <  1:186198/150‑1 (MQ=255)                                                                                                                    ccAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTCGCCGCCGATTGCCGCGCCTTGGATCACGCGCATCAGCAAcagca            <  2:439395/150‑1 (MQ=255)                                                                                                                              gccggcggtcaggccggcggtcagGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTCGCCGCCGATTGCCGCGCCTTGGATCACGCGCATCAGCAACAGCAGGATCGGTGc  <  1:61502/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                               AGGTACACCGAGAACAGGCCGAACACACCGCCAAGCAGGAACGGGATTCGCCAGGCGTAGTCAGCCACTTCTTCCGCGCTGTACACGCTGTTGATCAAGGTCGCTACCAGCGAGCCCAGCAGAATGCCGGCGGTCAGGCCGGCGGTCAGGGTGCCGCAGGCATAGCCGGTGTTGCGAGCCGGCACGTGTTCAGAGACGAACACCCATGCGCCCGGTACTTCGCCGCCGATTGCCGCGCCTTGGATCACGCGCATCAGCAACAGCAGGATCGGTGC  >  NC_002947/6078641‑6078915

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: ‑