Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	863,379	A→T	8.1%	A170A (GCA→GCT)	hemH →	ferrochelatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	863,379	0	A	T	8.1%	148.6 / 3.6	49	A170A (GCA→GCT)	hemH	ferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (18/27);  new base T (2/2);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGT  >  NC_002947/863239‑863517                                                                                                                                             |                                                                                                                                           cTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCTAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGctgccagttct                                                                                                                                   <  2:1243634/150‑5 (MQ=255)          caGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTAc                                                                                                                          <  1:1171753/150‑1 (MQ=255)             ttCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTg                                                                                                                       <  2:455311/150‑1 (MQ=255)                    aCAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAgg                                                                                                                <  2:517710/150‑1 (MQ=255)                    aCAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAgg                                                                                                                <  2:261972/150‑1 (MQ=255)                    aCAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAgg                                                                                                                <  2:247280/150‑1 (MQ=255)                    aCAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACACCCCGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAgg                                                                                                                <  2:916339/150‑1 (MQ=255)                           ggTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGa                                                                                                         >  2:1200332/1‑150 (MQ=255)                             tGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGacc                                                                                                       >  2:816396/1‑150 (MQ=255)                                  aCTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTg                                                                                                  >  2:254582/1‑150 (MQ=255)                                    tgtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCt                                                                                                >  1:838654/1‑150 (MQ=255)                                     gtgATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTa                                                                                               <  1:816396/150‑1 (MQ=255)                                                   aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATgg                                                                                 >  1:344741/1‑150 (MQ=255)                                                   aaaaCAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATgg                                                                                 >  1:140693/1‑150 (MQ=255)                                                         gACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGcc                                                                           <  1:1200332/150‑1 (MQ=255)                                                             gTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTAGCAGCCAGCGCCCGGCCGAAGCTTAGGCGGGATTACGATGACCTGCTACTGAGCTTCCATGGTTTGCCGGaa                                                                       <  1:958418/150‑1 (MQ=255)                                                             gTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGAGCTATGAGAACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGaa                                                                       >  2:1119096/1‑150 (MQ=255)                                                              tCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAAc                                                                      <  2:275486/150‑1 (MQ=255)                                                                 gCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCt                                                                                        <  2:478566/129‑1 (MQ=255)                                                                 gCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCt                                                                                        >  1:478566/1‑129 (MQ=255)                                                                   gAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCAc                                                                 <  1:729785/150‑1 (MQ=255)                                                                      cACCAGTTACCGTTGCAGATGCGCGTGCTGCGGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTg                                                              <  2:492351/150‑1 (MQ=255)                                                                           gTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTgaagaa                                                         >  2:1141434/1‑150 (MQ=255)                                                                                 gTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGtt                                                   <  2:838654/150‑1 (MQ=255)                                                                                  ttGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGttt                                                  <  2:140693/150‑1 (MQ=255)                                                                                   tGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTc                                                 <  1:845139/150‑1 (MQ=255)                                                                                           cgcgTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAgggg                                         <  1:254582/150‑1 (MQ=255)                                                                                              gTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGa                                                                                                         >  2:807060/1‑83 (MQ=255)                                                                                              gTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGa                                                                                                         <  1:807060/83‑1 (MQ=255)                                                                                                 ctgcAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCa                                                                  >  1:90421/1‑119 (MQ=255)                                                                                                 ctgcAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCa                                                                  <  2:90421/119‑1 (MQ=255)                                                                                                      gCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTa                                                                                               >  1:365306/1‑85 (MQ=255)                                                                                                      gCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTa                                                                                               <  2:365306/85‑1 (MQ=255)                                                                                                          ttCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTc                                      <  1:329243/138‑1 (MQ=255)                                                                                                          ttCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTc                                      >  2:329243/1‑138 (MQ=255)                                                                                                             tATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgc                       >  2:1427278/1‑150 (MQ=255)                                                                                                                   cACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGcc                 <  2:889534/150‑1 (MQ=255)                                                                                                                      ccGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAt              >  2:205195/1‑150 (MQ=255)                                                                                                                      ccGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGAt              <  1:1427278/150‑1 (MQ=255)                                                                                                                        ggCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACTGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtg            <  2:344741/150‑1 (MQ=255)                                                                                                                          ggTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATtgct          <  2:554237/149‑1 (MQ=255)                                                                                                                               aTCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCAc     >  1:1184114/1‑150 (MQ=255)                                                                                                                                tCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACg    >  1:1039073/1‑150 (MQ=255)                                                                                                                                  gAGGCGCGTGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGt  <  2:411760/150‑1 (MQ=255)                                                                                                                                   gggCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGttt                                                                              <  1:1260299/72‑1 (MQ=255)                                                                                                                                   gggCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGttt                                                                              >  2:1260299/2‑73 (MQ=255)                                                                                                                                   gggCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgcg                      <  1:1277822/128‑1 (MQ=255)                                                                                                                                   gggCGCTGGCTGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTgcgcg                      >  2:1277822/2‑129 (MQ=255)                                                                                                                                           gCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGc                  >  2:455732/1‑125 (MQ=255)                                                                                                                                           gCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGc                  <  1:455732/125‑1 (MQ=255)                                                                                                                                             |                                                                                                                                           CTTTATCCACAGTTCGCCGACAGCACGGTGACCACTGTGATCGAGCAGGCAAAACAGACGGTCAGCGAGCACCAGTTACCGTTGCAGATGCGCGTGCTGCAGCCGTTCTATGAGCACCCGGCGTACATCGAGGCGCTGGCAGCCAGCGCCCGGCCGTACCTGGAGCAGGGTTACGACCACCTGCTACTGAGCTTCCATGGTTTGCCGGAACGGCACCTGAAGAAACTGTTTCCGAAGGGGGTCAAGCACGACCTGCGCGCTGCCGATTGCTGCCACGGT  >  NC_002947/863239‑863517

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑