Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,765,203	A→T	Q119L (CAG→CTG)	yibJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,765,203	0	A	T	100.0%	132.1 / NA	44	Q119L (CAG→CTG)	yibJ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (22/22);  total (22/22)

CGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATC  >  NC_000913/3765184‑3765292                    |                                                                                          cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                  >  2:101330/1‑93 (MQ=9) cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                  >  2:5011093/1‑93 (MQ=9) cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                  <  2:1922262/93‑1 (MQ=9) cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                  >  2:2978704/1‑93 (MQ=9) cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                  >  1:314678/1‑93 (MQ=9) cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                  <  2:2076712/93‑1 (MQ=9) cGCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCg                  >  2:2106843/1‑93 (MQ=9)  gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                 <  1:1000704/93‑1 (MQ=12)  gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                 >  2:1252463/1‑93 (MQ=12)  gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                 <  1:5321531/93‑1 (MQ=12)  gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                 >  1:5275385/1‑93 (MQ=12)   cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                <  1:5116736/93‑1 (MQ=16)   cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                >  2:4397269/1‑93 (MQ=255)   cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                <  2:1204193/93‑1 (MQ=16)   cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                <  1:968374/93‑1 (MQ=16)    tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc               <  2:279119/93‑1 (MQ=255)    tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc               >  2:381217/1‑93 (MQ=255)    tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc               >  2:4075427/1‑93 (MQ=255)    tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc               >  2:4869776/1‑93 (MQ=255)    tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc               <  1:1131550/93‑1 (MQ=255)     attatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCa              >  1:3108470/1‑93 (MQ=255)      ttatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCAt             <  1:5810369/93‑1 (MQ=255)       tatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATa            >  1:2596157/1‑93 (MQ=255)       tatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATa            <  2:421572/93‑1 (MQ=255)       tatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATa            >  1:3031740/1‑93 (MQ=255)       tatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATa            <  2:3368983/93‑1 (MQ=255)        atGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaa           <  1:5880168/93‑1 (MQ=255)        atGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaa           <  1:3649111/93‑1 (MQ=255)         tGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa          >  2:2627383/1‑93 (MQ=255)         tGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa          <  2:2746975/93‑1 (MQ=255)          gACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt         >  2:790520/1‑93 (MQ=255)          gACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt         >  1:2882229/1‑93 (MQ=255)          gACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt         >  1:4831727/1‑93 (MQ=255)           aCCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATg        <  2:6014806/93‑1 (MQ=255)            cccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGt       <  1:2547539/93‑1 (MQ=255)            cccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGt       >  2:6061632/1‑93 (MQ=255)             ccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa      <  2:1726226/93‑1 (MQ=255)             ccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa      >  1:2849262/1‑93 (MQ=255)              cGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAg     >  1:776707/1‑93 (MQ=255)               gctgctGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa    <  2:2481334/93‑1 (MQ=255)               gctgctGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa    <  2:1102295/93‑1 (MQ=255)               gctgctGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa    <  2:5343596/93‑1 (MQ=255)                 tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc  >  1:3343845/1‑93 (MQ=255)                 tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc  <  1:3740832/93‑1 (MQ=255)                 tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc  <  2:1395627/93‑1 (MQ=255)                 tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc  <  2:6081903/93‑1 (MQ=255)                 tGCTGGGGCGGTATATCACTAAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAACTGTAGATc  <  2:401097/93‑1 (MQ=17)                    |                                                                                          CGCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATC  >  NC_000913/3765184‑3765292

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