Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,157,770	A→G	T41A (ACT→GCT)	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,770	0	A	G	100.0%	164.2 / NA	45	T41A (ACT→GCT)	yqhG	DUF3828 domain‑containing protein YqhG
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/23);  total (22/23)

GAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTT  >  NC_000913/3157636‑3157896                                                                                                                                       |                                                                                                                               gagGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGcc                                                                                                                           <  2:797111/140‑1 (MQ=255)     tGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCtt                                                                                                                       >  2:481462/1‑140 (MQ=255)     tGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCtt                                                                                                                       >  2:1190893/1‑140 (MQ=255)          aaCAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGCCCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTAtttt                                                                                                                  <  2:1085161/140‑1 (MQ=255)              gATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTg                                                                                                              <  2:934060/140‑1 (MQ=255)              gATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTg                                                                                                              <  2:590723/140‑1 (MQ=255)              gATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTg                                                                                                              >  1:590723/1‑140 (MQ=255)               aTGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCt                                                                                                                              >  1:858101/1‑123 (MQ=255)               aTGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGa                                                                                                             <  2:818845/140‑1 (MQ=255)                         taCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGagag                                                                                                   <  1:271945/140‑1 (MQ=255)                         taCTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGagag                                                                                                   >  2:804096/1‑140 (MQ=255)                            ttCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCgg                                                                                                <  1:964981/140‑1 (MQ=255)                                 tttttAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGAt                                                                                           >  2:928134/1‑140 (MQ=255)                                          gcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTgc                                                                                  <  1:1047109/140‑1 (MQ=255)                                          gcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTgc                                                                                  <  1:1190893/140‑1 (MQ=255)                                          gcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTgc                                                                                  <  1:1202011/140‑1 (MQ=255)                                          gcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTgc                                                                                  <  2:44061/140‑1 (MQ=255)                                          gcCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTgc                                                                                  <  1:804096/140‑1 (MQ=255)                                           cccTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTgcg                                                                                 <  2:58188/140‑1 (MQ=255)                                                gCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGcgcg                                                                               >  1:141988/1‑137 (MQ=255)                                                gCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGcgcg                                                                               <  2:141988/137‑1 (MQ=255)                                                     ttttACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCTGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAAc                                                                       >  2:528227/1‑140 (MQ=255)                                                     ttttACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGCCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAAc                                                                       >  2:237598/1‑140 (MQ=255)                                                                acGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTAcc                                                            <  1:481462/140‑1 (MQ=255)                                                                       gCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACg                                                     >  2:571381/1‑140 (MQ=255)                                                                        ccccccTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGa                                                    >  1:290745/1‑140 (MQ=255)                                                                        ccccccTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGa                                                    >  1:24662/1‑140 (MQ=255)                                                                                  aGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTAccc                                                           >  1:620291/1‑123 (MQ=255)                                                                                  aGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTAccc                                                           <  2:620291/123‑1 (MQ=255)                                                                                   gACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGc                                         >  1:263911/1‑140 (MQ=255)                                                                                    aCCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCt                                        >  1:656461/1‑140 (MQ=255)                                                                                      cGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACATGCGCTTACCCTGAACGACAATCTTACGCTGc                                      >  1:1221797/1‑140 (MQ=255)                                                                                                     ccGGCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggct                       >  2:674229/1‑140 (MQ=255)                                                                                                       ggCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACATGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggctgg                     <  2:1221797/140‑1 (MQ=255)                                                                                                       ggCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggctgg                     <  1:655341/140‑1 (MQ=255)                                                                                                       ggCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggctgg                     <  1:687785/140‑1 (MQ=255)                                                                                                        gCATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATtggctggc                    <  2:554496/140‑1 (MQ=255)                                                                                                         cATATTTATCAGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTggctggct                   <  2:656461/140‑1 (MQ=255)                                                                                                                   aGAACTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGAt         <  2:24662/140‑1 (MQ=255)                                                                                                                      aCTATAAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCg      >  1:852873/1‑140 (MQ=255)                                                                                                                        tataAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGt    >  1:148434/1‑140 (MQ=255)                                                                                                                          taAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGttt  >  2:531235/1‑140 (MQ=255)                                                                                                                          taAATCAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGttt  >  2:1071248/1‑140 (MQ=255)                                                                                                                               cAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGAtt              >  2:664180/1‑123 (MQ=255)                                                                                                                               cAGATGCCGCTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGAtt              <  1:664180/123‑1 (MQ=255)                                                                                                                                       |                                                                                                                               GAGGTGAAAAACAGATGAAAATAATACTTCTGTTTTTAGCAGCCCTGGCAAGTTTTACCGTACACGCACAGCCCCCCTCACAGACCGTAGAACAAACAGTCCGGCATATTTATCAGAACTATAAATCAGATGCCACTGCCCCTTATTTTGGTGAAACCGGAGAGCGGGCGATAACTTCTGCGCGTATTCAACAGGCGCTTACCCTGAACGACAATCTTACGCTGCCGGGCAATATTGGCTGGCTGGATTATGATCCGGTTT  >  NC_000913/3157636‑3157896

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