Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,268,241	G→A	8.9%	L354L (CTG→TTG)	PP_3740 ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,268,241	0	G	A	8.9%	121.8 / 2.4	45	L354L (CTG→TTG)	PP_3740	MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (18/23);  new base A (2/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

TCGAGCAGCACGCCAGCCACGGTCGGCGACAGAATCGCGCCCATGCGGCCGATGCCGATCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTGGCCAGCGCGAAGCCAGGCCGCCGATCAGCGCC  >  NC_002947/4268093‑4268386                                                                                                                                                      |                                                                                                                                                  tCGAGCAGCACGCCAGCCACGGTCGGCGACAGAATCGCGCCCATGCGGCCGATGCCGATCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAg                                                                                                                                                   >  2:53746/1‑149 (MQ=255) tCGAGCAGCACGCCAGCCACGGTCGGCGACAGAATCGCGCCCATGCGGCCGATGCCGATCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAg                                                                                                                                                   >  2:240334/1‑149 (MQ=255)            gccagccaCGGTCGGCGACAGAATCGCGCCCATGCGGCCGATGCCGATCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATcag                                                                                                                                          >  6:307013/1‑147 (MQ=255)                gccaCGGTCGGCGACAGAATCGCGCCCATGCGGCCGATGCCGATCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGcc                                                                                                                                    >  8:177766/1‑149 (MQ=255)                gccaCGGTCGGCGACAGAATCGCGCCCATGCGGCCGATGCCGATCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGcc                                                                                                                                    <  8:230168/149‑1 (MQ=255)                            gACAGAATCGCGCCCATGCGGCCGATGCCGATCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAACAGGCCCAACGCCAgcgc                                                                                                                        >  8:213685/1‑149 (MQ=255)                                  aTCGCGCCCATGCGGCCGATGCCGATCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCCGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGt                                                                                                                  >  5:287553/1‑149 (MQ=255)                                   tCGCGCCCATGCGGCCGATGCCGATCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTg                                                                                                                 >  8:186999/1‑149 (MQ=255)                                         ccATGCGGCCGATGCCGATCCCCCAACCCACACCGGTGGCGAGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTCGAAAACAGCCCGATCAGCAGGCCCAACGCCa                                                                                                                            >  7:177364/1‑132 (MQ=255)                                                  cgatgccgatCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTg                                                                                                  <  7:286918/149‑1 (MQ=255)                                                  cgatgccgatCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTg                                                                                                  <  3:282542/149‑1 (MQ=255)                                                   gatgccgatCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTgg                                                                                                 <  7:177766/149‑1 (MQ=255)                                                        cgatCCCCCAACCCACACCGCTGGCCCGTACCGATGCGTCATAGACCACCGGCGATAACGCATACAGCCCAGCCACGCAGC‑CGATGGAAAACAGCCCCAACAGCAGGCCCAACGACAGCCCCGCAGAGACCGAGGCACCATTGGCCAc                                                                                             >  2:158129/1‑148 (MQ=255)                                                         gatCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCa                                                                                                                            >  6:226113/1‑116 (MQ=255)                                                               ccAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGTCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCa                                                                                     <  7:324718/149‑1 (MQ=255)                                                               ccAACCCACACAGGTTGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCa                                                                                     <  7:319717/149‑1 (MQ=255)                                                                    ccACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACcag                                                                                 >  3:125345/1‑148 (MQ=255)                                                                     cacaCCGGTGGCGCGTAGCGATGCGTCATAGACCACCGGCGACAAGGCATTCAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACcagca                                                                               <  8:66323/149‑1 (MQ=255)                                                                        acCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGCGGGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCagcagc                                                                             >  1:302650/1‑149 (MQ=255)                                                                             tGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCgcggcg                                                                       >  7:86448/1‑149 (MQ=255)                                                                                 gcgTACCGATGCGTCATACACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGtga                                                                   <  7:338676/149‑1 (MQ=255)                                                                                  cgTACCGATGCGTCATAGACCACCGGCGATAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGtgat                                                                  <  1:158129/149‑1 (MQ=255)                                                                                  cgTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGtgat                                                                  <  1:240334/149‑1 (MQ=255)                                                                                        gATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCAt                                                            <  2:8093/149‑1 (MQ=255)                                                                                         aTGCGTCATAGACGACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGTGGCACCATTGGCCACGAACAGCACCagcagc                                                                             <  1:334008/131‑1 (MQ=255)                                                                                                aTAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGaacaac                                                     <  4:125345/148‑1 (MQ=255)                                                                                                aTAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGaacaac                                                     <  8:86448/148‑1 (MQ=255)                                                                                                 tAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGc                                                   <  1:53746/149‑1 (MQ=255)                                                                                                     ccaccGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAg                                                >  6:60829/1‑148 (MQ=255)                                                                                                      caccGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCa                                              >  6:20142/1‑149 (MQ=255)                                                                                                           gCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCAcgc                                           >  8:321218/1‑147 (MQ=255)                                                                                                           gCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGa                                         <  6:59133/149‑1 (MQ=255)                                                                                                                 aGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAACCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGt                                   >  4:112950/1‑149 (MQ=255)                                                                                                                 aGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAACCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGt                                   >  4:107631/1‑149 (MQ=255)                                                                                                                  ggCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTg                                  <  5:307013/149‑1 (MQ=255)                                                                                                                    cATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTGGc                                <  7:86914/149‑1 (MQ=255)                                                                                                                     aTACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTGGc                                >  6:186375/1‑148 (MQ=255)                                                                                                                     aTACAGCCCAGCCACGCAGC‑CGATGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGACGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTGGcc                               >  2:80208/1‑149 (MQ=255)                                                                                                                          gcccagccACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCAt                                                            <  1:232524/115‑1 (MQ=255)                                                                                                                          gcccagccACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCAt                                                            >  2:232524/1‑115 (MQ=255)                                                                                                                                   cGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACCAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACACCGCCAGCACGCGAGTCAGTGGCCAGCGCGAAGCCCgg                 >  8:331253/1‑149 (MQ=255)                                                                                                                                    gCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTGGCCAGCGCGAAGCCAGgc                <  7:321218/149‑1 (MQ=255)                                                                                                                                       gC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTGGCCAGCGCGAAGCCAGgccgc             <  6:287553/149‑1 (MQ=255)                                                                                                                                              ggAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTGGCCAGCGCGAAGCCAGGCCGCCGATCa       <  2:302650/149‑1 (MQ=255)                                                                                                                                                aaaaCAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGGCCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTGGCCAGCGCGAAGCCAGGCCGCCGATCAgc     <  7:186999/149‑1 (MQ=255)                                                                                                                                                   aCAACCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTGGCCAGCGCGAAGCCAGGCCGCCGATCAgcgcc  <  3:107631/149‑1 (MQ=255)                                                                                                                                                   aCAACCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTGGCCAGCGCGAAGCCAGGCCGCCGATCAgcgcc  <  3:112950/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  TCGAGCAGCACGCCAGCCACGGTCGGCGACAGAATCGCGCCCATGCGGCCGATGCCGATCCCCCAACCCACACCGGTGGCGCGTACCGATGCGTCATAGACCACCGGCGACAAGGCATACAGCCCAGCCACGCAGC‑CGTTGGAAAACAGCCCGATCAGCAGGCCCAACGCCAGCGCCGCAGTGACCGAGGCACCATTGGCCACGAACAGCACCAGCAGCGCGGCGGTGATGAGCATGAACAACGCCAGCACGCGAGTCAGTGGCCAGCGCGAAGCCAGGCCGCCGATCAGCGCC  >  NC_002947/4268093‑4268386

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