Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	819,067	(T)7→8	coding (20/246 nt)	moaD →	molybdopterin synthase sulfur carrier subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	819,060	1	.	T	100.0%	58.1 / NA	17	L5F (CTT→TTT)	moaD	molybdopterin synthase sulfur carrier subunit
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (5/12);  total (5/12)

GTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTC‑TTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCC  >  NC_000913/818981‑819124                                                                                 |                                                                 gTGATTGGTCCGGTACGTT‑‑‑TGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTC‑ttt                                                               <  2:156478/80‑1 (MQ=255)   gATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTC‑ttt                                                               >  2:156529/1‑81 (MQ=255)      tGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTC‑tttttt                                                            >  2:156513/1‑81 (MQ=255)      tGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTC‑tttttt                                                            >  1:156526/1‑81 (MQ=255)       ggTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTC‑ttttttt                                                           >  1:156516/1‑81 (MQ=255)        gTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTC‑ttttttt                                                           <  1:156479/80‑1 (MQ=255)             gTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGccc                                                      <  1:156480/81‑1 (MQ=255)               aCGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAg                                                    >  2:156522/1‑81 (MQ=255)                     gCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTgcg                                               <  2:156481/80‑1 (MQ=255)                         gCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGa                                            <  1:156482/79‑1 (MQ=255)                          cGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGtt                                         <  1:156483/81‑1 (MQ=255)                             agagCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTtggt                                      >  2:156519/1‑81 (MQ=255)                                gcggcgGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTTGGTggg                                   <  1:156484/81‑1 (MQ=255)                                    cggcAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACa                               <  2:156485/81‑1 (MQ=255)                                          gTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCa                         <  2:156486/81‑1 (MQ=255)                                                 gACTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAg                  >  2:156532/1‑81 (MQ=255)                                                  aCTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAg                  <  1:156487/80‑1 (MQ=255)                                                  aCTTTAAGGTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGt                 <  1:156488/81‑1 (MQ=255)                                                         ggTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCg          <  1:156489/81‑1 (MQ=255)                                                          gTGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCg          >  1:156523/1‑80 (MQ=255)                                                           tGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGa        >  2:156518/1‑81 (MQ=255)                                                           tGGAAGCGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGa        <  1:156490/81‑1 (MQ=255)                                                                 cGGATGATTAAAGTTCTTTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTccc  <  1:156491/81‑1 (MQ=255)                                                                                 |                                                                 GTGATTGGTCCGGTACGTTTGCTGGCGAAGAGCGGCGGCAAGTCGGGTGACTTTAAGGTGGAAGCGGATGATTAAAGTTC‑TTTTTTTCGCCCAGGTGCGCGAGTTGGTGGGAACAGATGCAACCGAAGTGGCTGCGGATTTCCC  >  NC_000913/818981‑819124

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: ‑