Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,397,789	T→C	A278A (GCA→GCG)	uspE ←	universal stress protein with a role cellular motility

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,397,789	0	T	C	100.0%	129.3 / NA	38	A278A (GCA→GCG)	uspE	universal stress protein with a role cellular motility
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/10);  total (28/10)

AAGAGTCGTTAATCGTCTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCC  >  NC_000913/1397664‑1397925                                                                                                                               |                                                                                                                                         aaGAGTCGTTAATCGTCTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGt                                                                                                                            <  2:1344102/140‑1 (MQ=255) aaGAGTCGTTAATCGTCTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGt                                                                                                                            <  1:1214394/140‑1 (MQ=255) aaGAGTCGTTAATCGTCTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGt                                                                                                                            <  2:286550/140‑1 (MQ=255)     gTCGTTAATCGTCTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGt                                                                                                                            >  1:619288/1‑136 (MQ=255)     gTCGTTAATCGTCTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGt                                                                                                                            <  2:619288/136‑1 (MQ=255)     gTCGTTAATCGTCTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGt                                                                                                                        >  2:415689/1‑140 (MQ=255)        gTTAATCGTCTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCt                                                                                                                     >  1:942196/1‑140 (MQ=255)         ttAATCGTCTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTa                                                                                                                    <  2:872677/140‑1 (MQ=255)           aaTCGTCTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTAcc                                                                                                                  <  1:1054434/140‑1 (MQ=255)                        gTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACca                                                                                                     <  2:815279/140‑1 (MQ=255)                        gTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACca                                                                                                     <  2:701219/140‑1 (MQ=255)                                gTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATcccc                                                                                             <  1:415689/140‑1 (MQ=255)                                 ttCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCg                                                                                            >  1:1304928/1‑140 (MQ=255)                                 ttCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCg                                                                                            >  1:220189/1‑140 (MQ=255)                                 ttCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCg                                                                                            >  1:270252/1‑140 (MQ=255)                                 ttCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCg                                                                                            >  2:460722/1‑140 (MQ=255)                                      cGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCg                                                                                            <  2:105940/135‑1 (MQ=255)                                      cGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCg                                                                                            >  1:105940/1‑135 (MQ=255)                                           gtCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTg                                                                                  >  1:291118/1‑140 (MQ=255)                                           gtCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTg                                                                                  >  1:1239007/1‑140 (MQ=255)                                                 tACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCCCCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATCCCGGGGCGTCCTCCCGGGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGc                                                                            <  2:158887/140‑1 (MQ=255)                                                      gTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAAt                                                                       >  1:849165/1‑140 (MQ=255)                                                      gTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAAt                                                                       >  1:1309796/1‑140 (MQ=255)                                                      gTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAAt                                                                       >  2:331214/1‑140 (MQ=255)                                                      gTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAAt                                                                       >  1:1281339/1‑140 (MQ=255)                                                         cggTTTAAATGCCCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTCCCGCCGTGTGGCCGAGGAAGGCCGCTGAAATACCGGTGCGTCCTCCCGTGCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCa                                                                     <  2:55356/139‑1 (MQ=255)                                                           gTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGa                                                                  >  1:219285/1‑140 (MQ=255)                                                                aaTGCCCAGCAGGTCGCAGGAAAGATGATAAATCCCCTGTCCCGCCGTGTTGCCGAGGAATGCCGCTAAAATCCCGGTGCGTCCTCCCGTGCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCa                                                              <  2:219285/140‑1 (MQ=255)                                                                aaTCACCAGAGGGGCGCGGCAAAGATGACAAATCCCCTTTTCCGCCGTGTGGCGGAGGAATCCCCCTAAAATCCGGGTGCTCCCCACCGTCCCAAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCa                                                              <  2:942196/140‑1 (MQ=255)                                                                       agcagGCCGCAGCGAAGATGATCAATCCCCTGTCCCGCGGTGTTGCCGAGGAATGCCGCTAAAATCCGGGTGCGTCCTACCGTCCCCAGACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCtcttc                                                       <  2:1281339/140‑1 (MQ=255)                                                                           ggTCGCGGCAAAGTTGACCAATCACCTGTCCCGCCGTGTTGCCGAGGAATCCCGCTAAAATACCGGTGCGTCCTCCCGTGCCCAGACCCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGc                                                   <  2:270252/140‑1 (MQ=255)                                                                           ggTCGCAGCGAAGATGATCAATCCCCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATCCCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGc                                                   <  2:453754/140‑1 (MQ=255)                                                                              cGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGa                                                >  2:199476/1‑140 (MQ=255)                                                                              cGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGa                                                >  2:80145/1‑140 (MQ=255)                                                                               gCAGCGAAGATGACCAATCCCTTGTCCCGCCGTGTTGCCGAGGAATGCCGCTAAAATCCCGGTGCGCCCTCCCGTCCCCAAACCCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGAc                                               <  2:291118/140‑1 (MQ=255)                                                                                           atcaatcaCCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGt                                   >  2:45313/1‑140 (MQ=255)                                                                                              aatccCCTGTCCCGCCGTGTTGCCGAGGAATGCCGCTGAAATCCCGGTGCGTCCTCCCGTGCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCg                                <  2:1304928/140‑1 (MQ=255)                                                                                                cccccTGTCCCGCCGTGTTGCCGAGGAAGGCCGCTGAAATCCCGGTGCGCCCTCCCGTCCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTc                              <  2:220189/137‑1 (MQ=255)                                                                                                  aCCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCAt                            >  1:876899/1‑140 (MQ=255)                                                                                                  aCCTGTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCAt                            >  1:1272354/1‑140 (MQ=255)                                                                                                      gTTCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGttt                        >  2:22968/1‑140 (MQ=255)                                                                                                       ttCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGtttt                       >  1:552865/1‑140 (MQ=255)                                                                                                       ttCCGCCGTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGtttt                       >  1:525835/1‑140 (MQ=255)                                                                                                          cgccgTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCAt                    >  1:664252/1‑140 (MQ=255)                                                                                                            ccgTGTTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTa                  >  2:220081/1‑140 (MQ=255)                                                                                                              gtgtTGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAAt                >  2:1017699/1‑140 (MQ=255)                                                                                                                tttGCCGAGGAATGCCGCTGAAATCCCGGTGCGTCCTACCGTGCCCAGAACCACAACCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGc              <  1:1160791/139‑1 (MQ=255)                                                                                                                 ttGCCGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATAAATGcc             >  2:201711/1‑140 (MQ=255)                                                                                                                     cGAGGAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATAAATGCCGAAt         <  1:201711/140‑1 (MQ=255)                                                                                                                         gAATGCCGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCt     >  2:477740/1‑140 (MQ=255)                                                                                                                            tgccgcTAAAATCCCGGTGCGTCCTCCCGTCCCCAGACCCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGcc  <  1:460722/140‑1 (MQ=255)                                                                                                                               |                                                                                                                                         AAGAGTCGTTAATCGTCTTCTTCGTCATCCAGTTCAACGGGTGTCTGATACTGGTCAGGTTT‑AATGACCAGCAGGTCGCAGCGAAGATGATCAATCACCTGTTCCGCCGTGTTGCCGAGGAATGCTGCTGAAATACCGGTGCGTCCTACCGTGCCCAGAACCACAATCCCCGCCTGTAAGTGCTCCGCCAAATCAGGAATCACCTCTTCTGGCAGACCTTTTTCTACGTGCGTCATGTTTTCATTAATGCCGAATTTCTGCC  >  NC_000913/1397664‑1397925

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