Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,418,015	G→A	9.7%	A352A (GCC→GCT)	PP_4758 ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,418,015	0	G	A	9.7%	109.4 / 2.5	41	A352A (GCC→GCT)	PP_4758	MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (17/20);  new base A (2/2);  total (19/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 = 1.81e-01

GTGGGGGCAACGAAACCGCCGAGGTTGCCGACCGAGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCCTGCTGGAACTTCCAGCGAGCGGAACCGGCGGCGAAGGCGTAC  >  NC_002947/5417872‑5418141                                                                                                                                                |                                                                                                                                gtggGGGCAACGAAACCGCCGAGGTTGCCGACCGAGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGcgac                                                                                                                            <  1:346378/149‑1 (MQ=255)     gggCAACGAAACCGCCGAGGTTGCCGACCGAGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAg                                                                                                                        <  8:194675/149‑1 (MQ=255)       gCAACGAAACCGCCGAGGTTGCCGACCGAGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAgg                                                                                                                       >  5:67741/1‑148 (MQ=255)         aaCGAAACCGCCGAGGTTGCCGACCGAGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCa                                                                                                                     >  3:232906/1‑148 (MQ=255)            gAAACCGCCGAGGTTGCCGACCGAGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAaca                                                                                                                 >  7:199246/1‑149 (MQ=255)                 cgccgAGGTTGCCGACCGAGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑Gc                                                                                                           <  2:300856/149‑1 (MQ=255)                 cgccgAGGTTGCCGACCGAGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑Gc                                                                                                           <  6:67741/149‑1 (MQ=255)                         ttGCCGACCGAGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCgt                                                                                                   <  1:42942/149‑1 (MQ=255)                           gCCGACCGAGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCgtg                                                                                                  <  4:75744/148‑1 (MQ=255)                             cgaccgaGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCgtggtg                                                                                               <  8:62117/149‑1 (MQ=255)                               accgaGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGa                                                                                             >  7:330004/1‑149 (MQ=255)                                   aGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGAc                                                                                            <  6:132842/146‑1 (MQ=255)                                   aGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGAc                                                                                            >  5:132842/1‑146 (MQ=255)                                     ttGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGaa                                                                                       <  3:332110/149‑1 (MQ=255)                                         ccagggctatctggacgcccgcatcgctgcggcggtgattgccctgGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATg                                                                                   >  2:266297‑M2/47‑149 (MQ=255)                                          cagggctatctggacgcccgcatcgctgcggcggtgattgccctgGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGc                                                                                  >  4:196880‑M2/46‑149 (MQ=255)                                           agggctatctggacgcccgcatcgctgcggcggtgattgccctgGGGGATTGGCCAGAACAGTGACGCCGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAACCACCG‑GCCCACGCGTGGGGGACATGAACATTcc                                                                                 >  8:63922‑M2/45‑149 (MQ=255)                                                     ccgccgCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAgg                                                                                                                       <  7:254299/102‑1 (MQ=255)                                                     ccgccgCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAgg                                                                                                                       >  8:254299/1‑102 (MQ=255)                                                             gCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACATTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACcagcag                                                               <  8:199246/149‑1 (MQ=255)                                                              catcgctgcggcggtgattgccctgGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAgg                                                              <  3:196880‑M2/124‑1 (MQ=255)                                                                      cggcggtgattgccctgGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCGTGCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCg                                                        <  1:266297‑M2/131‑1 (MQ=255)                                                                          ccAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCAcc                                                  >  3:19678/1‑149 (MQ=255)                                                                           cAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCAcc                                                  <  4:232906/148‑1 (MQ=255)                                                                            aGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACcagc                                                                 <  8:61804/132‑1 (MQ=255)                                                                            aGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACcagc                                                                 >  7:61804/1‑132 (MQ=255)                                                                                 tccctgGGGTATTGGCAATAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACAGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATTAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCct                                           <  6:224361‑M2/143‑1 (MQ=255)                                                                                   cctgGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACAGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCTGCCACCCAGGCct                                           >  5:224361‑M2/5‑147 (MQ=255)                                                                                   cctgGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACAGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCct                                           <  6:229533‑M2/143‑1 (MQ=255)                                                                                   cctgGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACAGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCct                                           >  5:229533‑M2/5‑147 (MQ=255)                                                                                          gATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCctgc                                         >  5:8747/1‑142 (MQ=255)                                                                                          gATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCctgc                                         <  6:8747/142‑1 (MQ=255)                                                                                            ttGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCCTGCTGGAACTTc                                <  2:112501/149‑1 (MQ=255)                                                                                                 cAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCCTGCTGGAACTTCCagc                            >  4:185372/1‑148 (MQ=255)                                                                                                     aCAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCCTGCTGGAACTTCCagcgagc                        <  4:19678/148‑1 (MQ=255)                                                                                                        gTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCCTGCTGGAACTTCCAGCGAGCGGaa                    >  5:358531/1‑149 (MQ=255)                                                                                                               cgCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCCTGCTGGAACTTCCAGCGAGCGGAACcggcgg             >  1:158901/1‑149 (MQ=255)                                                                                                                    aCTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCCTGCTGGAACTTCCAGCGAGCGGAACCGGCGGCGAAg        <  6:358531/149‑1 (MQ=255)                                                                                                                        gAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCCTGCTGGAACTTCCAGCGAGCGGAACCGGCGGCGAAGGCg     >  3:219074/1‑148 (MQ=255)                                                                                                                                              acGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCCTGCTGGAACTTCCAGCGAGCGGAACCGGCGGCGAAGGCGTAc  <  6:327964/129‑1 (MQ=255)                                                                                                                                              acGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCCTGCTGGAACTTCCAGCGAGCGGAACCGGCGGCGAAGGCGTAc  >  5:327964/1‑129 (MQ=255)                                                                                                                                                |                                                                                                                                GTGGGGGCAACGAAACCGCCGAGGTTGCCGACCGAGTTGATCAGGGCAATCACCGCCGCAGCGATGCGGGCGTCCAGATAGCCCTGGGGGATTGGCCAGAACAGTGACGACGCTGACTTGAAGCCGATCGCGGCGAAGCACACGGCGACGAAGGCAAACACCG‑GCCCACCCGTGGTGGACATGAACATGCCCGTGGCAGCGACCAGCAGGGCGCCGGCCACCCAGGCCTGCTGGAACTTCCAGCGAGCGGAACCGGCGGCGAAGGCGTAC  >  NC_002947/5417872‑5418141

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