Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,499,007	C→T	L394L (CTG→TTG)	cysG →	siroheme synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,499,007	0	C	T	100.0%	116.8 / NA	34	L394L (CTG→TTG)	cysG	siroheme synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/17);  total (17/17)

AAAACCTGGCGGCAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGAT  >  NC_000913/3498933‑3499083                                                                           |                                                                             aaaaCCTGGCGGCAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGAtt                                                                         >  2:619435/1‑80 (MQ=255) aaaaCCTGGCGGCAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTg                                                                        <  1:619388/81‑1 (MQ=255)     ccTGGCGGCAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAaca                                                                    >  1:619416/1‑81 (MQ=255)        ggcggcAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACg                                                                  >  1:619420/1‑80 (MQ=255)           ggcAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAAt                                                              <  1:619389/81‑1 (MQ=255)             cAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGc                                                            >  2:619430/1‑81 (MQ=255)                 aaaaCAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAgg                                                        <  2:619390/81‑1 (MQ=255)                      aGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGaaa                                                   <  2:619391/81‑1 (MQ=255)                       gACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAAt                                                  <  2:619392/81‑1 (MQ=255)                       gACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAAt                                                  >  1:619443/1‑81 (MQ=255)                          gCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGcc                                               >  2:619441/1‑81 (MQ=255)                          gCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGcc                                               <  2:619393/81‑1 (MQ=255)                          gCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGcc                                               >  1:619428/1‑81 (MQ=255)                           cTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCg                                              <  1:619394/81‑1 (MQ=255)                              gtgtTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTg                                           >  1:619429/1‑81 (MQ=255)                              gtgtTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTg                                           <  2:619395/81‑1 (MQ=255)                               tgtTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTgg                                          <  1:619396/81‑1 (MQ=255)                                gtTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTgg                                          >  1:619453/1‑80 (MQ=255)                                 ttCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCa                                        <  2:619397/81‑1 (MQ=255)                                 ttCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCa                                        <  1:619398/81‑1 (MQ=255)                                     atatGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTg                                    >  1:619442/1‑81 (MQ=255)                                       atGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTc                                  >  1:619431/1‑81 (MQ=255)                                           ggTTGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGaaa                              <  2:619399/81‑1 (MQ=255)                                              tGAATCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACg                           >  1:619444/1‑81 (MQ=255)                                                  tCAGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTAc                       <  2:619400/81‑1 (MQ=255)                                                    aGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACgg                     >  2:619432/1‑81 (MQ=255)                                                    aGGCCGCGACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACgg                     <  1:619401/81‑1 (MQ=255)                                                         gcgACTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGt                 >  2:619448/1‑80 (MQ=255)                                                            aCTATTCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACg             <  2:619402/81‑1 (MQ=255)                                                                ttCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAgc         <  1:619403/81‑1 (MQ=255)                                                                 tCAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAgcg        >  1:619437/1‑81 (MQ=255)                                                                  cAGCAAAAGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAgcgc       >  2:619447/1‑81 (MQ=255)                                                                       aaaGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGAt  >  1:619445/1‑81 (MQ=255)                                                                       aaaGTTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGAt  <  1:619404/81‑1 (MQ=255)                                                                           |                                                                             AAAACCTGGCGGCAGAAAAACAGACGCTGGTGTTCTATATGGGGTTGAATCAGGCCGCGACTATTCAGCAAAAGCTGATTGAACACGGAATGCCAGGCGAAATGCCGGTGGCAATTGTCGAAAACGGTACGGCAGTCACGCAGCGCGTGAT  >  NC_000913/3498933‑3499083

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