Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,640,281	T→G	D57E (GAT→GAG)	uspA →	universal stress global response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,640,281	0	T	G	100.0%	159.5 / NA	42	D57E (GAT→GAG)	uspA	universal stress global response regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/20);  total (22/20)

CGGTCGACCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGCAATCAAGAAATACGATATGGAT  >  NC_000913/3640136‑3640425                                                                                                                                                  |                                                                                                                                                 cgGTCGACCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCa                                                                                                                                             <  1:148728/151‑1 (MQ=255)  ggTCGACCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGc                                                                                                                                              >  1:266542/1‑149 (MQ=255)       aCCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACg                                                                                                                                       >  1:516192/1‑151 (MQ=255)        ccTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGc                                                                                                                                      <  2:362116/151‑1 (MQ=255)           ctcCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCa                                                                                                                                     <  1:230642/149‑1 (MQ=255)           ctcCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCAtc                                                                                                                                   >  2:269089/1‑151 (MQ=255)              cccGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATctct                                                                                                                                <  2:279098/151‑1 (MQ=255)                                   agagAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTg                                                                                                           <  1:717259/151‑1 (MQ=255)                                               ctctATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTcc                                                                                               >  1:784413/1‑151 (MQ=255)                                                   aTGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTa                                                                                           <  1:5632/151‑1 (MQ=255)                                                     ggCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTaa                                                                                          <  2:903499/150‑1 (MQ=255)                                                     ggCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTaa                                                                                          <  1:713352/150‑1 (MQ=255)                                                     ggCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAAt                                                                                         >  1:101674/1‑151 (MQ=255)                                                          gCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTAGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAgg                                                                                    >  1:69021/1‑151 (MQ=255)                                                                      cGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCAc                                                                        >  2:488219/1‑151 (MQ=255)                                                                        aaaGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACt                                                                       >  2:843161/1‑150 (MQ=255)                                                                           gTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGaa                                                                    >  1:778322/1‑150 (MQ=255)                                                                             ttCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAAc                                                                  <  2:441740/150‑1 (MQ=255)                                                                                tctGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTg                                                              >  1:437189/1‑151 (MQ=255)                                                                                  tGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTg                                                              >  1:497378/1‑149 (MQ=255)                                                                                   gATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGag                                                            <  1:68828/150‑1 (MQ=255)                                                                                      ccACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggc                                                        >  1:759050/1‑151 (MQ=255)                                                                                      ccACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggc                                                        <  2:83159/151‑1 (MQ=255)                                                                                        aCGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcag                                                      <  1:829366/151‑1 (MQ=255)                                                                                        aCGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggc                                                        >  2:632234/1‑149 (MQ=255)                                                                                         cGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGagcggcagc                                                     <  1:161393/151‑1 (MQ=255)                                                                                             gATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAgcggcg                                                 <  1:128719/151‑1 (MQ=255)                                                                                              aTGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGa                                                <  1:406629/151‑1 (MQ=255)                                                                                                 tAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCt                                             >  1:120817/1‑151 (MQ=255)                                                                                                   aaCTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTgg                                           <  1:286567/151‑1 (MQ=255)                                                                                                           tGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGtt                                   >  2:346798/1‑151 (MQ=255)                                                                                                              ccTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTg                                <  1:450508/151‑1 (MQ=255)                                                                                                               cTATACACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTgg                               >  1:773205/1‑151 (MQ=255)                                                                                                                tataCACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGt                              >  2:398162/1‑151 (MQ=255)                                                                                                                 ataCACCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTc                             >  2:517105/1‑151 (MQ=255)                                                                                                                    cacCGGGCTTATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGAt                          <  1:294236/151‑1 (MQ=255)                                                                                                                             tATTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGCAATCa                   >  1:151867/1‑149 (MQ=255)                                                                                                                              aTTGATGTGAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGCAATCa                   >  1:84844/1‑148 (MQ=255)                                                                                                                                    gtgAATCTGGGTGAGATGCAGAAACGCATCTCTGAAGAGACAGATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCGGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGACGATGCAATCAAGAAAt             <  2:353273/148‑1 (MQ=255)                                                                                                                                           tGGGTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGCAATCAAGAAATACGATATgg    >  2:451022/1‑150 (MQ=255)                                                                                                                                            gggTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGCAATCAAGAAATACGatat      <  1:468966/147‑1 (MQ=255)                                                                                                                                            gggTGAGATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGCAATCAAGAAATACGATATGGAt  >  1:288135/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CGGTCGACCTCTCCCCGGAAAGCAAAGTTCTGGTAGAGAAAGCAGTCTCTATGGCTCGCCCCTACAATGCGAAAGTTTCTCTGATCCACGTAGATGTAAACTACTCTGACCTATACACCGGGCTTATTGATGTGAATCTGGGTGATATGCAGAAACGCATCTCTGAAGAGACACATCATGCACTGACCGAGCTTTCCACTAATGCAGGCTACCCAATCACTGAAACCCTGAGCGGCAGCGGCGACCTGGGCCAGGTTCTGGTCGATGCAATCAAGAAATACGATATGGAT  >  NC_000913/3640136‑3640425

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