Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,033,217	G→A	G25R (GGG→AGG)	trkH →	potassium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,033,217	0	G	A	100.0%	187.1 / NA	48	G25R (GGG→AGG)	trkH	potassium transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/27);  total (21/27)

GTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAGGGGTTTCTGATAGTGG  >  NC_000913/4033072‑4033364                                                                                                                                                  |                                                                                                                                                    gTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTg                                                                                                                                                >  2:511815/1‑151 (MQ=255)    gTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTa                                                                                                                                             <  1:15769/151‑1 (MQ=255)    gTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTa                                                                                                                                             <  2:613162/151‑1 (MQ=255)       cATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGc                                                                                                                                           <  2:1004130/150‑1 (MQ=255)           gCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCa                                                                                                                                      <  1:581379/151‑1 (MQ=255)                   tAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCgg                                                                                                                              >  2:1002680/1‑151 (MQ=255)                    aGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCgg                                                                                                                              <  1:917051/150‑1 (MQ=255)                       gATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATg                                                                                                                          >  1:237935/1‑151 (MQ=255)                       gATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATg                                                                                                                          <  2:317020/151‑1 (MQ=255)                          tGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGa                                                                                                                        <  2:661625/150‑1 (MQ=255)                           gaagaaTAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGc                                                                                                                      <  2:630357/151‑1 (MQ=255)                           gaagaaTAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGc                                                                                                                      >  2:361935/1‑151 (MQ=255)                                  aatCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCcgc                                                                                                               >  1:135924/1‑151 (MQ=255)                                  aatCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCcgc                                                                                                               <  1:828776/151‑1 (MQ=255)                                            tCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTAc                                                                                                       <  2:432154/149‑1 (MQ=255)                                               ttttGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTAccc                                                                                                     >  2:661059/1‑148 (MQ=255)                                                     aGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCtttt                                                                                            >  1:569395/1‑151 (MQ=255)                                                           aaGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCg                                                                                      <  2:596259/151‑1 (MQ=255)                                                           aaGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCAGGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGt                                                                                        <  2:678047/149‑1 (MQ=255)                                                                      agagATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCAtt                                                                           >  2:71811/1‑151 (MQ=255)                                                                          aTGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGct                                                                       <  1:491214/151‑1 (MQ=255)                                                                          aTGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGc                                                                        >  1:952614/1‑150 (MQ=255)                                                                            gCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGctct                                                                     >  1:243152/1‑151 (MQ=255)                                                                             cATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTa                                                                    <  2:560485/151‑1 (MQ=255)                                                                              aTTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTAt                                                                   <  1:890649/151‑1 (MQ=255)                                                                                 ttCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCt                                                                >  2:273757/1‑151 (MQ=255)                                                                                    gcgcCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCtgtg                                                             >  1:187771/1‑151 (MQ=255)                                                                                            aCCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCg                                                      >  1:879251/1‑150 (MQ=255)                                                                                            aCCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGa                                                     <  1:951456/151‑1 (MQ=255)                                                                                              ccGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCTTTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAAc                                                   >  2:641794/1‑151 (MQ=255)                                                                                                gAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCg                                                 <  1:673154/151‑1 (MQ=255)                                                                                                gAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCg                                                 >  1:37273/1‑151 (MQ=255)                                                                                                       tGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAaga                                           <  1:557101/150‑1 (MQ=255)                                                                                                             actGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGc                                    <  1:287464/151‑1 (MQ=255)                                                                                                                        tATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTa                             <  2:92881/147‑1 (MQ=255)                                                                                                                        tATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTa                             <  2:925648/147‑1 (MQ=255)                                                                                                                         aTTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATcc                        >  1:714789/1‑151 (MQ=255)                                                                                                                         aTTTTCAGGGACCATGATTATCCCTAGGCTGGTAACACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGTCCTTGTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAAt                          <  1:256064/149‑1 (MQ=255)                                                                                                                          ttttCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATccc                       >  2:945484/1‑151 (MQ=255)                                                                                                                          ttttCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATccc                       >  2:930892/1‑151 (MQ=255)                                                                                                                           tttCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATccc                       <  2:578610/150‑1 (MQ=255)                                                                                                                              cAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTg                    >  2:314756/1‑150 (MQ=255)                                                                                                                              cAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGa                   <  2:159645/151‑1 (MQ=255)                                                                                                                                gggACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAgg                 <  1:450337/151‑1 (MQ=255)                                                                                                                                  gACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAgggg               <  2:610246/151‑1 (MQ=255)                                                                                                                                            aTCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAGGGGTTTCTGATa      <  1:749422/150‑1 (MQ=255)                                                                                                                                             tCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAGGGGTTTCTGATAg     >  1:225962/1‑150 (MQ=255)                                                                                                                                               ccTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAGGGGTTTCTGATAgtgg  >  2:684459/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAGGGGTTTCTGATAGTGG  >  NC_000913/4033072‑4033364

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