Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	4,111,596	G→A	L179L (CTG→TTG)	yiiQ ←	DUF1454 family putative periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	4,111,596	0	G	A	100.0%	104.9 / NA	34	L179L (CTG→TTG)	yiiQ	DUF1454 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/13);  total (21/13)

AGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCG  >  NC_000913_3_pae_tpiA/4111453‑4111737                                                                                                                                                |                                                                                                                                              aGGAAGAGGGCCGCCCCGCAGGGGAGCAGGTCTGTGAAACAGTATAGAGATTCATTGGCACAAAGGCTTTGCTTTTTGCCATTTATTCAAACCTTCGAGCGATTCAGATAGGGCCAGCTTAATCTGTTCAACAGCGAAGGTCAACCCCt                                                                                                                                          >  1:586978/1‑149 (MQ=255)         ggCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGTAGGTCAACCCCTTTTCGCCg                                                                                                                                  <  2:362638/149‑1 (MQ=255)                ccGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTACAACAGCGACGGTCAACCCCTTTTCGCCGTTGTCcg                                                                                                                           >  2:244924/1‑149 (MQ=255)                        cagcagGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGacaaca                                                                                                                   <  1:616350/149‑1 (MQ=255)                                         aGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTc                                                                                                  >  1:7727/1‑149 (MQ=255)                                          gTATAGAGATTCATCGTCACAAAGGCTTTGCTTTTTGGCATTTATTCAAACCTTCAAGCAATTTAGATAGCGCCAGCATAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGCCAACATAACGCAGTGCACCTTcc                                                                                                 <  2:586978/149‑2 (MQ=255)                                           tataGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCTACCCCTTTTCGCCGTTGTCCGCGCCAACATAACGCAGTGCACCTTCTg                                                                                                >  2:93293/1‑149 (MQ=255)                                               gagaTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtctc                                                                                            <  2:536484/149‑1 (MQ=255)                                               gagaTTCAGCGGCACAAAGGCTTTGCGTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCCACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGtctc                                                                                            <  1:501246/149‑1 (MQ=255)                                                  aTTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                         >  1:303995/1‑149 (MQ=255)                                                  aTTCACGGGCACAAAGGCTTTGCTTTTTGTCATTTTTGCAAATCGTCAAGGGATTCAGATAGCGCCAGCTTAAGCGGTTCAACAGCGAAGGTCAACCCCTTTGCGGCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                         <  2:65662/149‑1 (MQ=255)                                                          gCACAAATGCCTAGCTTTTTGTCATTTATACAAACCTTCAAGCGATTCTTATATCACCAGCTTAATCGGTTCAACAGCGAAGGTCAACACCGTTTCGAAATTTTCTGCTACAACATAACGCCGTGCACCTTCTGTCTCGGTGTAATAAc                                                                                 <  2:163633/149‑1 (MQ=255)                                                              aaaGGATTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTACTAACgttt                                                                             >  2:551918/1‑149 (MQ=255)                                                                       gctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTccc                                                                    >  2:468576/1‑149 (MQ=255)                                                                               tCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTta                                                            >  1:264232/1‑149 (MQ=255)                                                                                 aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                         >  1:14354/1‑118 (MQ=255)                                                                                 aTTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGgt                                                                                         <  2:14354/118‑1 (MQ=255)                                                                                    tATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGc                                                       <  1:42931/149‑1 (MQ=255)                                                                                    tATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGc                                                       <  1:43149/149‑1 (MQ=255)                                                                                       tCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                    >  1:402760/1‑149 (MQ=255)                                                                                         aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAAGGTTTGTTTTTCCCCGCCGTTCGTAGCGACTg                                                  >  2:332557/1‑149 (MQ=255)                                                                                         aaaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTg                                                  >  1:509404/1‑149 (MQ=255)                                                                                          aaCCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGGCGTTAGTCGCGACTga                                                 >  2:205644/1‑148 (MQ=255)                                                                                             cTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGt                                              >  1:501747/1‑149 (MQ=255)                                                                                              ttCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt                                             >  1:487821/1‑149 (MQ=255)                                                                                                 aaGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTACCCGCCGTTAGTAGCGACTGCAGttttt                                          >  2:340488/1‑149 (MQ=255)                                                                                                           gATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTTTTTTTCCCCGCCGTTACTCGCGACTGCAGTTTTTTCTCGCTTtg                                >  2:84678/1‑149 (MQ=255)                                                                                                           gATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtg                                >  1:496445/1‑149 (MQ=255)                                                                                                            aTAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgt                               >  1:458269/1‑149 (MQ=255)                                                                                                                  ccAGCTTAATCGGTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTg                         <  1:468576/149‑1 (MQ=255)                                                                                                                              gTTCAACAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg             <  2:487821/149‑1 (MQ=255)                                                                                                                                  aaCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGCCTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGa         >  1:192587/1‑149 (MQ=255)                                                                                                                                  aaCAGCGAAGGTCAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTACTGGGGTGa         >  2:185613/1‑149 (MQ=255)                                                                                                                                         aaGGTAAACCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCg  <  2:509404/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                              AGGAAGATGGCCGCCCCGCAGGGCAGCAGGTCTGTGAAACAGTATAGAGATTCATCGGCACAAAGGCTTTGCTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCG  >  NC_000913_3_pae_tpiA/4111453‑4111737

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑