Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,510,886	G→A	*419* (TAG→TAA)	yfeO →	putative ion channel protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,510,886	0	G	A	100.0%	128.2 / NA	41	*419* (TAG→TAA)	yfeO	putative ion channel protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/19);  total (22/19)

TTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAGCAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCAGGTCACTGTTTCT  >  NC_000913/2510759‑2511020                                                                                                                                |                                                                                                                                       ttGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             <  2:762752/139‑1 (MQ=255)        cAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCttt                                                                                                                      >  1:124732/1‑139 (MQ=255)                caTTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                                     >  2:330868/1‑132 (MQ=255)                caTTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                                     <  1:330868/132‑1 (MQ=255)                caTTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTa                                                                                                              >  2:671721/1‑139 (MQ=255)                     cTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGt                                                                                                         >  2:603832/1‑139 (MQ=255)                       aCCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTa                                                                                                       >  1:402099/1‑139 (MQ=255)                             gctCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCga                                                                                                 <  1:128839/139‑1 (MQ=255)                              ctctGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCgat                                                                                                >  2:384323/1‑139 (MQ=255)                                  gTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATa                                                                                            <  1:603832/139‑1 (MQ=255)                                        tCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             >  2:551946/1‑100 (MQ=255)                                        tCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             <  1:551946/100‑1 (MQ=255)                                           tGCTTCCGGCATGGCTGTTATTTGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGACACCCCACGATAACGTTTTACAAAACATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCttt                                                                                   >  1:663290/1‑139 (MQ=255)                                             cTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTAc                                                                                 >  1:748214/1‑139 (MQ=255)                                             cTTACGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAACAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTAc                                                                                 >  1:723071/1‑139 (MQ=255)                                               tCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCt                                                                               <  2:127105/139‑1 (MQ=255)                                                     aTGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATAc                                                                                           >  2:480763/1‑121 (MQ=255)                                                     aTGGCTGTTATTAGCAGGTAAGCAGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATAc                                                                                           <  1:480763/121‑1 (MQ=255)                                                       ggCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAGAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCTCCGTaa                                                                       <  1:671721/139‑1 (MQ=255)                                                        gCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAAAGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAat                                                                      >  1:354448/1‑139 (MQ=255)                                                                      gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             <  1:602558/70‑1 (MQ=255)                                                                      gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             >  2:602558/1‑70 (MQ=255)                                                                      gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATAc                                                                                           <  2:448233/104‑1 (MQ=255)                                                                      gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATAc                                                                                           >  1:448233/1‑104 (MQ=255)                                                                      gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAgg                                                        <  2:402099/139‑1 (MQ=255)                                                                      gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAgg                                                        >  2:206189/1‑139 (MQ=255)                                                                                     tgGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCg                                         >  1:775061/1‑139 (MQ=255)                                                                                      ggTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTAc                                                                                 >  2:147647/1‑98 (MQ=255)                                                                                      ggTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTAc                                                                                 <  1:147647/98‑1 (MQ=255)                                                                                              gTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAg                                <  2:354448/139‑1 (MQ=255)                                                                                              gTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAg                                <  2:124732/139‑1 (MQ=255)                                                                                              gTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAg                                <  1:384323/139‑1 (MQ=255)                                                                                              gTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAg                                <  1:206189/139‑1 (MQ=255)                                                                                                 cGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGgag                             >  2:220027/1‑139 (MQ=255)                                                                                                              aCCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTc                >  1:149831/1‑139 (MQ=255)                                                                                                                    cGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCAGGTCa          <  2:149831/139‑1 (MQ=255)                                                                                                                      aTAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGaa                           <  2:258350/120‑1 (MQ=255)                                                                                                                      aTAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGaa                           >  1:258350/1‑120 (MQ=255)                                                                                                                          cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTcccc                                                                           <  1:783696/68‑1 (MQ=255)                                                                                                                          cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTcccc                                                                           >  2:783696/1‑68 (MQ=255)                                                                                                                            tttAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCAGGTCACTGTTTCt  >  1:581487/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                       TTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAGCAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCAGGTCACTGTTTCT  >  NC_000913/2510759‑2511020

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