Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,686,441	A→G	12.4%	Q283Q (CAA→CAG)	PP_3251 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,686,441	0	A	G	12.4%	99.3 / 5.1	40	Q283Q (CAA→CAG)	PP_3251	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/18);  new base G (2/3);  total (19/21)
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.00e+00

GGCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCC  >  NC_002947/3686295‑3686561                                                                                                                                                   |                                                                                                                         ggCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAgg                                                                                                                        <  7:68424/149‑1 (MQ=255)     aCGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGa                                                                                                                     >  2:262038/1‑148 (MQ=255)       gcccgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACt                                                                                                                  >  1:18927/1‑149 (MQ=255)          cgccTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGCCTTGCAGGTCGAACTggc                                                                                                               <  5:34320/149‑3 (MQ=255)                gCAATTCAGGCGCTGAACCAGGTAAGCCGCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCa                                                                                                           <  7:216914/147‑1 (MQ=255)                gCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGCCTTGCAGGTCGAACTg                                                                                                                 <  5:201656/141‑1 (MQ=255)                gCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGCCTTGCAGGTCGAACTg                                                                                                                 >  6:201656/1‑141 (MQ=255)                            cTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGc                                                                                               >  3:37374/1‑147 (MQ=255)                               aaCCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCa                                                                                                           <  2:245806/132‑1 (MQ=255)                               aaCCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCa                                                                                                           >  1:245806/1‑132 (MQ=255)                                     cTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGGCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGCCTTGCAGGTCGAACTg                                                                                                                 <  8:8861/120‑1 (MQ=255)                                     cTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGCCTTGCAGGTCGAACTg                                                                                                                 >  7:8861/1‑120 (MQ=255)                                      tAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTAt                                                                                   <  1:262038/149‑1 (MQ=255)                                      tAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTAt                                                                                   >  4:2679/1‑149 (MQ=255)                                                gcATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGa                                                                         >  2:33273/1‑149 (MQ=255)                                                gcATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACATCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGa                                                                         >  5:71164/1‑149 (MQ=255)                                                 cATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGat                                                                        <  3:175544/149‑1 (MQ=255)                                                          cAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGa                                                               <  1:33273/149‑1 (MQ=255)                                                                ttCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTg                                                         <  6:71164/149‑1 (MQ=255)                                                                   tACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCa                                                                                                           <  4:227069/96‑1 (MQ=255)                                                                   tACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCa                                                                                                           >  3:227069/1‑96 (MQ=255)                                                                      gCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCt                                                   <  3:2679/149‑1 (MQ=255)                                                                      gCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCt                                                   >  6:96808/1‑149 (MQ=255)                                                                           ggCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTgcag                                               >  1:135642/1‑148 (MQ=255)                                                                              gcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCagcagc                                           <  4:37374/149‑1 (MQ=255)                                                                              gcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCagcagc                                           <  2:135642/149‑1 (MQ=255)                                                                              gcgcTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCagcagc                                           <  2:74295/149‑1 (MQ=255)                                                                                     ggCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGa                                     <  2:18927/148‑1 (MQ=255)                                                                                                   tGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAAtggctgg                      <  8:159143/149‑1 (MQ=255)                                                                                                              tCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTg                                                  <  7:1220/110‑1 (MQ=255)                                                                                                              tCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTg                                                  >  8:1220/1‑110 (MQ=255)                                                                                                                gccacgccaGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGgcgc         >  6:242592/1‑149 (MQ=255)                                                                                                                 ccacgccaGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCt                                                   >  6:131466/1‑106 (MQ=255)                                                                                                                 ccacgccaGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCt                                                   <  5:131466/106‑1 (MQ=255)                                                                                                                      ccaGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGcc   >  1:56009/1‑149 (MQ=255)                                                                                                                       caGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGccc  >  2:101891/1‑149 (MQ=255)                                                                                                                                 ccGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTgc                                                        <  3:30201/85‑1 (MQ=255)                                                                                                                                 ccGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTgc                                                        >  4:30201/1‑85 (MQ=255)                                                                                                                                                gCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGcc        >  2:140826/1‑118 (MQ=255)                                                                                                                                                gCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGcc        <  1:140826/118‑1 (MQ=255)                                                                                                                                                   |                                                                                                                         GGCCACGCCCGCCTGGCAATTCAGGCGCTGAACCAGCTAAGCCCCGCGCATGATCCGCAGGCGTTCTACGCCAGGGCGCGCTTAGGCTACCGCCTGAATGACCGGGGTATCGCCACGCCAGACCTCATCCGCCAGTTCGACCTGCAAGGCGAACTGTTTGCAGCCCTTGAACGCAAGCGGGTCTATGGCCAATTGATGCATGCCGACCATTGCCGCCTGCAGCAGCGCACGATCAATCAATGGCTGGCCGAACCCGGCGCCATGCCC  >  NC_002947/3686295‑3686561

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