Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,222,371	A→G	A88A (GCT→GCC)	ymgD ←	PF16456 family protein YmgD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,371	0	A	G	100.0%	122.0 / NA	34	A88A (GCT→GCC)	ymgD	PF16456 family protein YmgD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/19);  total (15/19)

ATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTG  >  NC_000913/1222236‑1222497                                                                                                                                        |                                                                                                                               atatAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGgcgc                                                                                                                            <  1:795097/140‑1 (MQ=255)    taAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCaaa                                                                                                                         <  2:141760/140‑1 (MQ=255)             aTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCa                                                                                                                >  2:13315/1‑140 (MQ=255)             aTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCa                                                                                                                >  1:693545/1‑140 (MQ=255)               cccGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGt                                                                                                              <  2:371572/140‑1 (MQ=255)                 cgccgTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGttt                                                                                                            <  1:590220/140‑1 (MQ=255)                                 aCCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGgcg                                                                                            <  2:258009/140‑1 (MQ=255)                                                    cgcgGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCt                                                                         >  1:589099/1‑140 (MQ=255)                                                    cgcgGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCt                                                                         <  1:548062/140‑1 (MQ=255)                                                                  aaTATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCa                                                           <  2:22845/140‑1 (MQ=255)                                                                     attattTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACgg                                                        >  1:618009/1‑140 (MQ=255)                                                                     attattTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACgg                                                        >  2:430448/1‑140 (MQ=255)                                                                     attattTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACgg                                                        >  1:827709/1‑140 (MQ=255)                                                                      ttattTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGc                                                       >  2:696608/1‑140 (MQ=255)                                                                          tttttGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGc                                                   <  2:290883/140‑1 (MQ=255)                                                                               gCCGATAGCACGCATGGTGTCAGCAATTGCCGTGATCGACAGCAGTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCCCTTTGGCGCTCAGTTCCGGAGGTTGCTCAATTTTTAAATCACGGCTGGCAACGc                                              <  2:357788/140‑1 (MQ=255)                                                                                  gATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTgg                                           <  2:865508/140‑1 (MQ=255)                                                                                  gATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTgg                                           <  1:610678/140‑1 (MQ=255)                                                                                    tAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGca                                         <  1:13315/140‑1 (MQ=255)                                                                                    tAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGca                                         <  2:589099/140‑1 (MQ=255)                                                                                    tAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGca                                         <  2:827709/140‑1 (MQ=255)                                                                                    tAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGca                                         <  1:430448/140‑1 (MQ=255)                                                                                    tAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGca                                         <  2:693545/140‑1 (MQ=255)                                                                                         cGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGcattacc                                    <  2:618009/140‑1 (MQ=255)                                                                                                tgtCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGcattacccattacc                             >  1:723260/1‑140 (MQ=255)                                                                                                        aTTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAAt                     <  1:696608/140‑1 (MQ=255)                                                                                                        aTTGCCGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAAt                     >  2:385649/1‑140 (MQ=255)                                                                                                            ccGTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAATTTCt                 >  2:524031/1‑140 (MQ=255)                                                                                                              gTGATCAACAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGc               <  2:320918/140‑1 (MQ=255)                                                                                                                      cAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGc               >  2:621559/1‑132 (MQ=255)                                                                                                                      cAGCATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGc               <  1:621559/132‑1 (MQ=255)                                                                                                                         cATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTgcgc    >  2:603784/1‑140 (MQ=255)                                                                                                                         cATTTGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTgcgc    >  2:603603/1‑140 (MQ=255)                                                                                                                           tttGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTg  >  2:659319/1‑140 (MQ=255)                                                                                                                           tttGCGGGTCTTTGGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTg  >  2:770364/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                               ATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTG  >  NC_000913/1222236‑1222497

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