Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,559,878	C→T	A41A (GCG→GCA)	yjiG ←	Gate family protein YjiG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,559,878	0	C	T	100.0%	136.7 / NA	39	A41A (GCG→GCA)	yjiG	Gate family protein YjiG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/18);  total (21/18)

TTCACCGGGTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGCGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAGCCACGGCGCGCGC  >  NC_000913/4559800‑4559954                                                                               |                                                                             ttCACCGGGTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAgtgc                                                                            >  2:840802/1‑81 (MQ=255)   cACCGGGTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAgtgc                                                                            <  1:840742/79‑1 (MQ=255)   cACCGGGTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCt                                                                          >  2:840793/1‑81 (MQ=255)   cACCGGGTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCt                                                                          >  2:840796/1‑81 (MQ=255)    aCCGGGTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCt                                                                          <  2:840745/80‑1 (MQ=255)    aCCGGGTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCt                                                                          >  2:840794/1‑80 (MQ=255)    aCCGGGTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCt                                                                          <  1:840744/80‑1 (MQ=255)    aCCGGGTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTg                                                                         <  2:840743/81‑1 (MQ=255)     ccGGGTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTg                                                                         >  1:840806/1‑80 (MQ=255)       gggTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAat                                                                      <  1:840746/81‑1 (MQ=255)            aTCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAAtgatga                                                                 >  1:840801/1‑81 (MQ=255)               cccACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACaa                                                              <  2:840747/81‑1 (MQ=255)               cccACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACa                                                               >  1:840812/1‑80 (MQ=255)                ccACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACaaa                                                             <  2:840748/81‑1 (MQ=255)                 cacaGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAAc                                                            <  1:840749/81‑1 (MQ=255)                    aGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGcc                                                         <  1:840750/81‑1 (MQ=255)                            tACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATcaccac                                                 >  2:840810/1‑81 (MQ=255)                               gggTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGtt                                              >  1:840807/1‑81 (MQ=255)                                 gTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTg                                            >  2:840804/1‑81 (MQ=255)                                     acaGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGgcag                                        >  1:840817/1‑81 (MQ=255)                                        gATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGcagcag                                     <  2:840751/81‑1 (MQ=255)                                            tGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGtt                                  <  1:840752/80‑1 (MQ=255)                                            tGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTg                                 <  1:840753/81‑1 (MQ=255)                                               cccacccaGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTc                              <  1:840754/81‑1 (MQ=255)                                                  acccaGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTg                           >  1:840830/1‑81 (MQ=255)                                                    ccaGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGc                         >  2:840811/1‑81 (MQ=255)                                                    ccaGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGc                         <  2:840756/81‑1 (MQ=255)                                                    ccaGTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGc                         <  1:840755/81‑1 (MQ=255)                                                       gTCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGa                       >  1:840818/1‑80 (MQ=255)                                                        tCGAGCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGa                       >  1:840809/1‑79 (MQ=255)                                                            gCAGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTa                  >  1:840823/1‑80 (MQ=255)                                                              aGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAg               >  1:840820/1‑81 (MQ=255)                                                              aGGCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAg               >  1:840814/1‑81 (MQ=255)                                                               ggCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAGc              <  2:840758/81‑1 (MQ=255)                                                               ggCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAGc              <  1:840757/81‑1 (MQ=255)                                                                gCCGGTAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAGc              >  2:840829/1‑80 (MQ=255)                                                                     tAATTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAGCCACGgc        <  1:840759/81‑1 (MQ=255)                                                                       aTTTTCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAGCCACGgcgc      >  2:840828/1‑81 (MQ=255)                                                                        ttttCAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAGCCACGgcgcg     >  1:840813/1‑81 (MQ=255)                                                                            cAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAGCCACGgcgcgcgc  <  1:840760/80‑1 (MQ=255)                                                                            cAGTGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAGCCACGgcgcgcgc  >  2:840819/1‑80 (MQ=255)                                                                               |                                                                             TTCACCGGGTAATCCCCACAGCGCCATTACGGGTTCACAGATATGCCCCACCCAGTCGAGCAGGCCGGTAATTTTCAGCGCCTGAATGATGACAAACGCCATCACCACGTTTGGCAGCAGGTTGGTCGTGGCGATAGTAAAGCCACGGCGCGCGC  >  NC_000913/4559800‑4559954

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