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%	130.8 / NA	38	Q119L (CAG→CTG)	yibJ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (21/17);  total (21/17)

GCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATC  >  NC_000913/3765185‑3765292                   |                                                                                          gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                 <  1:903718/93‑1 (MQ=12) gCTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCgg                 >  1:2705620/1‑93 (MQ=12)  cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                >  1:1406105/1‑93 (MQ=255)  cTATTATGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGt                >  1:1455581/1‑93 (MQ=255)   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc               >  1:1102085/1‑93 (MQ=255)   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc               >  2:588439/1‑93 (MQ=255)   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc               <  2:3361153/93‑1 (MQ=255)   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc               >  1:724210/1‑93 (MQ=255)   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc               >  1:385517/1‑93 (MQ=255)   tattatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTc               <  1:2100136/93‑1 (MQ=255)    attatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCa              >  1:2422607/1‑93 (MQ=255)    attatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCa              >  1:1046655/1‑93 (MQ=255)     ttatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGCTGGAATTTTTATCAGTATCCGTTGAATCCTGTCAt             >  1:2924818/1‑93 (MQ=11)     ttatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCAt             >  2:3903504/1‑93 (MQ=255)     ttatGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCAt             >  2:3175831/1‑93 (MQ=255)       atGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaa           >  2:2723656/1‑93 (MQ=255)       atGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaa           <  2:594338/93‑1 (MQ=255)       atGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaa           >  1:3839617/1‑93 (MQ=255)        tGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGTGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa          <  1:1138441/93‑1 (MQ=255)        tGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa          <  2:3876510/93‑1 (MQ=255)        tGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa          <  1:4667936/93‑1 (MQ=255)        tGACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATaaa          >  2:1936527/1‑93 (MQ=255)         gACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt         <  1:3746574/93‑1 (MQ=255)         gACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt         <  2:3094229/93‑1 (MQ=255)         gACCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAAt         <  1:1332711/93‑1 (MQ=255)          aCCCGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATg        >  1:3332718/1‑93 (MQ=255)           cccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGt       <  2:377367/93‑1 (MQ=255)           cccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGt       <  2:3977309/93‑1 (MQ=255)            ccGCTTCTGGTGCGGTATATCACTCTGGCTCTGATTGAGCTTAAGGGATGATGGAATTTTTATCAGTATCCTTTGAATCCGGTCATAAATGTa      <  2:3362157/93‑1 (MQ=14)            ccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa      <  2:3259728/93‑1 (MQ=255)            ccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa      <  2:709071/93‑1 (MQ=255)            ccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa      >  2:4629564/1‑93 (MQ=255)            ccGCTGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa      >  2:4104390/1‑93 (MQ=255)            ccGCTGCTGGGGCGCTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTa      <  2:129729/93‑1 (MQ=255)              gctgctGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa    >  1:1926207/1‑93 (MQ=255)              gctgctGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa    >  1:2211283/1‑93 (MQ=255)              gctgctGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGa    <  2:2716884/93‑1 (MQ=255)               ctgctgGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGAt   <  1:2216024/93‑1 (MQ=255)               ctgctgGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGAt   <  2:3852599/93‑1 (MQ=255)                tGCTGGGGCGGTATATCACTCATGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc  <  2:372187/93‑1 (MQ=255)                tGCTGGGGCGGTATATCACTCAGGATCCGATTGGATTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc  <  1:3764544/93‑1 (MQ=255)                tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc  >  1:3217372/1‑93 (MQ=255)                tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc  >  2:3856122/1‑93 (MQ=255)                tGCTGGGGCGGTATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc  >  1:402558/1‑93 (MQ=255)                tGCTGGGGCGGTATATCACTCAGGATCCGATTAGACTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATc  <  1:2664625/93‑1 (MQ=255)                   |                                                                                          GCTATTATGACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCGGTCATAAATGTAGATC  >  NC_000913/3765185‑3765292

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