Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,013,075	G→A	P389S (CCG→TCG)	yqeB ←	XdhC‑CoxI family protein YqeB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,013,075	0	G	A	100.0%	176.3 / NA	52	P389S (CCG→TCG)	yqeB	XdhC‑CoxI family protein YqeB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/28);  total (24/28)

GTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGGCCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCGGGATCTACCATCACAGGAATGAA  >  NC_000913/3012945‑3013210                                                                                                                                   |                                                                                                                                        gtATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAgcgc                                                                                                                                >  2:778755/1‑140 (MQ=255) gtATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAgcgc                                                                                                                                >  2:524359/1‑140 (MQ=255)           taataTTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACt                                                                                                                      <  1:778755/140‑1 (MQ=255)           taataTTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACt                                                                                                                      <  2:510335/140‑1 (MQ=255)           taataTTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCACGACCCAAGCGCGATTGTTACt                                                                                                                      <  2:864192/140‑1 (MQ=255)             atatTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTgg                                                                                                                    <  1:649132/140‑1 (MQ=255)               atTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTg                                                                                                                  <  1:186732/140‑1 (MQ=255)               atTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTg                                                                                                                  <  2:860156/140‑1 (MQ=255)                       gAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCt                                                                                                          >  1:376767/1‑140 (MQ=255)                             cGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCcgc                                                                                                    >  1:614502/1‑140 (MQ=255)                                       cTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAAt                                                                                          >  1:699869/1‑140 (MQ=255)                                       cTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAAt                                                                                          >  2:379083/1‑140 (MQ=255)                                         gcgcACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTc                                                                                        >  1:661914/1‑140 (MQ=255)                                          cgcACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCt                                                                                       >  1:213831/1‑140 (MQ=255)                                            cacaACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGt                                                                                     <  1:861370/140‑1 (MQ=255)                                            cacaACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGt                                                                                     <  1:280757/140‑1 (MQ=255)                                            cacaACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGt                                                                                     >  2:498673/1‑140 (MQ=255)                                               aaCCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTc                                                                                  >  2:1117805/1‑140 (MQ=255)                                                    gagTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAg                                                                             >  1:135456/1‑140 (MQ=255)                                                     agTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGa                                                                            <  2:135456/140‑1 (MQ=255)                                                                 gACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTccacc                                                                >  2:518095/1‑140 (MQ=255)                                                                    cGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGa                                                                            >  2:225789/1‑125 (MQ=255)                                                                    cGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGa                                                                            <  1:225789/125‑1 (MQ=255)                                                                      aGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCa                                                           >  1:1114964/1‑140 (MQ=255)                                                                       gCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCaa                                                          <  2:951880/140‑1 (MQ=255)                                                                        ccAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCATCCACCACGCaaa                                                         <  2:213831/140‑1 (MQ=255)                                                                         cAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAg                                                        <  1:524359/140‑1 (MQ=255)                                                                          aGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGc                                                       <  1:659646/140‑1 (MQ=255)                                                                            tGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCgg                                                     >  1:459299/1‑140 (MQ=255)                                                                            tGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCgg                                                     >  2:502921/1‑140 (MQ=255)                                                                              cccGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGtt                                                   <  1:393986/140‑1 (MQ=255)                                                                                 gcgATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCa                                                <  1:518095/140‑1 (MQ=255)                                                                                             aaTTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCa                                    <  2:336262/140‑1 (MQ=255)                                                                                               ttACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAAt                                  <  2:579164/140‑1 (MQ=255)                                                                                                  ccGCATGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGaa                               <  2:1114964/140‑1 (MQ=255)                                                                                                   cGCATGACAATCCTTACCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTtcatca                                        >  2:933604/1‑130 (MQ=255)                                                                                                   cGCATGACAATCCTTACCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTtcatca                                        <  1:933604/130‑1 (MQ=255)                                                                                                      aTGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAgg                           <  1:559147/140‑1 (MQ=255)                                                                                                       tGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGc                          >  1:286884/1‑140 (MQ=255)                                                                                                       tGACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGc                          <  2:661914/140‑1 (MQ=255)                                                                                                        gACAATCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCg                         <  1:502921/140‑1 (MQ=255)                                                                                                           aaTCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCGGGa                      <  2:614502/140‑1 (MQ=255)                                                                                                           aaTCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCGGGa                      <  2:459299/140‑1 (MQ=255)                                                                                                            aTCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCTCCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCGGGAt                     <  2:699869/140‑1 (MQ=255)                                                                                                            aTCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCGGGAt                     >  2:114914/1‑140 (MQ=255)                                                                                                             tCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCGGGAt                     >  2:333713/1‑139 (MQ=255)                                                                                                             tCCTTCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCGGGAt                     <  1:333713/139‑1 (MQ=255)                                                                                                                ttCCCAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCGGGATCTAc                 >  1:882110/1‑140 (MQ=255)                                                                                                                  cccAGCAGTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCGGGATCTACCa               >  2:1072044/1‑140 (MQ=255)                                                                                                                     agcagTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCgg                        >  2:1038107/1‑128 (MQ=255)                                                                                                                     agcagTAAAGCCCGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCgg                        <  1:1038107/128‑1 (MQ=255)                                                                                                                               cccGACCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCGGGATCTACCATCACAGgaatgaa  >  2:284275/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                        GTATGCCCCATAATATTGCCAGGAACACCGGTATTCTCCTGCGCACAACCAGAGTAAATCACCTGACCGAGCCAGTGCCCGCGATTTGTTTCAATTACCGCATGACAATCCTTCCCAGCAGTAAAGCCCGGCCCAAGCGCGATTGTTACTGGTGCCATATCTGCCCGCGTTCCCAAATTCTGTTTCGCCAGAATAGCGTCCACCACGCAAAGCGGTTTCAGTTCATCAAGCAATGAACAGGCGGGATCTACCATCACAGGAATGAA  >  NC_000913/3012945‑3013210

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