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

GCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAGCAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCA  >  NC_000913/2510752‑2511007                                                                                                                                       |                                                                                                                          gcAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACaaa                                                                                                                       >  2:925708/1‑139 (MQ=255)              ccAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCtt                                                                                                          <  2:659088/139‑1 (MQ=255)                aaTACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                        <  2:1046842/139‑1 (MQ=255)                 aTACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTAt                                                                                                       <  1:925708/139‑1 (MQ=255)                     cacaTTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGAtt                                                                                                   >  1:1166391/1‑139 (MQ=255)                        aTTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTaa                                                                                                >  1:175815/1‑139 (MQ=255)                            cTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                        >  1:929218/1‑127 (MQ=255)                            cTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                        <  2:929218/127‑1 (MQ=255)                            cTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGt                                                                                            >  1:748596/1‑139 (MQ=255)                            cTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGt                                                                                            >  2:451669/1‑139 (MQ=255)                              aCCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                <  2:231154/117‑1 (MQ=255)                              aCCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                >  1:231154/1‑117 (MQ=255)                                  ctgctCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGtt                                                                                           >  2:366260/1‑134 (MQ=255)                                  ctgctCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGtt                                                                                           <  1:366260/134‑1 (MQ=255)                                   tgctCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCg                                                                                     <  1:451669/139‑1 (MQ=255)                                     ctctGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCgat                                                                                   >  1:988786/1‑139 (MQ=255)                                          tATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                        <  1:294775/113‑1 (MQ=255)                                          tATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                        >  2:294775/1‑113 (MQ=255)                                          tATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATAc                                                                              <  2:903366/139‑1 (MQ=255)                                               tCATGCTTCCGGCATGGCTGTTATTAGCAGGTTAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCTCCCCACGATAACGTTTAACaaa                                                                                                                       >  2:1008565/1‑93 (MQ=255)                                               tCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACGGCCACCCCACGATAACGTTTAACaaa                                                                                                                       <  1:1008565/93‑1 (MQ=255)                                                    cTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTAc                                                                    >  1:568613/1‑139 (MQ=255)                                                      tCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCt                                                                  <  2:552281/139‑1 (MQ=255)                                                      tCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCt                                                                  <  2:175815/139‑1 (MQ=255)                                                       ccGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                <  1:511417/92‑1 (MQ=255)                                                       ccGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                >  2:511417/1‑92 (MQ=255)                                                                          cAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                <  2:342355/73‑1 (MQ=255)                                                                          cAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                >  1:342355/1‑73 (MQ=255)                                                                          cAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTAc                                                                    >  1:746616/1‑117 (MQ=255)                                                                          cAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTAc                                                                    <  2:746616/117‑1 (MQ=255)                                                                          cAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAAc                                              <  2:568613/139‑1 (MQ=255)                                                                             gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCga                                                                                    <  1:710611/98‑1 (MQ=255)                                                                             gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCga                                                                                    >  2:710611/1‑98 (MQ=255)                                                                             gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCga                                                                                    <  2:260191/98‑1 (MQ=255)                                                                             gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCga                                                                                    >  1:260191/1‑98 (MQ=255)                                                                             gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTcc                                                                >  2:1353714/1‑118 (MQ=255)                                                                             gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTcc                                                                <  1:1353714/118‑1 (MQ=255)                                                                             gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAgg                                           >  2:879250/1‑139 (MQ=255)                                                                             gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAgg                                           >  1:698991/1‑139 (MQ=255)                                                                                aGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGctc                                                                                                               >  2:1147023/1‑68 (MQ=255)                                                                                aGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGctc                                                                                                               <  1:1147023/68‑1 (MQ=255)                                                                                       gatgatgGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACgg                                 >  1:223015/1‑139 (MQ=255)                                                                                            tgGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCg                            >  1:1176341/1‑139 (MQ=255)                                                                                               tCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGTTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAata                                                        >  2:500191/1‑108 (MQ=255)                                                                                               tCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATa                                                                               >  1:43481/1‑85 (MQ=255)                                                                                               tCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATa                                                                               <  2:43481/85‑1 (MQ=255)                                                                                               tCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAata                                                        <  1:500191/108‑1 (MQ=255)                                                                                                  aTCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                        <  2:845790/57‑1 (MQ=255)                                                                                                  aTCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTa                                                                                                        >  1:845790/1‑57 (MQ=255)                                                                                                  aTCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCt                                                                  >  1:1225769/1‑95 (MQ=255)                                                                                                  aTCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCt                                                                  <  2:1225769/95‑1 (MQ=255)                                                                                                     gttcGAATCAACAGCCACTCCATGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGTTTTACCTCCCCATAAAATTGCGTTAACAGGCCGGTGACGGCTTCGAGCGTGAAg                   <  2:1158624/136‑1 (MQ=255)                                                                                                     gTCCGAAGCAACAGCCATCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGATATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGaac                   <  1:879250/139‑2 (MQ=255)                                                                                                     gTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAg                   <  2:1176341/139‑1 (MQ=255)                                                                                                              aaCAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAAc          <  2:223015/139‑1 (MQ=255)                                                                                                                cAGCCACCCCACGATAACGTTTAACAACAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATaa                                                                              <  1:878686/69‑2 (MQ=255)                                                                                                                cAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATAc                                                                              >  2:878686/1‑69 (MQ=255)                                                                                                                cAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACaa        >  1:70534/1‑139 (MQ=255)                                                                                                                    cACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATAc                                                                              <  1:105523/65‑1 (MQ=255)                                                                                                                    cACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATAc                                                                              >  2:105523/1‑65 (MQ=255)                                                                                                                     aCCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGa                     <  2:276672/121‑1 (MQ=255)                                                                                                                     aCCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGa                     >  1:276672/1‑121 (MQ=255)                                                                                                                      ccccACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTAt                                                                                         <  2:331777/52‑1 (MQ=255)                                                                                                                      ccccACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTAt                                                                                         >  1:331777/1‑52 (MQ=255)                                                                                                                      ccccACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCa  <  2:988786/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                          GCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAGCAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCA  >  NC_000913/2510752‑2511007

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: ‑