Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	799,844	T→A	15.6%	D211V (GAC→GTC)	ispB ←	octaprenyl diphosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	799,844	0	T	A	15.6%	110.5 / 13.7	45	D211V (GAC→GTC)	ispB	octaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (19/19);  new base A (5/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

GATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCAT  >  NC_002947/799699‑799986                                                                                                                                                  |                                                                                                                                               gATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGAc                                                                                                                                            <  1:673249/150‑1 (MQ=255)       cTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCTGGTCGACGACctgctg                                                                                                                                      <  1:35362/150‑1 (MQ=255)              aCCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGcc                                                                                                                               <  2:34379/150‑1 (MQ=255)                cAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGgtcgt                                                                                                                                                <  2:389701/131‑1 (MQ=255)                cAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGgtcgt                                                                                                                                                >  1:389701/1‑131 (MQ=255)                                           gcgcATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAg                                                                                                  <  1:518356/150‑1 (MQ=255)                                                 ggTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTc                                                                                            <  2:57739/150‑1 (MQ=255)                                                  gtgtAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTcg                                                                                           >  1:455460/1‑150 (MQ=255)                                                     tAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCgcgc                                                                                        >  1:523064/1‑150 (MQ=255)                                                      aGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCgct                                                                                       <  1:482979/150‑1 (MQ=255)                                                         tCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGCTCTGCGAGTCACCCTTGTAGTCCAGCAGGCCGTCGACCAGCTGGAAGGCCACGCCACGATGGTCGCCAAAGGTACGCAGGGCTTCGCgctgct                                                                                    >  2:87821/1‑150 (MQ=255)                                                                   ggTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGACCAGAGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTGGCGCGGGTCGTggggggg                                                                          >  2:632241/1‑149 (MQ=255)                                                                         tttGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACTCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGc                                                                    >  1:31314/1‑150 (MQ=255)                                                                         tttGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGc                                                                    >  2:844868/1‑150 (MQ=255)                                                                              ccTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAgcg                                                               >  1:753227/1‑150 (MQ=255)                                                                                      aGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCgcagca                                                       >  1:717562/1‑150 (MQ=255)                                                                                           gtcACCGACGTTCTTGCCCAGGGTCTGGGTCTCCATGTTGTAGTCCAGCAGGCCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgtg                                                  <  2:566558/150‑1 (MQ=255)                                                                                             cACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTggg                                                <  1:844868/150‑1 (MQ=255)                                                                                              aCCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTg                                                                       >  1:726975/1‑126 (MQ=255)                                                                                              aCCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTg                                                                       <  2:726975/126‑1 (MQ=255)                                                                                                      tCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCtt                                       <  1:632241/150‑1 (MQ=255)                                                                                                       cTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                      >  2:505230/1‑150 (MQ=255)                                                                                                            ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                      <  1:421626/145‑1 (MQ=255)                                                                                                            ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                      >  2:421626/1‑145 (MQ=255)                                                                                                             cAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAg                                >  1:421343/1‑150 (MQ=255)                                                                                                              aGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCACAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGc                               <  2:421343/150‑1 (MQ=255)                                                                                                                    tGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGcc                         <  2:572482/150‑1 (MQ=255)                                                                                                                      cGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCGTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGt                       >  1:401684/1‑150 (MQ=255)                                                                                                                          tCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTg                   <  2:455460/150‑1 (MQ=255)                                                                                                                            aCCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGcgccg                                                            >  1:158865/1‑107 (MQ=255)                                                                                                                            aCCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGcgccg                                                            <  2:158865/107‑1 (MQ=255)                                                                                                                            aCCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGcc                 >  1:206671/1‑150 (MQ=255)                                                                                                                                ttGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGc                  <  1:709360/145‑1 (MQ=255)                                                                                                                                ttGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGc                  >  2:709360/1‑145 (MQ=255)                                                                                                                                 tGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGaa            <  2:206671/150‑1 (MQ=255)                                                                                                                                  gccttCCAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAAt           >  1:264841/5‑150 (MQ=255)                                                                                                                                   tAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGaa            >  2:360012/1‑148 (MQ=255)                                                                                                                                   tAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGaa            <  1:360012/148‑1 (MQ=255)                                                                                                                                    cttCCAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGa         >  1:168759/3‑150 (MQ=255)                                                                                                                                     ttCCAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGAc        >  2:162205/2‑150 (MQ=255)                                                                                                                                     gTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAGGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGAc        <  1:87821/150‑1 (MQ=255)                                                                                                                                      tCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGAcc       >  1:2443/1‑150 (MQ=255)                                                                                                                                        cagctgGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTggg                                                <  2:660793/107‑1 (MQ=255)                                                                                                                                        cagctgGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTggg                                                >  1:660793/1‑107 (MQ=255)                                                                                                                                        cagctgGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTc     >  2:740941/1‑150 (MQ=255)                                                                                                                                        cagcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGGACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGGCTTGCCGCGAATGACGTc     >  2:844551/1‑150 (MQ=255)                                                                                                                                           cagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCAt  <  2:2443/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                               GATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCAT  >  NC_002947/799699‑799986

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