Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,268,885	T→A	10.0%	E139V (GAA→GTA) ‡	PP_3740 ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,268,885	0	T	A	10.0%	133.1 / 5.8	50	E139V (GAA→GTA) ‡	PP_3740	MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (21/24);  new base A (3/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

ACCAGTACGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAGCGCCAGA  >  NC_002947/4268741‑4269028                                                                                                                                                 |                                                                                                                                                aCCAGTACGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCt                                                                                                                                             <  6:216136/149‑1 (MQ=255)     gTACGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCa                                                                                                                                         <  1:314070/149‑1 (MQ=255)      tACGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGAACGCTGGcgc                                                                                                                                        <  7:312215/149‑3 (MQ=255)        cGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGAACGCTGGcgcg                                                                                                                                       <  1:266164/148‑4 (MQ=255)                    cgccAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAAt                                                                                                                           >  1:131983/1‑148 (MQ=255)                    cgccAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATc                                                                                                                          <  7:240691/149‑1 (MQ=255)                    cgccAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATc                                                                                                                          <  1:1444/149‑1 (MQ=255)                        acaCAGACCTACCCCTCGCCACCCCCACCGCCCCAGCACCCCTACCGCCACCATACCGCCAAGTGTCGCACCCAGCGCATAGCCACTGGACTGCCGGCGCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACAt                                                                                                                      <  4:293765/147‑1 (MQ=255)                                aTACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGctgctg                                                                                                              <  6:381184/149‑1 (MQ=255)                                    cGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGACGCGCCAGCGTTCCCTGGCATATTCCATGGCAATCACAATGctgctgcgcc                                                                                                          >  6:4883/1‑145 (MQ=255)                                     gATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAg                                                                                                         >  8:127490/1‑149 (MQ=255)                                     gATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAg                                                                                                         <  5:207719/149‑1 (MQ=255)                                       tCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGACAATCACATTGCTGCTGGCCAGGa                                                                                                       <  3:218074/149‑1 (MQ=255)                                                 cAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCCTTCGCTGGCATATTCACTGGCAATCACATAGCTGCAGGCCAGGATGCCAACg                                                                                               >  7:216085/1‑147 (MQ=255)                                                  aaTGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGAt                                                                                             <  3:267267/148‑1 (MQ=255)                                                  aaTGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGAt                                                                                             >  4:267267/1‑148 (MQ=255)                                                      ccccAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGcc                                                                                          >  8:211047/1‑147 (MQ=255)                                                               ccATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGt                                                                                 >  1:141012/1‑147 (MQ=255)                                                                 aTACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAgg                                                                             >  8:340075/1‑149 (MQ=255)                                                                     cGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCcgc                                                                         >  7:68308/1‑149 (MQ=255)                                                                             aaCCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCAAATTCACAGGCAATCACATTGCTGCTGGCCACGATGCCACCCATGCCCCGCCCGGTCCGGCCGCGCagcagc                                                                 >  6:12607/1‑149 (MQ=255)                                                                              aCCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAgg                                                                             >  7:158404/1‑136 (MQ=255)                                                                              aCCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAgg                                                                             <  8:158404/136‑1 (MQ=255)                                                                               ccGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAgcgc                                                               >  6:33341/1‑149 (MQ=255)                                                                                cGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAgcgc                                                               >  3:274374/1‑148 (MQ=255)                                                                                 gCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAgcgc                                                               >  4:57327/1‑147 (MQ=255)                                                                                  ccAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAg                                                            <  6:161107/149‑1 (MQ=255)                                                                                    aaGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGtt                                                          >  8:360288/1‑149 (MQ=255)                                                                                               ccAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGtt                                                          >  3:349597/1‑138 (MQ=255)                                                                                               ccAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGtt                                                          <  4:349597/138‑1 (MQ=255)                                                                                                  gcgcATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTgg                                            <  7:126785/149‑1 (MQ=255)                                                                                                  gcgcATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTgg                                            <  2:131983/149‑1 (MQ=255)                                                                                                      aTAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCt                                          >  8:121775/1‑147 (MQ=255)                                                                                                                   tGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGagcag                           >  5:85364/1‑149 (MQ=255)                                                                                                                     cAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGagcag                           <  6:79510/147‑1 (MQ=255)                                                                                                                     cAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGagcag                           <  7:127490/147‑1 (MQ=255)                                                                                                                        gCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGc                      <  7:362105/149‑1 (MQ=255)                                                                                                                        gCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGc                      >  3:205485/1‑149 (MQ=255)                                                                                                                          tCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGcc                     >  4:359400/1‑148 (MQ=255)                                                                                                                            aCCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCg                    <  7:211047/147‑1 (MQ=255)                                                                                                                              cGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATAc                <  4:205485/149‑1 (MQ=255)                                                                                                                                  aaGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGa            <  6:85364/149‑1 (MQ=255)                                                                                                                                    gCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACa          <  1:218245/149‑1 (MQ=255)                                                                                                                                     atatgcCAGCGAACGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAg         >  8:179146/5‑149 (MQ=255)                                                                                                                                      tatgcCAGCGAACGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCGGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAgc        >  7:204426/4‑149 (MQ=255)                                                                                                                                         gcCAGCGAACGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAgcgc      >  6:67048/1‑148 (MQ=255)                                                                                                                                         gacAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAGCGcc     <  8:216085/147‑1 (MQ=255)                                                                                                                                           cAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGAATGCCGATACCGGACAgcgc      >  3:218824/1‑146 (MQ=255)                                                                                                                                            aGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAGCGCCAGa  <  7:360288/149‑1 (MQ=255)                                                                                                                                            aGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAGCGCCAGa  <  2:141012/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                ACCAGTACGGTGACGACGCCGCCAAACAGAAATACCGATCGCCAGCCCCAATGCCCCAGCAGCCATACCGCCAGCAAACCGCCAAGGGTCGCACCCAGCGCATAGCCAGTGGACTGCAGGCTCACCGCCAAGCCGCGCCAGCGTTCGCTGGCATATTCACTGGCAATCACATTGCTGCTGGCCAGGATGCCACCGATGCCCAGCCCGGTCAGGCCGCGCAGCAGCGCCAGTTGCAGCGGGCTCTGGCTCAGTGCAGAGAGCAGCATGCCGATACCGGACAGCGCCAGA  >  NC_002947/4268741‑4269028

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