Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,040,801	T→C	L107L (TTA→CTA)	msrQ →	periplasmic protein‑L‑methionine sulfoxide reductase heme binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,040,801	0	T	C	100.0%	125.5 / NA	32	L107L (TTA→CTA)	msrQ	periplasmic protein‑L‑methionine sulfoxide reductase heme binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/16);  total (16/16)

ATTTAACAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGTTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAGGTAT  >  NC_000913/2040754‑2040847                                                |                                               aTTTAACAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGcta                                              >  1:943521/1‑50 (MQ=255)     aaCAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAg                                          >  1:1636900/1‑50 (MQ=255)       cAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAgg                                         >  1:1599303/1‑49 (MQ=255)        aaGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGt                                        >  1:2052211/1‑49 (MQ=255)           ttACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAgg                                    <  1:2800476/50‑1 (MQ=255)             aCGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGa                                   <  1:866725/49‑1 (MQ=255)             aCGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAg                                  >  1:1967243/1‑50 (MQ=255)              cGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGt                                 <  1:474842/50‑1 (MQ=255)               gCATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGt                                 >  1:2233746/1‑49 (MQ=255)                cATTGCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTa                               <  1:1168122/50‑1 (MQ=255)                    gCTGGAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTa                           >  1:2734209/1‑50 (MQ=255)                       ggAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTAccc                        >  1:2472414/1‑50 (MQ=255)                        gAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTAccc                        <  1:1905179/49‑1 (MQ=255)                        gAGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCg                       <  1:830391/50‑1 (MQ=255)                         aGCTGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCgg                      <  1:1437428/50‑1 (MQ=255)                            tGGGCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCt                   >  1:2210134/1‑50 (MQ=255)                             gggCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCt                   >  1:1250238/1‑49 (MQ=255)                             gggCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCt                   >  1:825374/1‑49 (MQ=255)                              ggCGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTa                 <  1:2165888/50‑1 (MQ=255)                                cGTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTAtt               <  1:2678600/50‑1 (MQ=255)                                 gTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTAttt              <  1:660802/50‑1 (MQ=255)                                 gTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTAtt               >  1:1057925/1‑49 (MQ=255)                                 gTGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTAtt               >  1:2804573/1‑49 (MQ=255)                                  tGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTAttt              >  1:1342158/1‑49 (MQ=255)                                  tGAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTa             <  1:2362087/50‑1 (MQ=255)                                   gAATAATCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTaa            >  1:188029/1‑50 (MQ=255)                                    aataatCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAAc           >  1:622737/1‑50 (MQ=255)                                    aataatCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAAc           <  1:613342/50‑1 (MQ=255)                                    aataatCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAAc           <  1:408828/50‑1 (MQ=255)                                      taatCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGt         <  1:2565757/50‑1 (MQ=255)                                        atCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGtt        <  1:1371397/49‑1 (MQ=255)                                        atCTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTa       <  1:221143/50‑1 (MQ=255)                                          cTGGCGCTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAgg     >  1:2017708/1‑50 (MQ=255)                                             gcgcTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAg      <  1:1194127/46‑1 (MQ=255)                                             gcgcTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAGGtat  <  1:2341873/50‑1 (MQ=255)                                             gcgcTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAGGta   <  1:1738436/49‑1 (MQ=255)                                             gcgcTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAGGt    <  1:2249272/48‑1 (MQ=255)                                             gcgcTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAGGt    <  1:1361186/48‑1 (MQ=255)                                                |                                               ATTTAACAAGTTACGCATTGCTGGAGCTGGGCGTGAATAATCTGGCGTTATTAGGTAAGGAGTTAATTACCCGGCCTTATTTAACGTTAGGTAT  >  NC_000913/2040754‑2040847

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