Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	336,297	A→G	T125A (ACC→GCC)	yahE →	DUF2877 domain‑containing protein YahE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	336,297	0	A	G	100.0%	126.9 / NA	35	T125A (ACC→GCC)	yahE	DUF2877 domain‑containing protein YahE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/19);  total (16/19)

TCAACATCTTCATATAGAGATGAGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAACCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACAAAATATCGCCTCAGCGGTCGCCAGTATGAT  >  NC_000913/336161‑336427                                                                                                                                         |                                                                                                                                   tCAACATCTTCATATAGAGATGAGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGTAAATGAAGccc                                                                                                                                 <  1:628599/140‑1 (MQ=255)          tCATATAGAGATGAGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAAtt                                                                                                                        >  2:527287/1‑140 (MQ=255)          tCATATAGAGATGAGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAAtt                                                                                                                        >  2:405734/1‑140 (MQ=255)           cATATAGAGATGAGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGccc                                                                                                                                 <  2:806890/130‑1 (MQ=255)           cATATAGAGATGAGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGccc                                                                                                                                 >  1:806890/1‑130 (MQ=255)                  agaTGAGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCgcggc                                                                                                                >  1:757709/1‑140 (MQ=255)                   gaTGAGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCgcggcg                                                                                                               >  2:4936/1‑140 (MQ=255)                    aTGAGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGa                                                                                                              >  1:622732/1‑140 (MQ=255)                      gAGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATa                                                                                                            <  1:111802/140‑1 (MQ=255)                       aGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATaa                                                                                                           >  2:605017/1‑140 (MQ=255)                                  gCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATc                                                                                                <  1:527287/140‑1 (MQ=255)                                  gCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATc                                                                                                <  1:405734/140‑1 (MQ=255)                                            tCCCCAACCTTGCAAATGACCGCCGTGAATTTTCACCTTATCGCCTTGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTtaaa                                                                                      >  2:793966/1‑140 (MQ=255)                                                cAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTtaaataaa                                                                                  <  2:757709/140‑1 (MQ=255)                                                   ccTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGaa                                                                               <  2:628665/140‑1 (MQ=255)                                                     ttGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAAtt                                                                             <  2:615316/140‑1 (MQ=255)                                                                cGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTGCATATTAAac                                                                  >  1:631756/1‑140 (MQ=255)                                                                    gtgAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAacgacg                                                              <  1:908065/140‑1 (MQ=255)                                                                    gtgAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAacgacg                                                              <  1:247286/140‑1 (MQ=255)                                                                      gAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACggg                                                            <  1:793966/140‑1 (MQ=255)                                                                       aaTTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGc                                                           <  1:337918/140‑1 (MQ=255)                                                                         ttttCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAg                                                         >  1:224677/1‑140 (MQ=255)                                                                              aCCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTAtt                                                    <  1:4936/140‑1 (MQ=255)                                                                                        cTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTaa                                          <  2:732825/140‑1 (MQ=255)                                                                                              aCAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATa                                    <  2:622732/140‑1 (MQ=255)                                                                                                  tGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACa                                >  1:103535/1‑140 (MQ=255)                                                                                                     cACGACATTATTCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACaaaa                             >  1:792076/1‑140 (MQ=255)                                                                                                           attattCATCAGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACAAAATATCGc                       >  2:559174/1‑140 (MQ=255)                                                                                                                tcatcaGCACCTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACAAAATATCGCCTCAg                  <  1:150741/140‑1 (MQ=255)                                                                                                                         ccTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACAAAATATCGCCTCAGCGGTCGCCa         >  1:527805/1‑140 (MQ=255)                                                                                                                         ccTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACAAAATATCGCCTCAGCGGTCGCCa         >  2:569415/1‑140 (MQ=255)                                                                                                                         ccTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACAAAATATCGCCTCAGCGGTCGCCa         >  1:809307/1‑140 (MQ=255)                                                                                                                         ccTGGGGGAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACAAAATATCGCCTCAGCGGTCGCCa         <  2:631756/140‑1 (MQ=255)                                                                                                                            gggggAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACAAAATATCGCCTCAGCGGTCGCCAGTa      >  1:892993/1‑140 (MQ=255)                                                                                                                               ggAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACAAAATATCGCCTCAGCGGTCGCCAGTatga   <  2:224677/140‑1 (MQ=255)                                                                                                                                gAAAATGAAGCCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACAAAATATCGCCTCAGCGGTCGCCAGTatgat  <  2:5916/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                   TCAACATCTTCATATAGAGATGAGTCGTTGTCGGCGTTGGCTGTCCCCAACCTTGCAAATGACCGCTGTGAATTTTCACCTTATCGCCTGGCTACAGTGGCACGACATTATTCATCAGCACCTGGGGGAAAATGAAACCCTGTTTAATTATCGCGGCGATAATCCGTTTTATCAGGCGTTAAATAAAGAATTACATATTAAACGACGGGCAGTTATTCAGGCCGTAAACGATAAACAAAATATCGCCTCAGCGGTCGCCAGTATGAT  >  NC_000913/336161‑336427

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