Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,991,035	C→T	intergenic (‑269/‑118)	hemC ← / → cyaA	hydroxymethylbilane synthase/adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,991,035	0	C	T	100.0%	146.4 / NA	40	intergenic (‑269/‑118)	hemC/cyaA	hydroxymethylbilane synthase/adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/16);  total (24/16)

TTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGC  >  NC_000913/3990960‑3991111                                                                            |                                                                             tttctttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGgc                                                                         >  1:783954/1‑81 (MQ=255)    ctttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGgcgc                                                                       <  2:783908/80‑1 (MQ=255)     tttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGcgcg                                                                      >  1:783957/1‑80 (MQ=255)       tACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGaa                                                                    >  2:783952/1‑80 (MQ=255)       tACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGaa                                                                    <  2:783909/80‑1 (MQ=255)        aCGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGaa                                                                    >  1:783961/1‑79 (MQ=255)             cAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTgg                                                              >  2:783948/1‑80 (MQ=255)              aaTCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTa                                                            <  1:783911/81‑1 (MQ=255)              aaTCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTa                                                            <  1:783910/81‑1 (MQ=255)                 cAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACg                                                         >  1:783938/1‑81 (MQ=255)                 cAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACg                                                         >  1:783975/1‑81 (MQ=255)                    cAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGtt                                                      <  1:783912/81‑1 (MQ=255)                    cAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGt                                                       >  1:783993/1‑80 (MQ=255)                     aaGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTa                                                     >  2:783943/1‑81 (MQ=255)                           ttAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCttt                                                <  2:783913/80‑1 (MQ=255)                           ttAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTg                                               <  2:783914/81‑1 (MQ=255)                           ttAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTg                                               <  1:783915/81‑1 (MQ=255)                              aaTTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGAc                                             <  1:783916/80‑1 (MQ=255)                                 tGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATAc                                         >  1:783972/1‑81 (MQ=255)                                   aTCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGa                                       <  1:783917/81‑1 (MQ=255)                                    tCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGa                                       >  1:783966/1‑80 (MQ=255)                                     cACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGaaa                                     >  2:783971/1‑81 (MQ=255)                                        gTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAatat                                  >  1:783958/1‑81 (MQ=255)                                        gTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAatat                                  <  1:783918/81‑1 (MQ=255)                                         ttttAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATc                                 >  2:783963/1‑81 (MQ=255)                                          tttAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATcc                                >  1:783941/1‑81 (MQ=255)                                           ttAGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATccc                               >  2:783955/1‑81 (MQ=255)                                             aGACCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCg                              >  2:783953/1‑80 (MQ=255)                                               aCCATTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAAt                           <  2:783919/81‑1 (MQ=255)                                                  aTTTTTTCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGcc                        <  2:783920/81‑1 (MQ=255)                                                    tttttCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCcgc                      >  1:783969/1‑81 (MQ=255)                                                      tttCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCgt                    >  1:783959/1‑81 (MQ=255)                                                      tttCGTCGTGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCgt                    >  1:783956/1‑81 (MQ=255)                                                             gtGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCg             >  1:783968/1‑81 (MQ=255)                                                              tGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGt            <  2:783921/81‑1 (MQ=255)                                                              tGAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGt            >  2:783977/1‑81 (MQ=255)                                                               gAAACTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGtt           >  1:783973/1‑81 (MQ=255)                                                                  aCTAAAAAAATCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGAt        <  1:783922/81‑1 (MQ=255)                                                                     aaaaaaaTCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGtt     >  1:783980/1‑81 (MQ=255)                                                                        aaaaTCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGc  <  1:783923/81‑1 (MQ=255)                                                                            |                                                                             TTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGC  >  NC_000913/3990960‑3991111

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