Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	151,912	G→A	9.5%	T65T (ACC→ACT)	PP_0143 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	151,912	0	G	A	9.5%	112.5 / 2.5	42	T65T (ACC→ACT)	PP_0143	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (16/22);  new base A (2/2);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

CGGTGAAGACCAGCTGACCGCTGTTCCAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGAGAGGAAACCGCC  >  NC_002947/151764‑152034                                                                                                                                                     |                                                                                                                           cGGTGAAGACCAGCTGACCGCTGTTCAAGGGGGGGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTACTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGAGGTGGAAAGTGCGAGCGCAAATCCCg                                                                                                                            >  6:78504/1‑149 (MQ=255)   gTGAAGACCAGCTGACCGCTGTTCCAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCgg                                                                                                                           <  5:344166/148‑1 (MQ=255)   gTGAAGACCAGCTGACCGCTGTTCCAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCgg                                                                                                                           >  6:344166/1‑148 (MQ=255)       aGACCAGCTGACCGCTGTTCCAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAgg                                                                                                                      <  8:181434/149‑1 (MQ=255)        gACCAGCTGACCGCTGTTCCAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAggg                                                                                                                     >  3:339636/1‑149 (MQ=255)             gCTGACCGCTGTTCCAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTtc                                                                                                                <  5:35157/149‑1 (MQ=255)                 aCCGCTGTTCCAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCa                                                                                                            <  5:314736/149‑1 (MQ=255)                   cGCTGTTCCAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAg                                                                                                           >  2:156150/1‑148 (MQ=255)                     cTGTTCCAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGgtg                                                                                                                   >  6:170372/1‑138 (MQ=255)                     cTGTTCCAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGgtg                                                                                                                   <  5:170372/138‑1 (MQ=255)                          ccAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGc                                                                                                    >  6:244712/1‑148 (MQ=255)                             ggtggtggtGATGAAAAAGAACGGCAGCGCGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGa                                                                                                >  3:351225/1‑149 (MQ=255)                                ggtggtgATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGAcc                                                                                              >  7:87109/1‑148 (MQ=255)                                   ggtgATGAAACAGAACGGCCGCACGAAGAACAGCGTTTCCTGGTTGATCATCTGCGTTGCGTATCGCAGCAGCGGATGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCATCATCAGCCAGACCCAGc                                                                                          <  5:78504/149‑1 (MQ=255)                                         gaAAAAGAACGGCAGCACGAAGAACAGCGTTTCATGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTCTTCTCCAGCATCAGCCAGACCCAGCTGACCa                                                                                    >  1:41916/1‑149 (MQ=255)                                         gaAAAAGAACGGCAGCAAGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCAGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGAc                                                                                      >  1:47181/1‑147 (MQ=255)                                            aaaGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCa                                                                                 >  4:25094/1‑149 (MQ=255)                                                 aCGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGc                                                                            <  3:235681/149‑1 (MQ=255)                                                     cagcaCGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACgg                                                                        <  1:156150/149‑1 (MQ=255)                                                               aaCAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCAGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGa                                                                     <  3:295787/142‑1 (MQ=255)                                                               aaCAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCAGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGa                                                                     >  4:295787/1‑142 (MQ=255)                                                                     gTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTggccag                                                        <  3:25094/149‑1 (MQ=255)                                                                       ttCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATc                                                                                                       >  4:304608/1‑100 (MQ=255)                                                                       ttCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATc                                                                                                       <  3:304608/100‑1 (MQ=255)                                                                         ccTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGTTGACCAGCATGAGCACGGCGACCCAGCTggccaggc                                                      >  2:281979/1‑147 (MQ=255)                                                                             gtggATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGGGCACGGCGACCCAGCTGGCCGGGCCCGcgg                                                >  4:241913/1‑147 (MQ=255)                                                                                gATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGctggcctggc                                             <  5:244712/149‑1 (MQ=255)                                                                                 atcatcTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCg                                            <  4:307132/149‑1 (MQ=255)                                                                                    atcTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTc                                         <  5:277019/149‑1 (MQ=255)                                                                                     tcTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCAGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGctggcctggccag                                          <  8:129135/147‑4 (MQ=255)                                                                                     tcTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCAGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGctggcctggccag                                          >  7:129135/1‑144 (MQ=255)                                                                                         cgttgcgtAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCa                                    <  4:351225/149‑1 (MQ=255)                                                                                              cgtAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGAt                                >  3:54894/1‑148 (MQ=255)                                                                                                     cagcagCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGc                        >  2:153458/1‑149 (MQ=255)                                                                                                     cagcagCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGc                        <  1:153458/149‑1 (MQ=255)                                                                                                     cagcagCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGc                        <  3:241913/149‑1 (MQ=255)                                                                                                      agcagcGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGc                                                  <  5:195520/122‑1 (MQ=255)                                                                                                      agcagcGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGc                                                  >  6:195520/1‑122 (MQ=255)                                                                                                       gcagcGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCg                      <  4:339636/149‑1 (MQ=255)                                                                                                                 aGGAAGCTCGCGGTTGAAAGTGCGCGCGAACCGCCCGGCGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGgaga            <  4:250866/149‑1 (MQ=255)                                                                                                                  ggAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGgaga            <  4:54894/148‑1 (MQ=255)                                                                                                                  ggAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGagag           <  8:87109/149‑1 (MQ=255)                                                                                                                    aaTCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGAGAgg          >  2:157088/1‑148 (MQ=255)                                                                                                                           cgGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCTACCAGGATGAAGCTGCCGATACCGGAGAGGAAAccgcc  <  2:41916/149‑1 (MQ=255)                                                                                                                           cgGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGAGAGGTAAccgcc  <  2:47181/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                           CGGTGAAGACCAGCTGACCGCTGTTCCAGGTGGTGGTGATGAAAAAGAACGGCAGCACGAAGAACAGCGTTTCCTGGTGGATCATCTGCGTTGCGTAGCGCAGCAGCGGCTGAGGAATCTCGCGGTTGAAAGTGCGCGCGAACAGCCCGGTGAGGGTGTTCTCCAGCATCAGCCAGACCCAGCTGACCAGCATGAGCACGGCGACCCAGCTGGCCAGGCCAGCCTGGCGGTCGACCAGGATGAAGCTGCCGATACCGGAGAGGAAACCGCC  >  NC_002947/151764‑152034

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