Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,436,697	G→A	E60E (GAG→GAA)	trkA →	NAD‑binding component of Trk potassium transporters

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,436,697	0	G	A	96.0%	142.1 / ‑1.8	51	E60E (GAG→GAA)	trkA	NAD‑binding component of Trk potassium transporters
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (23/25);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

GACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGGCGCATT  >  NC_000913/3436461‑3436945                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         gACGGCTTCTTTTACGCAAAGCTAATCAAAAAGCGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGGTACCCACGGTGAGCGTCTGCGGACCTTACAGGGTAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCCTCCACGCGGATTGCGGGAAGCAggggacga                                                                                                                                                                                                                                               >  2:29119/1‑242 (MQ=255)  aCGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCGCATCCACGCGTATTGCGGGAAGCAGGTGCcgtcg                                                                                                                                                                                                                                             <  2:63992/249‑4 (MQ=255)  aCGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTAAAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCcgacg                                                                                                                                                                                                                                             >  1:63992/1‑249 (MQ=255)    ggCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGc                                                                                                                                                                                                                                         >  1:187799/1‑251 (MQ=255)     gCTTCTTTTACGCAAAGACAATCAAAAAGTGATGAGATAACGGGTTGCGACTGCTGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGCAACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGAGCCGACGACGcc                                                                                                                                                                                                                                        >  1:335174/1‑251 (MQ=255)         cttTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGat                                                                                                                                                                                                                                     >  1:295127/1‑250 (MQ=255)            ttACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGAAGAAAATTATCATTCTGGGAGCCGGCCAGGTTGGCGGAACACTGGCGGAAAACCCGGTTGGCGAGAACAACGGTATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTAAAGGATAAATTTTGCCTGCGGGTCGTGCAGGGGCATGGCGCTCATCCAAGCGTATTGCGGGCAGCAGGGGCCGACCACGCCggcatgc                                                                                                                                                                                                                                 >  2:73641/1‑245 (MQ=255)                   aagctaaTCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTTATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGtt                                                                                                                                                                                                                            <  2:163681/249‑1 (MQ=255)                              aaaGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGtt                                                                                                                                                                                                               <  2:94029/251‑1 (MQ=255)                                        gaTAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAAc                                                                                                                                                                                                     <  1:413527/251‑1 (MQ=255)                                         aTAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCAGCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAAcc                                                                                                                                                                                                    <  1:11737/251‑1 (MQ=255)                                             cGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAata                                                                                                                                                                                                <  1:307097/251‑1 (MQ=255)                                             cGGGGCGCGACTGATGAAAATTATCATTCTGGGGGCCGGCCAGGTGGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAata                                                                                                                                                                                                <  2:40674/251‑1 (MQ=255)                                                  agctaCTGCTGAAAATTATCATTCTGGGTCCTGTCCATGTTGGCGGCACACTGACGGAAAACCTTGTTGTCGAGAACAACGATATTACTGTTTTCGATACCAACGGTGAGCGTCTTCTGACCTTACAGGATAAATTTGACCTGCGTGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGTGAAGCATGTGCCAACGACTCCGATATGCTGGTTTCTGTAACCAGTTCATATTAAACCAATATgtt                                                                                                                                                                                            <  2:209119/246‑3 (MQ=255)                                                     gACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTtgcc                                                                                                                                                                                        <  1:188929/251‑1 (MQ=255)                                                       cTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTtgcctg                                                                                                                                                                                      >  1:167460/1‑251 (MQ=255)                                                         gatgaAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTtgcctg                                                                                                                                                                                      <  1:15505/249‑1 (MQ=255)                                                                atgaTCATGCTGGGGGACGGCAAGGTTGGCTGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCATGGGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGATCAGATGAAACCAATATGGTTGCCTGCCAGGTAGc                                                                                                                                                                             <  2:191120/248‑1 (MQ=255)                                                                            ggTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTc                                                                                                                                                                                                              >  1:25715/1‑206 (MQ=255)                                                                            ggTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTc                                                                                                                                                                                                              <  2:25715/206‑1 (MQ=255)                                                                                 cGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCGGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAaca                                                                                                                                                            <  2:167460/251‑1 (MQ=255)                                                                                     cAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACAcccc                                                                                                                                                        >  1:295009/1‑251 (MQ=255)                                                                                           ggcggcACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAACTTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTCCTCACTTTTCAACACCCCTAatcg                                                                                                                                                  >  1:164087/1‑251 (MQ=255)                                                                                           ggcggcACACTGGCGGAAAACCTGGGTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAatcg                                                                                                                                                  >  2:61147/1‑251 (MQ=255)                                                                                                          gAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAAcc                                                                                                                                                                                                                   <  2:394988/171‑1 (MQ=255)                                                                                                          gAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAAcc                                                                                                                                                                                                                   >  1:394988/1‑171 (MQ=255)                                                                                                               ccTGGTTGGCGAGAACAACGATATTCCTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGGGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCCGATGAAACCAATATGGTTGCCTGCAAGGTAGCCTACTCACTTTTCAACACCCCTAATCGAATCGCTCGTATCCGCTCAc                                                                                                                              >  2:340795/1‑251 (MQ=255)                                                                                                                              caacGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGAGATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGcgcg                                                                                                               <  2:353199/251‑1 (MQ=255)                                                                                                                                     attgcTGTTGTCGATACCAACGGTGAGCTTCTGCGGACCTTACAGGAGAAATTCGACCTGCGGGTCGTGCAGGTGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTGGCCTGCCAGGTAGCTTACTCACTTTTCAACACCCCGAAGCGCATCGCGCGTATCCGCTCACCAGACTACGTGCGCGATGCCg                                                                                                         <  2:295127/246‑1 (MQ=255)                                                                                                                                           gttgtCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACTCCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGc                                                                                                   <  1:245545/250‑1 (MQ=255)                                                                                                                                                          ggagaaagTCTGCGGACCTTACAGGCTAAATTTGACCTGCGGGTCGTGAAAGGGCATGGCTCACAACCACGCGTATTCCGAGCAGCAGGTGCCGAGGACGCCGATATGCTCGGAGCTGTAAACAGATCCGATGAAAACAATAGGATAGCCTGCCAGGTAGCCTACTCACAATTCAACACCCCTAATCGCATCGCTCGTATCCGCCCACCAGACTACGTGCGCGATCACAATAAACTATttaattccgaggc                                                                                   >  2:443443/8‑245 (MQ=255)                                                                                                                                                           gggggCGTCTGCGGACCTTACAGGATAAATTTGCCCTGCGGGTCTTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCt                                                                                  <  2:187799/247‑1 (MQ=255)                                                                                                                                                                cGTCTGCGGACCTTACAGGATCAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACACCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCCGATGCtgtgac                                                                             >  2:467120/1‑249 (MQ=255)                                                                                                                                                                cGTCTGCGGACCTTACAGGATAAATTTGTCCGGCGGGGCGGGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGGAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCGCGTATCCGCCCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGcc                                                                             <  2:126078/251‑1 (MQ=255)                                                                                                                                                                     gCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGAtt                                                                         <  1:155528/250‑1 (MQ=255)                                                                                                                                                                          ccTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGatcat                                                                   >  1:58749/1‑251 (MQ=255)                                                                                                                                                                             tACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTg                                                                >  1:248185/1‑251 (MQ=255)                                                                                                                                                                               cAGGAGAAATTTGACCTGCGGTTCGTGCAGGGGCATGGCTCTCATCCCCTCGTATTGCGGTATGCAGGGGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATTGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCTCCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGAt                                                              <  2:238345/251‑1 (MQ=255)                                                                                                                                                                                               tGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGCTGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTg                                              >  1:142268/1‑251 (MQ=255)                                                                                                                                                                                                   ggTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGCTGCAGTGCCGATTGATCATCTGATCGCACCAGACCAGTGGGttt                                          >  2:138646/1‑250 (MQ=255)                                                                                                                                                                                                    gTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTAt                                         <  1:61147/251‑1 (MQ=255)                                                                                                                                                                                                    gTCGTGCAGGGGCAGGGCTCGCATCCACGCGTAGTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTTCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGAGTGATCATCTGATCGCACCAGAGCAGTGGGTTAt                                         <  2:295009/251‑1 (MQ=255)                                                                                                                                                                                                     gcGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATc                                        <  1:276684/250‑1 (MQ=255)                                                                                                                                                                                                      cGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTTCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTAGCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCg                                       <  2:124244/251‑1 (MQ=255)                                                                                                                                                                                                      cGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCg                                       <  1:118201/251‑1 (MQ=255)                                                                                                                                                                                                      cGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCg                                       >  2:181865/1‑251 (MQ=255)                                                                                                                                                                                                       gTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGCTGCCGATAAGCTATTTCATTCAGATGCTGTGCCGTTTGATCATCTGATCGCACCAGAGCAGTTGGTTCTCgt                                      >  2:309413/1‑250 (MQ=255)                                                                                                                                                                                                                      actcATCCACGCGTATTGCGGGAAGAAGGTTCCGACGACGCCGATCTGCTGGTTGCTGTAACCAGGTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCCTCGCTCGTATCCGCTCACCAGACTACGTGCACGATGCCGATAAGCTATTTCCTTCAGATGCTGTGCCGAATGATCATCTGATCGCACCATAGCTGTGGGTAATCTAGAATTTTTACCTACt                       >  2:328062/2‑251 (MQ=255)                                                                                                                                                                                                                       ctcATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTg                      >  1:156973/1‑251 (MQ=255)                                                                                                                                                                                                                                       gCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGgcg      <  1:206445/251‑1 (MQ=255)                                                                                                                                                                                                                                           gAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATTGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGGGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGGCGCAtt  <  2:131396/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         GACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGGCGCATT  >  NC_000913/3436461‑3436945

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