Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,899,468	G→A	N370N (AAC→AAT)	yieL ←	putative hydrolase YieL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,899,468	0	G	A	100.0%	143.2 / NA	41	N370N (AAC→AAT)	yieL	putative hydrolase YieL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/21);  total (20/21)

AGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGGTTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGCTGG  >  NC_000913/3899332‑3899600                                                                                                                                         |                                                                                                                                     aGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTc                                                                                                                                   >  1:421785/1‑140 (MQ=255)    aGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGa                                                                                                                                <  1:715644/140‑1 (MQ=255)      aTGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGAcc                                                                                                                              >  1:120818/1‑140 (MQ=255)      aTGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGAcc                                                                                                                              >  1:42728/1‑140 (MQ=255)          tgtgCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGa                                                                                                                          <  2:42728/140‑1 (MQ=255)             gCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGatat                                                                                                                       <  2:120818/140‑1 (MQ=255)                gcgACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTcc                                                                                                                    >  2:405813/1‑140 (MQ=255)                 cgACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCt                                                                                                                   >  1:255491/1‑140 (MQ=255)                   acacGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTgg                                                                                                                 <  2:369123/140‑1 (MQ=255)                         cTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGat                                                                                                                         >  1:296641/1‑126 (MQ=255)                         cTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGat                                                                                                                         <  2:296641/126‑1 (MQ=255)                                   ggTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAAttt                                                                                                 >  2:100168/1‑140 (MQ=255)                                             ggTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTc                                                                                       >  1:644157/1‑140 (MQ=255)                                                  tGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCt                                                                                  <  1:34846/140‑1 (MQ=255)                                                  tGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCt                                                                                  >  2:473616/1‑140 (MQ=255)                                                    aTTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCa                                                                                >  1:178078/1‑140 (MQ=255)                                                              tCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGcc                                                                      <  1:462497/140‑1 (MQ=255)                                                                    tGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCGGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGa                                                                >  2:823992/1‑140 (MQ=255)                                                                          attttatttAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCt                                                          >  1:509111/1‑140 (MQ=255)                                                                           ttttatttAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCtt                                                         <  1:100168/140‑1 (MQ=255)                                                                            tttatttaAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTg                                                        <  1:183906/140‑1 (MQ=255)                                                                              tatttaAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTg                                                        <  1:387630/138‑1 (MQ=255)                                                                              tatttaAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTg                                                        >  2:387630/1‑138 (MQ=255)                                                                                 ttaAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGt                                                   <  2:178078/140‑1 (MQ=255)                                                                                     ataatTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACg                                               >  2:481988/1‑140 (MQ=255)                                                                                      taatTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGa                                              <  1:473616/140‑1 (MQ=255)                                                                                        atTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACa                                            <  1:656227/140‑1 (MQ=255)                                                                                               tACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTAt                                     <  1:600539/140‑1 (MQ=255)                                                                                               tACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATACCCTTGCCGGTTACGACATCTTTAt                                     <  1:481988/140‑1 (MQ=255)                                                                                                   aaGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATcccc                                 <  2:255491/140‑1 (MQ=255)                                                                                                         gCATAGGCTGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAAc                           >  2:252838/1‑140 (MQ=255)                                                                                                                 tGGTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTaa                   >  1:714920/1‑140 (MQ=255)                                                                                                                  ggTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTGATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTaaa                  >  2:124973/1‑140 (MQ=255)                                                                                                                  ggTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATccccc                                <  1:323677/126‑1 (MQ=255)                                                                                                                  ggTCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATccccc                                >  2:323677/1‑126 (MQ=255)                                                                                                                    tCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAAt                >  1:518668/1‑140 (MQ=255)                                                                                                                    tCGCCAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTGTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAAt                >  2:628766/1‑140 (MQ=255)                                                                                                                        cAGACATCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACg            <  1:291235/140‑1 (MQ=255)                                                                                                                              tCCATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTg      <  2:644157/140‑1 (MQ=255)                                                                                                                               ccATTTCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGc     <  1:252838/140‑1 (MQ=255)                                                                                                                                  tttCATGATTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGCTgg  <  2:714920/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                     AGCAGATGGTGTGCCGCGACACGGCTCATTTCCTGGTAATCTTCGGTAATGATTAATTTCATCTTAATGCCTTATTTTATTTAAATAATTTCTGTACAAAGGCTGCATAGGCTGGTCGCCAGACATCCATTTCATGGTTCAGACCCGGATATTCCTGGTAATCAAAGTTAATTTTTTTCTGCTCAAGCTCAGTTTTCAGCCCGGCGATATCCTTGCCGGTTACGACATCTTTATCCCCCACAACCACAGTAAAATTACGTAGTTGCTGG  >  NC_000913/3899332‑3899600

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: ‑