Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,988,081	T→C	D329G (GAC→GGC)	hemX ←	PF04375 family protein HemX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,988,081	0	T	C	100.0%	163.6 / NA	45	D329G (GAC→GGC)	hemX	PF04375 family protein HemX
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/25);  total (20/25)

CGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGA  >  NC_000913/3987948‑3988212                                                                                                                                      |                                                                                                                                    cGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAg                                                                                                                                 <  1:67619/138‑1 (MQ=255) cGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAg                                                                                                                                 >  2:67619/1‑138 (MQ=255) cGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGa                                                                                                                               >  2:774913/1‑140 (MQ=255)    tGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACg                                                                                                                            <  2:277767/140‑1 (MQ=255)       ccccGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTAGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCtt                                                                                                                         <  1:217656/140‑1 (MQ=255)              gCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTtggtggt                                                                                                                  <  2:958598/140‑1 (MQ=255)               ccGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTtggtggtg                                                                                                                 <  1:763813/140‑1 (MQ=255)                 ggTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGc                                                                                                               <  2:21492/140‑1 (MQ=255)                 ggTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGc                                                                                                               <  2:925556/140‑1 (MQ=255)                    tgtgCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGcatc                                                                                                            <  2:903576/140‑1 (MQ=255)                        ccAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCAtcatca                                                                                                        <  2:230932/140‑1 (MQ=255)                            cagGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTAt                                                                                                    >  2:203616/1‑140 (MQ=255)                            cagGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTAt                                                                                                    <  1:776527/140‑1 (MQ=255)                            cagGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTAt                                                                                                    <  1:774913/140‑1 (MQ=255)                                ttACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTAt                                                                                                    >  1:944118/1‑136 (MQ=255)                                ttACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTAt                                                                                                    <  2:944118/136‑1 (MQ=255)                                  acgcacgcGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCgtagt                                                                                              <  1:696431/140‑1 (MQ=255)                                  acgcacgcGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTCTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCgtagt                                                                                              <  2:149727/140‑1 (MQ=255)                                   cgcacgcGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCgtagta                                                                                             >  1:235050/1‑140 (MQ=255)                                         cgAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGc                                                                                          <  1:835755/137‑1 (MQ=255)                                         cgAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGc                                                                                          >  2:835755/1‑137 (MQ=255)                                         cgAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACg                                                                                       <  2:368275/140‑1 (MQ=255)                                         cgAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACg                                                                                       >  2:968942/1‑140 (MQ=255)                                            gTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTAc                                                                                    <  1:306958/140‑1 (MQ=255)                                             tCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTAcc                                                                                   <  1:788044/140‑1 (MQ=255)                                                      gTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGgag                                                                          <  1:446944/140‑1 (MQ=255)                                                      gTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAACACGTACCCAGGTGgag                                                                          >  1:871055/1‑140 (MQ=255)                                                                cATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCa                                                                <  2:1007958/140‑1 (MQ=255)                                                                cATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCa                                                                >  1:871418/1‑140 (MQ=255)                                                                     ccTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTgg                                                                            >  2:606512/1‑123 (MQ=255)                                                                     ccTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTgg                                                                            <  1:606512/123‑1 (MQ=255)                                                                     ccTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGcc                                                           >  2:1053156/1‑140 (MQ=255)                                                                              gCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGAt                                                    <  2:755966/138‑1 (MQ=255)                                                                              gCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGAt                                                    >  1:755966/1‑138 (MQ=255)                                                                                             gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGa                                   >  2:245577/1‑140 (MQ=255)                                                                                              aTCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAc                                  >  2:55029/1‑140 (MQ=255)                                                                                              aTCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGGGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAc                                  >  2:373359/1‑140 (MQ=255)                                                                                               tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAcg                                 <  1:207320/140‑1 (MQ=255)                                                                                               tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAcg                                 <  1:975574/140‑1 (MQ=255)                                                                                                        atatTTTGCTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGc                        >  1:125021/1‑140 (MQ=255)                                                                                                                cTGGCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgctg                <  2:871418/140‑1 (MQ=255)                                                                                                                   gCTTAACTGGTCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGa             >  2:870617/1‑140 (MQ=255)                                                                                                                             tCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCg   >  2:751704/1‑140 (MQ=255)                                                                                                                             tCCACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCg   >  2:697123/1‑140 (MQ=255)                                                                                                                              ccACCTCGCCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCga  >  2:732256/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                    CGTTGTCCCCGCTGCCGGTTGTGCCAGCAGGTTACGCACGCGAGTCTGCATCAGTTTTTCCAGCATCGCCTGGCTTTGCAGGGTTTCCGGAAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGA  >  NC_000913/3987948‑3988212

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