Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,157,770	A→G	T41A (ACT→GCT)	yqhG →	DUF3828 family putative periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,157,770	0	A	G	100.0%	97.7 / NA	28	T41A (ACT→GCT)	yqhG	DUF3828 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/10);  total (18/10)

TGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGG  >  NC_000913/3157640‑3157915                                                                                                                                   |                                                                                                                                                  tGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGa                                                                                                                                >  2:713834/1‑150 (MQ=255) tGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGa                                                                                                                                >  2:387849/1‑150 (MQ=255)  gAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGaa                                                                                                                               >  1:463932/1‑150 (MQ=255)       aCAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCgg                                                                                                                          >  1:613342/1‑150 (MQ=255)             gAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCg                                                                                                                    <  2:608137/150‑1 (MQ=255)                    ataCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATaa                                                                                                            <  1:412725/151‑1 (MQ=255)                    ataCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATa                                                                                                             >  2:109794/1‑150 (MQ=255)                          cTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTg                                                                                                      >  2:456411/1‑151 (MQ=255)                                    cagcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCa                                                                                            <  1:109796/151‑1 (MQ=255)                                     agcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCaa                                                                                           <  2:258432/151‑1 (MQ=255)                                           ggCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAgg                                                                                       <  2:575887/149‑1 (MQ=255)                                                       cGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAAc                                                                         <  2:578038/151‑1 (MQ=255)                                                          acacGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGAc                                                                      <  2:129839/151‑1 (MQ=255)                                                                              aGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCa                                                   >  2:93955/1‑150 (MQ=255)                                                                                aCCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAat                                                 >  2:502790/1‑150 (MQ=255)                                                                                 ccGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAatat                                               >  1:544254/1‑151 (MQ=255)                                                                                  cGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATAtt                                              >  1:283030/1‑151 (MQ=255)                                                                                        acaaacaGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggctgg                                        <  2:430229/151‑1 (MQ=255)                                                                                             caGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGa                                   <  2:493836/151‑1 (MQ=255)                                                                                                  cGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATg                              >  1:477970/1‑151 (MQ=255)                                                                                                      atatTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATcc                          >  1:336494/1‑151 (MQ=255)                                                                                                         tttATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGt                       >  1:753526/1‑151 (MQ=255)                                                                                                              cAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTtgtg                  >  2:568776/1‑151 (MQ=255)                                                                                                                    tataAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTc            >  2:196599/1‑151 (MQ=255)                                                                                                                     ataAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTc            >  1:177812/1‑150 (MQ=255)                                                                                                                         aTCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTCCGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGAt       <  1:497685/151‑1 (MQ=255)                                                                                                                            aGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGAtttt    >  2:511602/1‑151 (MQ=255)                                                                                                                              aTGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTgg  >  2:261962/1‑151 (MQ=255)                                                                                                                                   |                                                                                                                                                  TGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTTGTGATTGTCAGGATTTTGG  >  NC_000913/3157640‑3157915

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