Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,009,043	G→A	D94N (GAC→AAC)	ygeY →	putative peptidase YgeY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,009,043	0	G	A	100.0%	109.3 / NA	33	D94N (GAC→AAC)	ygeY	putative peptidase YgeY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/20);  total (13/20)

AGTCGGCTTCGATAAAGTTGAAATCGACCCGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGGACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCGGTAAAATCATTAAAGACCTCGGTCTGGAAGA  >  NC_000913/3008916‑3009170                                                                                                                                |                                                                                                                                aGTCGGCTTCGATAAAGTTGAAATCGACCCGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTc                                                                                                                            >  1:591543/1‑133 (MQ=255)   tCGGCTTCGATAAAGTTGAAATCGACCCGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGa                                                                                                                          >  2:453222/1‑133 (MQ=255)    cGGCTTCGATAAAGTTGAAATCGACCCGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGAt                                                                                                                         <  2:358833/133‑1 (MQ=255)    cGGCTTCGATAAAGTTGAAATCGACCCGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTAGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGAt                                                                                                                         <  1:456019/133‑1 (MQ=255)                     aaaTCGACCCGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGa                                                                                                        <  1:92567/133‑1 (MQ=255)                     aaaTCGACCCGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGa                                                                                                        <  1:358026/133‑1 (MQ=255)                       aTCGACCCGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGaaa                                                                                                      <  1:327341/133‑1 (MQ=255)                           aCCCGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACtga                                                                                                  >  1:493546/1‑133 (MQ=255)                            cccGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACtgat                                                                                                 >  2:338524/1‑133 (MQ=255)                            cccGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACtgat                                                                                                 <  1:517892/133‑1 (MQ=255)                                aTGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAAc                                                                                             <  1:381120/133‑1 (MQ=255)                                aTGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAAc                                                                                             <  2:493546/133‑1 (MQ=255)                                 tGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACt                                                                                            <  2:117515/133‑1 (MQ=255)                                        cGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATtggt                                                                                     <  1:338524/133‑1 (MQ=255)                                               ggTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTcgcg                                                                              <  1:453222/133‑1 (MQ=255)                                                         gCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGa                                                                    >  1:458038/1‑133 (MQ=255)                                                                 ccGCGTCTGGCGGCAATGGACGCTTATATCGATACCGTCGGAATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAg                                                            <  2:179735/133‑1 (MQ=255)                                                                        tggtggCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCAt                                                     <  2:92479/133‑1 (MQ=255)                                                                             gCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCct                                                >  1:513591/1‑133 (MQ=255)                                                                               aaTGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCctct                                              <  2:513591/133‑1 (MQ=255)                                                                                   gACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATgg                                          >  2:539012/1‑133 (MQ=255)                                                                                       cTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAAGTGGCAAGGCCTCTATGGTTTa                                      >  1:305292/1‑133 (MQ=255)                                                                                            atCGATACCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCg                                 >  2:335075/1‑133 (MQ=255)                                                                                                  aCCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCGGTaaaa                           >  2:371485/1‑133 (MQ=255)                                                                                                  aCCGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCGGTaaaa                           >  2:98089/1‑133 (MQ=255)                                                                                                    cGTCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCGGTAAAATc                         >  2:47033/1‑133 (MQ=255)                                                                                                      tCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTTTGGCTTATGCCGGCAAAATCAt                       >  1:257103/1‑133 (MQ=255)                                                                                                      tCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCg                                 <  1:566311/123‑1 (MQ=255)                                                                                                      tCGGCATTGGCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCg                                 >  2:566311/1‑123 (MQ=255)                                                                                                               gCAACATCAAAAACTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCGGTAAAATCATTAAAGACCt              <  2:591543/133‑1 (MQ=255)                                                                                                                       aaaaaCTGGACCTTCCATCCGTACGAAGCCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCGGTAAAATCATTAAAGACCTCGGTCTgg      <  2:257103/133‑1 (MQ=255)                                                                                                                          aaCTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCGGTAAAATCATTAAAGACCTCGGTCTGGAAg   <  1:539012/133‑1 (MQ=255)                                                                                                                           aCTGGAACTTCGATCCGTACGAATTCATCGAATCTGATGAAGTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCGGTAAAATCATTAAAGACCTCGGTCTGGAAGa  <  2:305292/133‑1 (MQ=255)                                                                                                                           aCTGGAACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCGGTAAAATCATTAAAGACCTCGGTCTGGAAGa  <  1:335075/133‑1 (MQ=255)                                                                                                                                |                                                                                                                                AGTCGGCTTCGATAAAGTTGAAATCGACCCGATGGGCAACGTTCTCGGTTATATCGGCCACGGCCCGCGTCTGGTGGCAATGGACGCTCATATCGATACCGTCGGCATTGGCAACATCAAAAACTGGGACTTCGATCCGTACGAAGGCATGGAAACTGATGAACTGATTGGTGGTCGCGGTACTTCCGACCAGGAAGGTGGCATGGCCTCTATGGTTTATGCCGGTAAAATCATTAAAGACCTCGGTCTGGAAGA  >  NC_000913/3008916‑3009170

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