Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,250,097	T→C	R7R (CGT→CGC)	yqjK →	conserved protein YqjK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,250,097	0	T	C	100.0%	126.9 / NA	38	R7R (CGT→CGC)	yqjK	conserved protein YqjK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/16);  total (22/16)

GCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGG  >  NC_000913/3249968‑3250232                                                                                                                                  |                                                                                                                                        gCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTa                                                                                                                               >  2:720134/1‑140 (MQ=255)       aTTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGgcgc                                                                                                                         >  1:712059/1‑140 (MQ=255)        ttGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCa                                                                                                                        <  1:315229/140‑1 (MQ=255)           gCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACt                                                                                                                     >  2:418178/1‑140 (MQ=255)                 tCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAg                                                                                                               >  1:363755/1‑140 (MQ=255)                     gACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAg                                                                                                           >  1:40198/1‑140 (MQ=255)                          tACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCa                                                                                                      <  1:259142/140‑1 (MQ=255)                           aCGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAg                                                                                                     >  2:710442/1‑140 (MQ=255)                           aCGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAg                                                                                                     >  1:122057/1‑140 (MQ=255)                               aaaTCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGcaac                                                                                                 <  1:720134/140‑1 (MQ=255)                                aaTCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGcaaca                                                                                                <  1:418178/140‑1 (MQ=255)                                   cgcgTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACg                                                                                             >  1:548053/1‑140 (MQ=255)                                      gTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAACAACAACGGCt                                                                                          >  2:507697/1‑140 (MQ=255)                                                    ctgcGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCTGCAGCAGAAAGGCCGAACGCGACCGACGCAAGGCGCAACTCCTTACCCAGATCCAGCAACAAGGGCTGAACCTTTCCCCCa                                                                            >  2:362843/1‑140 (MQ=255)                                                    ctgcGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCa                                                                            >  1:579164/1‑140 (MQ=255)                                                     tgcGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAg                                                                           <  1:379813/140‑1 (MQ=255)                                                               acGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGaa                                                                 <  1:513154/140‑1 (MQ=255)                                                                 gCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAg                                                                           >  2:33369/1‑128 (MQ=255)                                                                 gCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAg                                                                           <  1:33369/128‑1 (MQ=255)                                                                   cATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGc                                                             >  2:618205/1‑140 (MQ=255)                                                                      gAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTgg                                                          >  1:634475/1‑140 (MQ=255)                                                                           aGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGaca                                                     >  2:630072/1‑140 (MQ=255)                                                                              aaaCGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGacaaca                                                  <  2:363755/140‑1 (MQ=255)                                                                                aCGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGACCGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAgg                                                <  2:634475/140‑1 (MQ=255)                                                                                aCGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAgg                                                <  1:618205/140‑1 (MQ=255)                                                                                 cGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGgc                                               >  1:720839/1‑140 (MQ=255)                                                                                   aTCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGgcgc                                             <  2:122057/140‑1 (MQ=255)                                                                                                    gaggagTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGa                            <  2:579164/140‑1 (MQ=255)                                                                                                      ggagTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAat                          <  2:40198/140‑1 (MQ=255)                                                                                                       gagTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAata                         >  2:684109/1‑140 (MQ=255)                                                                                                        agTCCCGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAatat                        <  2:720839/140‑1 (MQ=255)                                                                                                             cGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATAtgctgc                   >  2:368050/1‑140 (MQ=255)                                                                                                             cGTGAGCAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATAtgctgc                   >  2:329757/1‑140 (MQ=255)                                                                                                                   cAGTAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTc             >  2:12226/1‑140 (MQ=255)                                                                                                                      tAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTgc          >  2:370383/1‑140 (MQ=255)                                                                                                                      tAAAGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTgc          <  1:329757/140‑1 (MQ=255)                                                                                                                       aaaGTCGAACGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTgcg         <  1:710442/140‑1 (MQ=255)                                                                                                                              aaCGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTgg  >  1:502728/1‑140 (MQ=255)                                                                                                                              aaCGCGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTgg  >  1:324658/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                                        GCACTGATTGGCGGTATCTGGACGCTACGTAAATCGCGTAAGTCTACGTTGCTGCGCCATACACGCCATGAGTTAGCAAACGATCGGCAGCTGCTCGAGGAGGAGTCCCGTGAGCAGTAAAGTCGAACGTGAACGACGTAAGGCGCAACTGCTTAGCCAGATCCAGCAACAACGGCTGGATCTTTCCGCCAGTCGTCGTGAATGGCTGGAGACAACAGGCGCTTACGATCGTCGCTGGAATATGCTGCTAAGTCTGCGCTCCTGG  >  NC_000913/3249968‑3250232

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