Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,157,798	A→G	E50G (GAG→GGG)	yqhG →	DUF3828 domain‑containing protein YqhG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,157,798	0	A	G	100.0%	164.0 / NA	48	E50G (GAG→GGG)	yqhG	DUF3828 domain‑containing protein YqhG
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/22);  total (26/22)

CTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGGCGATCTGGTGCTAGA  >  NC_000913/3157665‑3157930                                                                                                                                      |                                                                                                                                     cTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAgggcgggc                                                                                                                                >  1:242303/1‑140 (MQ=255) cTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAgggcgggc                                                                                                                                >  1:562220/1‑140 (MQ=255)         agcagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACtt                                                                                                                        >  1:506596/1‑140 (MQ=255)           cagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCt                                                                                                                      >  1:63491/1‑140 (MQ=255)              cccTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCgg                                                                                                                                  >  2:583052/1‑125 (MQ=255)              cccTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCgg                                                                                                                                  <  1:583052/125‑1 (MQ=255)              cccTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTgcg                                                                                                                   >  1:771490/1‑140 (MQ=255)              cccTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGCCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTCTTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTgcg                                                                                                                   >  2:447006/1‑140 (MQ=255)                  ggCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACGATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAGCCGGAGGGCGGGCGATAACTTCTGCGCGTa                                                                                                               <  2:316798/140‑1 (MQ=255)                        ttttACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAAc                                                                                                         >  1:81365/1‑140 (MQ=255)                              cGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGgcgc                                                                                                   <  1:253121/140‑1 (MQ=255)                                  cacGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTAc                                                                                               >  2:833728/1‑140 (MQ=255)                                      cacaGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTg                                                                                           <  2:731018/140‑1 (MQ=255)                                          gCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACg                                                                                       >  2:38397/1‑140 (MQ=255)                                          gCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACg                                                                                       >  1:743122/1‑140 (MQ=255)                                           ccccccTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGa                                                                                      >  2:40063/1‑140 (MQ=255)                                           ccccccTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGa                                                                                      >  2:823224/1‑140 (MQ=255)                                           ccccccTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGa                                                                                      <  1:833728/140‑1 (MQ=255)                                           ccccccTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGa                                                                                      <  2:506596/140‑1 (MQ=255)                                               ccTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAAt                                                                                  >  1:613986/1‑140 (MQ=255)                                                        ccGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTg                                                                         <  1:64711/140‑1 (MQ=255)                                                             gAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCggg                                                                    <  2:675085/140‑1 (MQ=255)                                                             gAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCggg                                                                    <  2:434333/140‑1 (MQ=255)                                                             gAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCggg                                                                    <  1:823224/140‑1 (MQ=255)                                                                       tCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggc                                                          >  1:393345/1‑140 (MQ=255)                                                                       tCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggc                                                          <  2:242303/140‑1 (MQ=255)                                                                          ggCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggctgg                                                       <  1:40063/140‑1 (MQ=255)                                                                          ggCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggctgg                                                       <  2:81365/140‑1 (MQ=255)                                                                          ggCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggctgg                                                       <  1:415450/140‑1 (MQ=255)                                                                              tatTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGctggctgg                                                   >  2:318285/1‑140 (MQ=255)                                                                               atTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGa                                                  <  1:111556/140‑1 (MQ=255)                                                                                 ttATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGAtt                                                >  2:68010/1‑140 (MQ=255)                                                                                      aGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGAt                                           >  2:647284/1‑140 (MQ=255)                                                                                       gAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATc                                          <  2:238038/140‑1 (MQ=255)                                                                                         aCTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCg                                        >  1:370854/1‑140 (MQ=255)                                                                                          cTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGa                                                  <  2:96608/129‑1 (MQ=255)                                                                                          cTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGa                                                  >  1:96608/1‑129 (MQ=255)                                                                                          cTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCgg                                       >  2:304565/1‑140 (MQ=255)                                                                                             taAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCg                                        <  1:764934/136‑1 (MQ=255)                                                                                             taAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCg                                        >  2:764934/1‑136 (MQ=255)                                                                                                     aTGCCACTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGt                            <  2:171535/140‑1 (MQ=255)                                                                                                           cTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGAt                      >  2:803649/1‑140 (MQ=255)                                                                                                           cTGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGAt                      >  2:861753/1‑140 (MQ=255)                                                                                                            tGCCCCTTATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGAtt                     >  1:325979/1‑140 (MQ=255)                                                                                                                   tATTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGGCGa              >  2:753928/1‑140 (MQ=255)                                                                                                                    aTTTTGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGGCGAt             <  2:325979/140‑1 (MQ=255)                                                                                                                     ttttGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATGCAGCAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGGCGATc            <  1:447006/140‑1 (MQ=255)                                                                                                                        tGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGGCGATCTgg         >  2:390868/1‑140 (MQ=255)                                                                                                                        tGGTGAAACCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTACCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGGCGATCTgg         >  2:374657/1‑140 (MQ=255)                                                                                                                               aCCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGGCGATCTGGTGCTAGa  <  1:753928/140‑1 (MQ=255)                                                                                                                               aCCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGGCGATCTGGTGCTAGa  >  1:313819/1‑140 (MQ=255)                                                                                                                               aCCGGAGGGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGGCGATCTGGTGCTAGa  <  1:304565/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     CTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGGCGATCTGGTGCTAGA  >  NC_000913/3157665‑3157930

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