Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,033,689	C→A	A182E (GCG→GAG)	trkH →	K(+) transporter TrkH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,033,689	0	C	A	100.0%	122.6 / NA	35	A182E (GCG→GAG)	trkH	K(+) transporter TrkH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/20);  total (15/20)

TTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCTTTGCGACTATCGCTATTGGCGGCTTCTCGACACAT  >  NC_000913/4033553‑4033822                                                                                                                                         |                                                                                                                                       tttGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGaa                                                                                                                                     <  1:369877/140‑1 (MQ=255)   tGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGaaaa                                                                                                                                   >  2:83355/1‑140 (MQ=255)          aTGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAAACCtgtg                                                                                                                            >  1:198149/1‑140 (MQ=255)           tGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTgg                                                                                                                           <  1:298218/140‑1 (MQ=255)                   atcGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTa                                                                                                                   >  2:234518/1‑140 (MQ=255)                       tgtTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTc                                                                                                               <  2:392202/140‑1 (MQ=255)                       tgtTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTc                                                                                                               <  1:83355/140‑1 (MQ=255)                           ggCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGc                                                                                                           <  1:305661/140‑1 (MQ=255)                                     aTACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑Gc                                                                                                <  1:320757/140‑1 (MQ=255)                                        cTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCtg                                                                                             <  1:234518/140‑1 (MQ=255)                                            cTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTgcg                                                                                         <  1:157950/140‑1 (MQ=255)                                             tATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTgcgc                                                                                        <  2:384568/140‑1 (MQ=255)                                                   cGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGt                                                                                  <  1:253538/140‑1 (MQ=255)                                                    ggCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGtt                                                                                 >  2:244779/1‑140 (MQ=255)                                                     gCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTtg                                                                                >  1:174158/1‑140 (MQ=255)                                                     gCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTtg                                                                                >  2:453708/1‑140 (MQ=255)                                                         gggtgggATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGt                                                                            >  1:352083/1‑140 (MQ=255)                                                         gggtgggATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGt                                                                            <  1:453708/140‑1 (MQ=255)                                                          ggtgggATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGtt                                                                           >  2:235994/1‑140 (MQ=255)                                                          ggggggATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGGTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGtt                                                                           >  2:439461/4‑140 (MQ=255)                                                                      ctctATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGa                                                               <  1:244779/140‑1 (MQ=255)                                                                         tATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGc                                                            >  1:36268/1‑140 (MQ=255)                                                                                    aaTGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCa                                                 <  2:136178/140‑1 (MQ=255)                                                                                    aaTGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCa                                                 <  2:352083/140‑1 (MQ=255)                                                                                        cccggcccgCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCgg                                             >  2:128793/1‑140 (MQ=255)                                                                                             cccgCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATa                                        >  2:296372/1‑140 (MQ=255)                                                                                                 cTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCtt                                    >  2:147557/1‑140 (MQ=255)                                                                                                  tGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCttt                                   >  2:261851/1‑140 (MQ=255)                                                                                                          aaCAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCTTTGCGACTAt                           <  2:459579/140‑1 (MQ=255)                                                                                                            cAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCTTTGCGACTATCg                         >  1:347779/1‑140 (MQ=255)                                                                                                            cAAAATGCGCCCGCGAATTGCGGAAACGGAGAAAACCCTGGGGTTGATTTATGGCTTGCTGGCCGGCGGCCGGTG‑GCCGGCGTTGGGGTTTGGCGGAATTGGTGCCTTTGATGCCATCGGGCATAGCTTTTGGACttacc                         >  2:142735/1‑136 (MQ=255)                                                                                                                  gcgcCCGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCTTTGCGACTATCGCTATTg                   >  2:30161/1‑140 (MQ=255)                                                                                                                      ccGCGAATTGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCTTTGCGACTATCGCTATTggcgg               <  1:235994/140‑1 (MQ=255)                                                                                                                              tGCGGAAACGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGAGCAGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCTTTGCGACTATCGCTATTGGCGGCTTCTCGa       <  2:172656/140‑1 (MQ=255)                                                                                                                                   aaaCGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCTTTGCGACTATCGCTATTGGCGGCTTCTCGACACat  <  2:347779/140‑1 (MQ=255)                                                                                                                                   aaaCGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCTTTGCGACTATCGCTATTGGCGGCTTCTCGACACat  <  2:198149/140‑1 (MQ=255)                                                                                                                                   aaaCGGAGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCTTTGCGACTATCGCTATTGGCGGCTTCTCGACACat  <  1:261851/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       TTTGGCGGGATGGGGATCATCGTGTTGGCGGTTGCGATACTGCCTATCCTCGGCGTGGGTGGGATGCAGCTCTATCGCGCAGAAATGCCCGGCCCGCTGAAAGATAACAAAATGCGCCCGCGAATTGCGGAAACGGCGAAAACCCTGTGGTTGATTTATGTCTTGCTGACCGTC‑GCCTGTGCGCTGGCGTTGTGGTTTGCCGGAATGGATGCCTTTGATGCCATCGGCCATAGCTTTGCGACTATCGCTATTGGCGGCTTCTCGACACAT  >  NC_000913/4033553‑4033822

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