Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	58,877	A→G	N135S (AAT→AGT)	yabP →	putative uncharacterized protein YabP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	58,877	0	A	G	100.0%	102.4 / NA	30	N135S (AAT→AGT)	yabP	putative uncharacterized protein YabP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/16);  total (14/16)

TGGAATTATTCAATAAACTTAAAACATGCGCAAAAAATCCGGATGAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAG  >  NC_000913/58748‑59006                                                                                                                                  |                                                                                                                                  tGGAATTATTCAATAAACTTAAAACATGCGCAAAAAATCCGGATGAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTaa                                                                                                                                <  1:531476/133‑1 (MQ=255)        attCAATAAACTTAAAACATGCGCAAAAAATCCGGATGAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCAAGTTAAGTGTCAGTTTGTTCATTGAAAAGAGTAAGATATGt                                                                                                                         >  2:165123/1‑133 (MQ=255)         ttCAATAAACTTAAAACATGCGCAAAAAATCCGGATGAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTa                                                                                                                        >  2:532169/1‑133 (MQ=255)                    tAAAACATGCGCAAATAATCCGGATGAGTTCGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTg                                                                                                             <  1:302169/133‑1 (MQ=255)                       aaCATGCGCAAAAAACCCCGATCAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAg                                                                                                          <  1:165123/133‑1 (MQ=255)                         cATGCGCAAAAAATCCGGATGAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCt                                                                                                        <  2:203874/133‑1 (MQ=255)                           tGCGCAAAAAATCCGGATGAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTGCATTGAAAAGAGGAAGATATGTACTAAGACGTTGCAGCTGc                                                                                                      >  2:246641/1‑133 (MQ=255)                           tGCGCAAAAAATCCAGATGAGTTGGTCGCAGAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGc                                                                                                      >  2:383734/1‑133 (MQ=255)                             cgcAAAAAATCCGGATGAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCAt                                                                                                    <  1:532169/133‑1 (MQ=255)                              gcAAAAAATCCGGATGAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATAAGTACTCAGACGTTGCAGCTGCATa                                                                                                   <  2:416064/133‑1 (MQ=255)                                  aaaaTCCGGATGAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAggg                                                                                               <  2:310320/133‑1 (MQ=255)                                  aaaaTCCGGATGAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAggg                                                                                               <  1:448115/133‑1 (MQ=255)                                        cGGATGAGTTGGTCGGAAAGGGTAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAgg                                                                                         <  1:217641/133‑1 (MQ=255)                                                   gTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCg                                                                              <  1:1176/133‑1 (MQ=255)                                                     cGCAAAGGGCAATATCCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGtt                                                                            <  1:246641/133‑1 (MQ=255)                                                          aGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCt                                                                       >  2:535153/1‑133 (MQ=255)                                                             gCAATATGTGCTGGGGCGCCCACACCTTCGCGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATAAGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAg                                                                    <  1:170234/133‑1 (MQ=255)                                                             gCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAg                                                                    <  2:446860/133‑1 (MQ=255)                                                             gCAATATGAGCTGGGCCGTCCAAACCTTCGCGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAg                                                                    <  1:257085/133‑1 (MQ=255)                                                                      gcTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGGTGATCTTAGAAAGGCTTa                                                           >  1:238705/1‑133 (MQ=255)                                                                          ggTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCtt                                                       >  1:500932/1‑133 (MQ=255)                                                                                         cGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTAAGAAAGGGTTAACTGCCCGGCGTTGACATg                                        >  1:183362/1‑133 (MQ=255)                                                                                         cGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCATTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATg                                        >  2:411871/1‑133 (MQ=255)                                                                                              atatCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGa                                   >  2:367902/1‑133 (MQ=255)                                                                                                        aaGTGTCACTTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAGAGGCTTACCTTCCCGGGGTTGACATGAGGGATTGTTACCtt                         >  2:164855/1‑133 (MQ=255)                                                                                                             tCACCTTTTTCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTaa                    >  1:524205/1‑133 (MQ=255)                                                                                                                    tttCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACa             <  1:367902/133‑1 (MQ=255)                                                                                                                     ttCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACaa            >  1:316355/1‑133 (MQ=255)                                                                                                                     ttCATTGAAAAGAGTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACaa            >  1:105974/1‑133 (MQ=255)                                                                                                                          tGAAAAGAGTAAGATGTGTACTCACACGTTGCCGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGaaa       <  2:238705/133‑1 (MQ=255)                                                                                                                               agagTAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAg  <  2:105974/133‑1 (MQ=255)                                                                                                                                  |                                                                                                                                  TGGAATTATTCAATAAACTTAAAACATGCGCAAAAAATCCGGATGAGTTGGTCGCAAAGGGCAATATGCGCTGGGTCGCCCAGACCTTCGGGGATATCGAGTTAAGTGTCACTTTTTTCATTGAAAAGAATAAGATATGTACTCAGACGTTGCAGCTGCATAAGGGCCAAGGTAACTTGGGCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAG  >  NC_000913/58748‑59006

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