Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	799,835	A→T	14.1%	L213Q (CTG→CAG)	ispB ←	octaprenyl diphosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	799,835	0	A	T	14.1%	98.2 / 6.9	41	L213Q (CTG→CAG)	ispB	octaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (16/19);  new base T (2/4);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCAT  >  NC_002947/799696‑799983                                                                                                                                            |                                                                                                                                                     gATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGgtcgt                                                                                                                                                <  8:204022/146‑1 (MQ=255)  aTCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGTCAGCGTAGTATAAAAAAAAAAAAAATCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGAc                                                                                                                                            >  2:236995/1‑149 (MQ=255)          gCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCTGGTCGACGACctgctg                                                                                                                                      <  1:247462/147‑1 (MQ=255)           cGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGaa                                                                                                                                   <  1:155568/149‑1 (MQ=255)                          tGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATg                                                                                                                    >  5:101166/1‑149 (MQ=255)                               ggCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCg                                                                                                                <  2:57700/148‑1 (MQ=255)                                   gTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCaaa                                                                                                           <  4:7846/149‑1 (MQ=255)                                      ccTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGt                                                                                                        <  8:116184/149‑1 (MQ=255)                                        ttCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGaa                                                                                                                                   <  6:8493/120‑1 (MQ=255)                                        ttCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGaa                                                                                                                                   >  5:8493/1‑120 (MQ=255)                                                  gtgtTGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGGCACCGACGTTCTTGGTCTGGGGCTGGGTGTTACCCTTGTCGTCCAGCAGGGCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTc                                                                                            <  1:257388/149‑1 (MQ=255)                                                  gtgtAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGGCGTCACCGACGTGCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCtt                                                                                             <  6:190881/148‑1 (MQ=255)                                                                cAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCgg                                                                              <  6:101166/149‑1 (MQ=255)                                                                    gTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGc                                                                          >  4:232148/1‑149 (MQ=255)                                                                        gTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGt                                                                                                                  >  2:287259/1‑105 (MQ=255)                                                                        gTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGt                                                                                                                  <  1:287259/105‑1 (MQ=255)                                                                         tttGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTggg                                                                     >  2:154438/1‑149 (MQ=255)                                                                             cccTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCa                                                                  <  3:232148/148‑1 (MQ=255)                                                                              ccTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCTGCTGGTGGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAgc                                                                <  1:236995/149‑1 (MQ=255)                                                                                         tcgtcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACt                                                     <  1:51265/149‑1 (MQ=255)                                                                                           gtcACCGACGTTCTTTCCCAGTGTGTGCGAGTCACCCTTTTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTAGGCAGGGCTTCGCGCTGCTCGTCGGTGGCGTGGGCCAGCGCCGCAGCACtgt                                                   <  4:116766/149‑1 (MQ=255)                                                                                            tcACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgtg                                                  >  2:175642/1‑149 (MQ=255)                                                                                              aCCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTgg                                                 <  3:115244/148‑1 (MQ=255)                                                                                                cGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGc                                                                                                                           <  6:14811/72‑1 (MQ=255)                                                                                                cGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGc                                                                                                                           >  5:14811/1‑72 (MQ=255)                                                                                                cGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTggg                                                >  4:271533/1‑147 (MQ=255)                                                                                                 gACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgt                                                   <  4:221069/143‑1 (MQ=255)                                                                                                 gACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACtgt                                                   >  3:221069/1‑143 (MQ=255)                                                                                                            ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTc                                                                                                                 <  1:283888/70‑1 (MQ=255)                                                                                                            ccAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTc                                                                                                                 >  2:283888/1‑70 (MQ=255)                                                                                                                ggTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCa                              >  3:135709/1‑149 (MQ=255)                                                                                                                            aCCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGc                  >  6:100648/1‑149 (MQ=255)                                                                                                                             cccTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGcc                 <  1:175642/149‑1 (MQ=255)                                                                                                                              gcTTGTAGTCCATCAGGTCGTCGATCAGCTGGAAGTCCACGCCCAGATGGTCGCCAAATGTACGCAGGGCTTCGCGCTGCTCGTCGGtt                                                                            <  5:63375/88‑2 (MQ=255)                                                                                                                              ccTTGTAGTCCAGCAGGTCTTCGATCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAACGTACGCAGGGCTTCGCGCAGCTCGTCGGTg                                                                            >  6:63375/1‑89 (MQ=255)                                                                                                                              ccTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCcg                >  4:44771/1‑149 (MQ=255)                                                                                                                                    cttCCAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                      >  8:113244/3‑121 (MQ=255)                                                                                                                                    cttCCAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTc                                      <  7:113244/119‑1 (MQ=255)                                                                                                                                       ccAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACa                                 >  2:207930/1‑123 (MQ=255)                                                                                                                                       ccAGCTGGTCGACGACCTGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACa                                 <  1:207930/123‑1 (MQ=255)                                                                                                                                       ccAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGAcc       >  4:121734/1‑149 (MQ=255)                                                                                                                                        cagcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCt      >  5:322082/1‑149 (MQ=255)                                                                                                                                          gcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTcc    <  3:121734/149‑1 (MQ=255)                                                                                                                                          gcagGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTc     >  8:271024/1‑148 (MQ=255)                                                                                                                                            agGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCAt  <  5:100648/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                     GATCGCCTGGCGTACCAGCGCCGCCTGCTCGGCAGTGCCTTCGCGCATGGTGTAGATCAGCGGCAGGGTAGGTTTGCCCTCGGCCAGGTCGTCACCGACGTTCTTGCCCAGGGTCTGCGAGTCACCCTTGTAGTCCAGCAGGTCGTCGACCAGCTGGAAGGCCACGCCCAGATGGTCGCCAAAGGTACGCAGGGCTTCGCGCTGCTCGTCGGTGGCCTGGGCCAGCGCCGCAGCACTGTGGGTAGAAGCTTCGAACAGCATCGCCGTCTTGCCGCGAATGACCTCCAT  >  NC_002947/799696‑799983

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