Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,757,890	G→A	8.3%	S68N (AGC→AAC)	hemBB →	porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,757,890	0	G	A	8.3%	137.7 / 2.3	48	S68N (AGC→AAC)	hemBB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base A (2/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

CTGACGAATGGACGCGTCGTCTGGTCGCTGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTGCCGCGCCGCC  >  NC_002947/3757758‑3758035                                                                                                                                     |                                                                                                                                                  cTGACGAATGGACGCGTCGTCTGGTCGCTGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCg                                                                                                                                   <  5:152032/149‑1 (MQ=255)            aCGCGTCGTCTGGTCGCTGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAgg                                                                                                                        <  5:55473/149‑1 (MQ=255)                 tcgtcTGGTCGCTGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCggg                                                                                                                   <  4:80522/149‑1 (MQ=255)                 tcgtcTGGTCGCTGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAACCTGGCCGAACTGCCGGATCa                                                                                                                              >  4:117145/1‑138 (MQ=255)                 tcgtcTGGTCGCTGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAACCTGGCCGAACTGCCGGATCa                                                                                                                              >  4:112506/1‑138 (MQ=255)                 tcgtcTGGTCGCTGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAACCTGGCCGAACTGCCGGATCa                                                                                                                              <  3:112506/138‑1 (MQ=255)                 tcgtcTGGTCGCTGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAACCTGGCCGAACTGCCGGATCa                                                                                                                              <  3:117145/138‑1 (MQ=255)                   gtcTGGTCGCTGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCg                                                                                                                 >  2:304837/1‑149 (MQ=255)                           gCTGAAAACAGCCTGTCGCGCGATGACCTGATCTCGCCGGTGTTTCTCAAGGCAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACt                                                                                                         <  5:80672/149‑1 (MQ=255)                            cTGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTg                                                                                                        >  1:264624/1‑149 (MQ=255)                             tGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAAc                                                                                                                                        >  4:268906/1‑116 (MQ=255)                             tGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAAc                                                                                                                                        <  3:268906/116‑1 (MQ=255)                                aaaCAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCa                                                                                                    >  4:80300/1‑149 (MQ=255)                                     gCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCg                                                                                               >  2:198965/1‑149 (MQ=255)                                      ccTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCgg                                                                                              >  1:133573/1‑149 (MQ=255)                                         gTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCa                                                                                           >  4:71621/1‑149 (MQ=255)                                            tcGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTg                                                                                        >  2:30358/1‑149 (MQ=255)                                              gcgATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCggg                                                                                                                   >  2:8698/1‑120 (MQ=255)                                              gcgATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCa                                                                                           <  7:254552/144‑1 (MQ=255)                                              gcgATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCa                                                                                           >  8:254552/1‑144 (MQ=255)                                              gcgATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTg                                                                                        <  1:30358/147‑1 (MQ=255)                                              gcgATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGcc                                                                                      <  3:71621/149‑1 (MQ=255)                                              gcgATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGcc                                                                                      <  4:103931/149‑1 (MQ=255)                                                              gCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTc                                                                      >  8:110350/1‑149 (MQ=255)                                                                  gTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGt                                                                  <  3:80300/149‑1 (MQ=255)                                                                      ttCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCa                                                               <  2:237774/148‑1 (MQ=255)                                                                         tcAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGTACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAgcg                                                            >  4:221703/1‑148 (MQ=255)                                                                         tcAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAgcg                                                            >  6:5324/1‑148 (MQ=255)                                                                          cAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTc                                                                                                 >  6:71071/1‑110 (MQ=255)                                                                          cAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTc                                                                                                 <  5:71071/110‑1 (MQ=255)                                                                            aggaaggCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCaaa                                                        <  7:110350/149‑1 (MQ=255)                                                                              gaaggCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGa                                                      >  3:318867/1‑149 (MQ=255)                                                                                aggCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCAAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACa                                                    <  6:109056/149‑1 (MQ=255)                                                                                    tCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGt                                                  >  6:274815/1‑147 (MQ=255)                                                                                      aGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGa                                                <  1:198965/147‑1 (MQ=255)                                                                                        gTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAggg                                             >  4:60661/1‑148 (MQ=255)                                                                                                    cGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGTACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTc                                <  3:221703/149‑1 (MQ=255)                                                                                                    cGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTc                                <  3:146550/149‑1 (MQ=255)                                                                                                     gATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTc                                >  7:219652/1‑148 (MQ=255)                                                                                                         aCCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGAcg                            <  3:60661/148‑1 (MQ=255)                                                                                                          ccGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGAcgcg                          >  6:161343/1‑149 (MQ=255)                                                                                                             ccATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCTGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGa                        >  3:46321/1‑148 (MQ=255)                                                                                                                    gggCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTACGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGa                >  5:98213/1‑149 (MQ=255)                                                                                                                       ccAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCt             >  2:125385/1‑149 (MQ=255)                                                                                                                          ggTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCTGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTacc          <  4:46321/149‑1 (MQ=255)                                                                                                                           gTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTgccg         >  4:57348/1‑149 (MQ=255)                                                                                                                                 gTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTgccgcgccgc   <  5:5324/149‑1 (MQ=255)                                                                                                                                  tCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTGCCGCgccgcc  <  3:57348/149‑1 (MQ=255)                                                                                                                                     |                                                                                                                                                  CTGACGAATGGACGCGTCGTCTGGTCGCTGAAAACAGCCTGTCTCGCGATGACCTGATCTGGCCGGTTTTTCTCAAGGAAGGCTCAGGTACACGCGAGGCGATCACCGCCATGCCGGGCCAGGTACGCGTCAGCCTGGCCGAACTGCCGGATCACGTAGGCCGGGCGGTCGAACTGGGCATTCCGGCCATTGCCCTGTTTCCCGCCACTCCAGTGCAGCGCAAAGACAGTGAGGGCACCCACGCCCTCGACGCGGACAACCTGATCTGCCGCGCCGCC  >  NC_002947/3757758‑3758035

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