Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,591,577	T→A	12.4%	L258* (TTG→TAG) ‡	PP_3169 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,591,577	0	T	A	12.4%	97.1 / 4.7	39	L258* (TTG→TAG) ‡	PP_3169	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (16/18);  new base A (3/2);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGGTGT  >  NC_002947/3591436‑3591724                                                                                                                                              |                                                                                                                                                    gggTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTgg                                                                                                                                              <  5:82166/149‑1 (MQ=255) gggTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTgg                                                                                                                                              <  4:108819/149‑1 (MQ=255)       tGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCAcc                                                                                                                                        >  7:101744/1‑149 (MQ=255)           cAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTg                                                                                                                                               <  3:218114/138‑1 (MQ=255)           cAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTg                                                                                                                                               >  4:218114/1‑138 (MQ=255)                   cGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCAc                                                                                                                                         >  6:263418/1‑136 (MQ=255)                         gCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTAt                                                                                                                      <  7:112275/149‑1 (MQ=255)                         gCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTAt                                                                                                                      >  4:147895/1‑149 (MQ=255)                             ttACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGAcc                                                                                                                  <  6:189027/149‑1 (MQ=255)                                 tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGCACTGGCGCAAGCACTCAACGATGACCCGTTGCACCAGCGCTAGCGCTGGTTACCCGATGACTGCCTGCGGTAAGACCAGAt                                                                                                              >  5:138034/1‑149 (MQ=255)                                 tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAt                                                                                                              >  7:86237/1‑149 (MQ=255)                                 tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGt                                                                                                                                             <  7:193808/118‑1 (MQ=255)                                 tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGt                                                                                                                                             >  8:193808/1‑118 (MQ=255)                                     aaGGCATCATCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGc                                                                                                          <  4:66845/149‑1 (MQ=255)                                      aGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGc                                                                                                          >  6:35195/1‑148 (MQ=255)                                      aGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCg                                                                                                         <  3:147895/149‑1 (MQ=255)                                         cATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGt                                                                                                                                             <  4:51622/110‑1 (MQ=255)                                         cATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGt                                                                                                                                             >  3:51622/1‑110 (MQ=255)                                             agcagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACt                                                                                                  <  5:35195/149‑1 (MQ=255)                                                 gcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggtg                                                                                              <  5:270086/149‑1 (MQ=255)                                                      cTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCa                                                                                          <  8:86237/148‑1 (MQ=255)                                                       tCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCAcc                                                                                        >  1:270899/1‑149 (MQ=255)                                                          aCCAGCTACCACCCGAAATGCAACGCTTGTCGCGCTACCTGGAACTGGCGCAAGCACTGATCGATGACCCGGTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGAc                                                                                     <  6:138034/149‑1 (MQ=255)                                                             aGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTg                                                                                  <  4:79549/149‑1 (MQ=255)                                                              gCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTgc                                                                                 >  1:115406/1‑149 (MQ=255)                                                                      ccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgggg                                                                         >  8:63329/1‑149 (MQ=255)                                                                      ccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATggg                                                                          >  2:151407/1‑148 (MQ=255)                                                                                           gcAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCttt                                                    <  6:112022/149‑1 (MQ=255)                                                                                              aCCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGc                                                 <  1:15130/149‑1 (MQ=255)                                                                                                 tGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCtt                                               >  7:239125/1‑148 (MQ=255)                                                                                                   gAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCa                                             >  8:179770/1‑148 (MQ=255)                                                                                                         gcgcAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATgg                                       <  8:101744/148‑1 (MQ=255)                                                                                                                 aCTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAgcgc                               >  6:97502/1‑148 (MQ=255)                                                                                                                     aaCGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGatc                           <  2:270899/148‑1 (MQ=255)                                                                                                                      aCGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGatcat                         >  5:183741/1‑149 (MQ=255)                                                                                                                                      ccAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATgg         >  8:23072/1‑149 (MQ=255)                                                                                                                                         gcgcTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCg      >  8:152011/1‑149 (MQ=255)                                                                                                                                         gcgcTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCg      <  8:239125/149‑1 (MQ=255)                                                                                                                                           gcTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGgt    <  5:97502/149‑1 (MQ=255)                                                                                                                                             ttGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGgtgt  >  8:24227/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                    GGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGTCGGTGT  >  NC_002947/3591436‑3591724

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