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%	186.2 / NA	58	*419* (TAG→TAA)	yfeO	putative ion channel protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/28);  total (30/28)

GGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAGCAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCAGGTCACTGTTTCT  >  NC_000913/2510751‑2511020                                                                                                                                        |                                                                                                                                       ggcAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACaa                                                                                                                                     >  1:1081576/1‑139 (MQ=255)         ttGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             <  2:281278/139‑1 (MQ=255)          tGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGAGGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGctc                                                                                                                            <  1:705927/139‑1 (MQ=255)              cccAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCt                                                                                                                        >  2:327124/1‑139 (MQ=255)                cAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCttt                                                                                                                      >  1:145376/1‑139 (MQ=255)                             cTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGt                                                                                                         >  2:406171/1‑139 (MQ=255)                               aCCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTa                                                                                                       >  1:514543/1‑139 (MQ=255)                                     gctCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCga                                                                                                 >  2:656948/1‑139 (MQ=255)                                          gTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATa                                                                                            <  1:898170/139‑1 (MQ=255)                                                 cATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             <  1:878355/99‑1 (MQ=255)                                                 cATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             >  2:878355/1‑99 (MQ=255)                                                     cTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             >  1:1063921/1‑95 (MQ=255)                                                     cTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             <  2:1063921/95‑1 (MQ=255)                                                     cTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTAc                                                                                 >  1:1166442/1‑139 (MQ=255)                                                     cTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTAc                                                                                 >  2:1098421/1‑139 (MQ=255)                                                       tCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCt                                                                               <  1:406171/139‑1 (MQ=255)                                                       tCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCt                                                                               <  2:1081576/139‑1 (MQ=255)                                                       tCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCt                                                                               <  1:1060504/139‑1 (MQ=255)                                                        ccGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTc                                                                              >  1:424998/1‑139 (MQ=255)                                                         cGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTcc                                                                             <  1:327124/139‑1 (MQ=255)                                                              tGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGtt                                                                                                        <  1:132565/107‑1 (MQ=255)                                                              tGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGtt                                                                                                        >  2:132565/1‑107 (MQ=255)                                                                gCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGctc                                                                                                                            <  2:707950/85‑1 (MQ=255)                                                                gCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGctc                                                                                                                            >  1:707950/1‑85 (MQ=255)                                                                gCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAat                                                                      >  1:196888/1‑139 (MQ=255)                                                                           cAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             >  2:502566/1‑73 (MQ=255)                                                                           cAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                             <  1:502566/73‑1 (MQ=255)                                                                              gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATg                                                                                                                               >  2:692904/1‑68 (MQ=255)                                                                              gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATg                                                                                                                               <  1:692904/68‑1 (MQ=255)                                                                                 aGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                                     >  2:678532/1‑75 (MQ=255)                                                                                 aGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                                     <  1:678532/75‑1 (MQ=255)                                                                                 aGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCg                                                     >  2:658377/1‑139 (MQ=255)                                                                                             tgGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                                     >  2:720246/1‑63 (MQ=255)                                                                                             tgGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                                     <  1:720246/63‑1 (MQ=255)                                                                                             tgGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTAt                                                                                                                    >  2:629975/1‑64 (MQ=255)                                                                                             tgGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTAt                                                                                                                    <  1:629975/64‑1 (MQ=255)                                                                                                   aTCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTg                                                  >  1:529809/1‑124 (MQ=255)                                                                                                   aTCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTg                                                  <  2:529809/124‑1 (MQ=255)                                                                                                    tCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGa                                  <  2:514543/139‑1 (MQ=255)                                                                                                       tCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAgg                               >  2:616294/1‑139 (MQ=255)                                                                                                        ccGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGga                              <  1:616294/139‑1 (MQ=255)                                                                                                         cGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                                     <  2:1141655/51‑1 (MQ=255)                                                                                                         cGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                                     >  1:1141655/1‑51 (MQ=255)                                                                                                         cGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGgag                             <  2:1166442/139‑1 (MQ=255)                                                                                                                      aCCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATa                                                                                            >  2:1173602/1‑63 (MQ=255)                                                                                                                      aCCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATa                                                                                            <  1:1173602/63‑1 (MQ=255)                                                                                                                           aCGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCAGGTc           <  1:656948/139‑1 (MQ=255)                                                                                                                              aTAGCGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTACTATTGCGTTAACAGGCCACTGACGGCTTCGAGCg                                     >  1:783464/1‑110 (MQ=255)                                                                                                                              aTAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCg                                     <  2:783464/110‑1 (MQ=255)                                                                                                                                  cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATAc                                                                                           <  2:675904/52‑1 (MQ=255)                                                                                                                                  cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTcc                                                                             >  2:152552/1‑66 (MQ=255)                                                                                                                                  cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTcc                                                                             <  1:152552/66‑1 (MQ=255)                                                                                                                                  cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAat                                                                      <  2:210396/73‑1 (MQ=255)                                                                                                                                  cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAat                                                                      >  1:210396/1‑73 (MQ=255)                                                                                                                                  cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTg                                                                  >  2:475487/1‑77 (MQ=255)                                                                                                                                  cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTg                                                                  <  1:475487/77‑1 (MQ=255)                                                                                                                                  cGTTTAACAAAAAATGCTATCCTTTATGATTTAAGAGTTATGGCGATGATAc                                                                                           >  1:675904/1‑52 (MQ=255)                                                                                                                                    tttAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCAGGTCACTGTTTCt  >  1:260503/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                       GGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAGCAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCAGGTCACTGTTTCT  >  NC_000913/2510751‑2511020

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