Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,147,968	G→T	8.2%	P34Q (CCG→CAG)	PP_1006 ←	heme receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,147,968	0	G	T	8.2%	148.5 / 3.6	49	P34Q (CCG→CAG)	PP_1006	heme receptor
Reads supporting (aligned to +/- strand):  ref base G (16/29);  new base T (2/2);  total (18/31)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

AAGCCGACCTGCCAGCCGGTGACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCTGTTCTC  >  NC_002947/1147839‑1148116                                                                                                                                  |                                                                                                                                                     aaGCCGACCTGCCAGCCGGTGACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTg                                                                                                                                  <  2:902069/150‑1 (MQ=255)               gCCGGTGACCTGGCTGAAGGCGTTGATGGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCACGCCCGGGCTGCTCGCTTGGCCTGGCGGGGTGGCCAGCGCCGTCAGGCTGCATTCGCCGCTGGCGAGCAGGGCAAg                                                                                                                    <  2:137474/150‑1 (MQ=255)                ccGGTGACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAgg                                                                                                                   <  2:1119454/150‑1 (MQ=255)                ccGGTGACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAgg                                                                                                                   <  1:940561/150‑1 (MQ=255)                     gACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCaa                                                                                                              <  1:931311/150‑1 (MQ=255)                      aCCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCTGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAAc                                                                                                             <  2:466301/150‑1 (MQ=255)                        cTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAg                                                                                                           >  1:927454/1‑150 (MQ=255)                        cTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAg                                                                                                           <  1:810934/150‑1 (MQ=255)                            cTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGAc                                                                                                       <  2:364018/150‑1 (MQ=255)                               aaGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGa                                                                                                                              <  2:273448/124‑1 (MQ=255)                               aaGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGa                                                                                                                              >  1:273448/1‑124 (MQ=255)                                  gCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCTGCAGCCTGccg                                                                                                                                           >  1:632848/1‑106 (MQ=255)                                  gCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCTGCAGCCTGccg                                                                                                                                           <  2:632848/108‑3 (MQ=255)                                   cGTTGATCGACGACACCACCGGCTTCTGAACGATGGTGAAGAGGTAGTCGCCCATTCGCGGGCTACTGGCCTGAGGCGGCTCGTCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTc                                                                                                <  2:73741/150‑1 (MQ=255)                                     ttGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCt                                                                                              <  1:887514/150‑1 (MQ=255)                                        aTCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGGCAACAGGGACAGTTGTCCTGTa                                                                                           <  2:418183/150‑1 (MQ=255)                                                 acCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGa                                                                                  >  1:330706/1‑150 (MQ=255)                                                  ccAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCCGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGaa                                                                                 <  1:378791/150‑1 (MQ=255)                                                   cAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAg                                                                                <  1:614969/150‑1 (MQ=255)                                                   cAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAg                                                                                >  2:1110838/1‑150 (MQ=255)                                                     gCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGGCAACAGGGGCAGTTGTCCTGTACGGCgggtcaggg                                                                              >  2:142470/1‑144 (MQ=255)                                                       ggCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGATTGGCGCGGCCCGGCAGCCGCCTCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTg                                                                            <  1:142470/150‑1 (MQ=255)                                                            ctgGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGCGCAAGGCAAGGGTCTAAAAGGCAAGTTGGCCGGGAGGGCGGGGAAGGTCGGAcca                                                                       >  2:197464/1‑148 (MQ=255)                                                                 cGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAg                                                                                                                          <  1:846942/94‑1 (MQ=255)                                                                 cGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAg                                                                                                                          >  2:846942/1‑94 (MQ=255)                                                                 cGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCggg                                                                  <  1:160818/150‑1 (MQ=255)                                                                   aTGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTa                                                                <  1:9347/150‑1 (MQ=255)                                                                      cTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGa                                                             >  2:1077201/1‑150 (MQ=255)                                                                       tGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAgg                                                                               >  1:630436/1‑131 (MQ=255)                                                                       tGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAgg                                                                               <  2:630436/131‑1 (MQ=255)                                                                        tcaGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAg                                                               >  1:948029/3‑146 (MQ=255)                                                                        gAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAg                                                               <  2:948029/146‑1 (MQ=255)                                                                              gTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGa                                                     >  2:1084562/1‑150 (MQ=255)                                                                              gTAGTCGCACATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGTGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGa                                                     >  2:565754/1‑150 (MQ=255)                                                                                     cccATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGa                                              >  2:1022113/1‑150 (MQ=255)                                                                                      ccATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGaa                                             <  2:229774/150‑1 (MQ=255)                                                                                       cATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTAcggcgg                                                                                     >  2:255387/1‑109 (MQ=255)                                                                                       cATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTAcggcgg                                                                                     <  1:255387/109‑1 (MQ=255)                                                                                          gcgcGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAgcgc                                         <  2:927454/150‑1 (MQ=255)                                                                                             cgGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTcc                                      <  1:1186148/150‑1 (MQ=255)                                                                                                  tgctgGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACTGCGGGTAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGAc                                 <  2:295816/150‑1 (MQ=255)                                                                                                                     ggcggcCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAgg                <  1:77311/148‑1 (MQ=255)                                                                                                                     ggcggcCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAgg                >  2:77311/1‑148 (MQ=255)                                                                                                                         gcCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATgg          <  1:910477/150‑1 (MQ=255)                                                                                                                          gcggcgcCTGCAGCCTGCCGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGc         >  1:272788/4‑150 (MQ=255)                                                                                                                           cAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCt        <  1:464928/150‑1 (MQ=255)                                                                                                                            ggcgcCTGCAGCCTGCCGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCTg       <  2:272788/149‑1 (MQ=255)                                                                                                                             gcgcCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCTgt      >  2:780552/1‑150 (MQ=255)                                                                                                                                 cGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCTgttctc  <  1:1110838/150‑1 (MQ=255)                                                                                                                                  |                                                                                                                                                     AAGCCGACCTGCCAGCCGGTGACCTGGCTGAAGGCGTTGATCGCCGACACCAGCGGCTGCTGGCCGATGCTGAAGCGGTAGTCGCCCATGCGCGGGCTGCTGGCCTGGGCCGGCTCGGCGGCCAGCGCCGGCAGGCTGCAGGCGCCGCTGGCGAGCAGGGCAAGGGTCAACAGGGACAGTTGTCCTGTACGGCGGGGAAGGTTGGACGAGCGGGTAGGACCAGTGGACATCGGAAGCGCTCCCAGACGCGCGGGCTGTTATAGGTGATGGCTGTTCTC  >  NC_002947/1147839‑1148116

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