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

TTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAGCAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCAGGTCACTGTTTC  >  NC_000913/2510759‑2511019                                                                                                                                |                                                                                                                                      ttGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            <  1:1300027/139‑1 (MQ=255) ttGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            <  2:617351/139‑1 (MQ=255)           tACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCG‑AGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGa                                                                                                                 >  1:1148430/1‑139 (MQ=255)                caTTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTa                                                                                                             >  2:793219/1‑139 (MQ=255)                             gctCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCga                                                                                                <  1:793219/139‑1 (MQ=255)                                  gTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATa                                                                                           <  2:102853/139‑1 (MQ=255)                                  gTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCA‑CCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATAc                                                                                          <  1:182282/139‑1 (MQ=255)                                             cTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            >  2:358302/1‑95 (MQ=255)                                             cTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            <  1:358302/95‑1 (MQ=255)                                                    cATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCg                                                                         <  2:347135/139‑1 (MQ=255)                                                        gCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            <  2:992212/84‑1 (MQ=255)                                                        gCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            >  1:992212/1‑84 (MQ=255)                                                              catcatcggTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTaa                                                                                                            <  2:98270/87‑1 (MQ=255)                                                              catcatcggTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTaa                                                                                                            >  1:98270/8‑94 (MQ=255)                                                                  gCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            <  2:659378/74‑1 (MQ=255)                                                                  gCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            >  1:659378/1‑74 (MQ=255)                                                                  gCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTaa                                                           >  1:519210/1‑139 (MQ=255)                                                                      gTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAgg                                                       >  1:623957/1‑139 (MQ=255)                                                                        aaGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGcc                                                     >  2:359780/1‑139 (MQ=255)                                                                         aGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            >  1:1428150/1‑67 (MQ=255)                                                                         aGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            <  2:1428150/67‑1 (MQ=255)                                                                         aGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATa                                                                                           >  1:611571/1‑100 (MQ=255)                                                                         aGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATa                                                                                           <  2:611571/100‑1 (MQ=255)                                                                           ccGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGCTAACAGGCCGCt                                                  >  2:336491/1‑139 (MQ=255)                                                                             gatgatgatGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            >  1:44454/1‑63 (MQ=255)                                                                             gatgatgatGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGct                                                                                                                            <  2:44454/63‑1 (MQ=255)                                                                                    atgGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTc                                         >  1:365375/1‑139 (MQ=255)                                                                                     tgGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCg                                        >  1:278841/1‑139 (MQ=255)                                                                                      ggTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGa                                       <  2:365375/139‑1 (MQ=255)                                                                                        tCAATCGTCCG‑AGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCg                                    <  2:1148430/139‑1 (MQ=255)                                                                                          aaTCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGc                                                                                                                             <  1:770754/49‑1 (MQ=255)                                                                                          aaTCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGc                                                                                                                             >  2:770754/1‑49 (MQ=255)                                                                                           aTCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATAtt                                                                  <  2:570617/107‑1 (MQ=255)                                                                                           aTCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCATTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATAtt                                                                  >  1:570617/1‑107 (MQ=255)                                                                                             cGTCCGAAGCAACAGCCACCCCACGGTAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCt                                                                              <  1:940791/93‑1 (MQ=255)                                                                                             agTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCt                                                                              >  2:940791/2‑93 (MQ=255)                                                                                              gTCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAg                               <  1:359780/139‑1 (MQ=255)                                                                                               tCCGAAGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAgg                              >  2:530920/1‑139 (MQ=255)                                                                                                   aaGCAACAGCCACCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGaa                          <  2:54711/139‑1 (MQ=255)                                                                                                              aCCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGtt                                                                                                       >  1:811288/1‑51 (MQ=255)                                                                                                              aCCCCACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGtt                                                                                                       <  2:811288/51‑1 (MQ=255)                                                                                                                  cACGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGCTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCAGGt           <  1:336491/139‑1 (MQ=255)                                                                                                                    cGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTAc                                                                                >  2:136617/1‑68 (MQ=255)                                                                                                                    cGATAACGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTAc                                                                                <  1:136617/68‑1 (MQ=255)                                                                                                                          cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCtt                                          >  1:59159/1‑100 (MQ=255)                                                                                                                          cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCtt                                          <  2:59159/100‑1 (MQ=255)                                                                                                                          cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGa                                       >  1:1476276/1‑103 (MQ=255)                                                                                                                          cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGa                                       <  2:1476276/103‑1 (MQ=255)                                                                                                                          cGTTTAACAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGG‑GAATAACAATGTTCAGGTCACTGTTTc  >  1:1395114/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                      TTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCAATCGTCCGAAGCAACAGCCACCCCACGATAACGTTTAGCAAAAAATGCTCTCCTTTATGATTTAAGAGTTATGGCGATGATACGGGGCTTTACCTCCCCGTAATATTGCGTTAACAGGCCGCTGACGGCTTCGAGCGTGAAGGAGAATAACAATGTTCAGGTCACTGTTTC  >  NC_000913/2510759‑2511019

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