Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,033,155	G→T	R4L (CGC→CTC)	trkH →	K(+) transporter TrkH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,033,155	0	G	T	100.0%	148.1 / NA	44	R4L (CGC→CTC)	trkH	K(+) transporter TrkH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (18/26);  total (18/26)

CGCTTCCGGCGGCGAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATT  >  NC_000913/4033018‑4033291                                                                                                                                          |                                                                                                                                         cgcTTCCGGCGGCGAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTctcg                                                                                                                                        >  1:774568/1‑140 (MQ=255) cgcTTCCGGCGGCGAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTctcg                                                                                                                                        >  1:168581/1‑140 (MQ=255)   cTTCCGGCGGGGAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGcc                                                                                                                                      <  2:27291/140‑1 (MQ=255)      ccGGCGGCGAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCAtt                                                                                                                                   >  1:618210/1‑140 (MQ=255)           ggcgAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCg                                                                                                                              <  2:966546/140‑1 (MQ=255)            gcgAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGa                                                                                                                             <  2:73843/140‑1 (MQ=255)              gAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAAt                                                                                                                           <  2:618210/140‑1 (MQ=255)              gAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAAt                                                                                                                           <  1:302978/140‑1 (MQ=255)              gAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAAt                                                                                                                           <  1:515570/140‑1 (MQ=255)              gAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAAt                                                                                                                           <  1:83827/140‑1 (MQ=255)                aaGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTGCCCGAATCg                                                                                                                         <  2:168581/140‑1 (MQ=255)                aaGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCg                                                                                                                         <  2:774568/140‑1 (MQ=255)                             ttCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGa                                                                                                                             >  2:937538/1‑123 (MQ=255)                             ttCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGa                                                                                                                             <  1:937538/123‑1 (MQ=255)                             ttCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGt                                                                                                            <  2:795154/140‑1 (MQ=255)                              tCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTc                                                                                                           >  2:392939/1‑140 (MQ=255)                                               ttttAGGCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGAc                                                                                          <  1:109033/140‑1 (MQ=255)                                               ttttAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGAc                                                                                          >  2:201982/1‑140 (MQ=255)                                                     ccGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGTACTACTGGTCATCTTATTTTCAGGGACCATGAt                                                                                    <  2:729820/140‑1 (MQ=255)                                                     ccGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGAt                                                                                    <  1:724977/140‑1 (MQ=255)                                                           ttCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATccc                                                                              <  2:274213/140‑1 (MQ=255)                                                            tCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCt                                                                             <  1:392939/140‑1 (MQ=255)                                                              aTTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTgg                                                                           <  2:537419/140‑1 (MQ=255)                                                               ttGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTggg                                                                          <  1:338233/140‑1 (MQ=255)                                                                     ttGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTa                                                                    >  2:152575/1‑140 (MQ=255)                                                                      tGTTAGCGATTGAAGAATAATCCCCCCTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAg                                                                   <  2:369738/140‑1 (MQ=255)                                                                          aGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACt                                                               >  2:679200/1‑140 (MQ=255)                                                                           gCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTc                                                              <  1:305144/140‑1 (MQ=255)                                                                                            cccACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGc                                             >  2:163261/1‑140 (MQ=255)                                                                                                  tCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCcg                                       <  1:304202/140‑1 (MQ=255)                                                                                                               cTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCcg                                       >  1:940978/1‑127 (MQ=255)                                                                                                               cTAAGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCcg                                       <  2:940978/127‑1 (MQ=255)                                                                                                                  aGGAAGCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGAcc                       >  2:932776/1‑140 (MQ=255)                                                                                                                       gCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCttttt                  >  2:116392/1‑140 (MQ=255)                                                                                                                       gCGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCttttt                  >  1:733932/1‑140 (MQ=255)                                                                                                                        cGGCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCtttttt                 >  2:340868/1‑140 (MQ=255)                                                                                                                         ggCAGAGATGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTg                >  1:21060/1‑140 (MQ=255)                                                                                                                              agaTGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGTACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGcc           >  2:642731/1‑140 (MQ=255)                                                                                                                              agaTGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGcc           >  1:134665/1‑140 (MQ=255)                                                                                                                                 tGCATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGTACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTc        <  1:642731/140‑1 (MQ=255)                                                                                                                                   cATTTTCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGc      >  1:885869/1‑140 (MQ=255)                                                                                                                                     ttttCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCa    >  1:246090/1‑140 (MQ=255)                                                                                                                                     ttttCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCa    >  1:181212/1‑140 (MQ=255)                                                                                                                                      tttCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCAt   <  1:445154/140‑1 (MQ=255)                                                                                                                                      tttCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCAt   <  2:885869/140‑1 (MQ=255)                                                                                                                                       ttCTCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCAtt  <  1:201982/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         CGCTTCCGGCGGCGAAAGTGGCTGAATTTTCCGCAAAGCTGGCGGATTTTAGCCGTGGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATT  >  NC_000913/4033018‑4033291

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