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	100.0%	135.9 / NA	40	E60E (GAG→GAA)	trkA	NAD‑binding component of Trk potassium transporters
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/23);  total (17/23)

CCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGGCGCA  >  NC_000913/3436450‑3436943                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       ccGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTGTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCa                                                                                                                                                                                                                                                     <  1:291232/251‑1 (MQ=255)           cGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCTTATTGCGGGAAGCAGGTGCcgacg                                                                                                                                                                                                                                           >  1:443641/1‑251 (MQ=255)                       ttACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATg                                                                                                                                                                                                                                <  1:342332/250‑1 (MQ=255)                                       aaaaaGGGATGAGATAACGGGGCGCGACTGATGAAAATTATCATTATGGGTGCCGGCCAGGTTGGCGGACCACTGGCGGAAAACCTGGTTGGCGAGCACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCCCGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAAcc                                                                                                                                                                                                                 >  2:319846/1‑249 (MQ=255)                                                    aTAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCAGGAAGCAGGGGCCGACGACGCCGATATTCTGGTTGCTGTAACCCGTTCAGATGAAAcc                                                                                                                                                                                                  >  2:238568/1‑251 (MQ=255)                                                              gcgACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGTCCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGAGATGCTGGTTGTTGTAACCAGTTCAGATGAAACCAATATGGTtg                                                                                                                                                                                        <  1:337667/251‑1 (MQ=255)                                                                                 aTTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCCTCCATGCGTATTGCGGGAAGCAGGTGCCGCCGACGCCGATCTGCTGGTTGCTGTAACCAGTTCAGATGCAACCCATATGGTTGCCTGCCAGGTCGCCTAc                                                                                                                                                                       >  2:82093/1‑249 (MQ=255)                                                                                 aTTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCATGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTAc                                                                                                                                                                       <  1:82093/249‑1 (MQ=255)                                                                                       ggTGCCGGCCAGGGTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACtttt                                                                                                                                                               <  2:186109/251‑1 (MQ=255)                                                                                             ggCCAGGTTGGCGGCCCACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCCGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGTTAGCCTCCTCACTTGTCAacac                                                                                                                                                         <  2:391144/251‑1 (MQ=255)                                                                                              gCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACAcc                                                                                                                                                        >  1:388989/1‑251 (MQ=255)                                                                                                      ggcggcACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAatcg                                                                                                                                                <  1:381561/251‑1 (MQ=255)                                                                                                        cggcACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAatcgc                                                                                                                                               >  1:19775/1‑250 (MQ=255)                                                                                                              ccTGGCGGAAAGCCGGGTGGGCGAGAACAACGATATTACGGTTGGCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGTCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGGAGCCTACTCTCTTTTCAACACCCCTAATCGCATCGCTc                                                                                                                                        <  2:47923/250‑1 (MQ=255)                                                                                                                                           acGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGAt                                                                                                           <  1:59850/251‑1 (MQ=255)                                                                                                                                            cGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATg                                                                                                          >  2:298561/1‑251 (MQ=255)                                                                                                                                              atatTACTGTTGTCGATACCAACGGTGAGCGTCTGCCGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTCCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGcc                                                                                                        <  2:347247/251‑1 (MQ=255)                                                                                                                                               tatTACTGTTGTCGATACCACCGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACCCCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCg                                                                                                       >  1:171887/1‑251 (MQ=255)                                                                                                                                                      gttgtCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCt                                                                                                >  2:306721/1‑251 (MQ=255)                                                                                                                                                           cGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGGCGACGCCGATATGCTGGTTGCTGGAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAatcgcatcg                                                                                                                                           <  1:189211/203‑1 (MQ=255)                                                                                                                                                            gATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCAAGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGAGGCCGATACACTATttc                                                                                           >  1:280291/1‑250 (MQ=255)                                                                                                                                                               aCCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATttcatt                                                                                        <  2:187271/250‑1 (MQ=255)                                                                                                                                                                   aCGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATttcattc                                                                                       >  2:212507/1‑247 (MQ=255)                                                                                                                                                                     ggTGAGCGTCTGCGGACCTTACAGGCTAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGc                                                                                 >  1:354894/1‑251 (MQ=255)                                                                                                                                                                                gCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGAtt                                                                       <  1:306721/250‑1 (MQ=255)                                                                                                                                                                                  ggACCTTACAGGATAAATGTGACCTGCGTGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGTGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTTCCAGTTAGCCTACTCACTGTTCAACACCCCCAAGCGCATCTCTCGTATCCGCTCACCAGACCACGTGCACGATGCCGATAAGCTATTTCATTCATATGCTGTGCCTATTGat                                                                    <  2:403788/251‑1 (MQ=255)                                                                                                                                                                                    aCCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATg                                                                                                                                                                                                                                <  1:12081/93‑1 (MQ=255)                                                                                                                                                                                    aCCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATg                                                                                                                                                                                                                                >  2:12081/1‑93 (MQ=255)                                                                                                                                                                                         aCAGGAGAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCGCGCGTATTGCGGGAAGCAGGTGCCGTCGTCGACGATATGCTGGTTGCGGTAACCAGGTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCGACTCTCTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGAGGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCtta                                                             <  2:201483/251‑3 (MQ=255)                                                                                                                                                                                           aGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATc                                                           <  1:37176/251‑1 (MQ=255)                                                                                                                                                                                           aGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATc                                                           <  1:298561/251‑1 (MQ=255)                                                                                                                                                                                                aaaTTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCAcc                                                      >  1:185921/1‑251 (MQ=255)                                                                                                                                                                                                  aTTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCTCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCag                                                    <  1:198579/251‑1 (MQ=255)                                                                                                                                                                                                                    gCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTAATCGCACCAGAGCAGTTGGTTATCGat                                   <  1:76463/250‑1 (MQ=255)                                                                                                                                                                                                                           cagggCTCGCATCCACGCGTATTGCGGTAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCTCTTTTCAACACCCCTAACCGCATCGCTCGTATCCGCTCCCCAGACTACGTGCGCGATGCCGATAAGCTATATCATTTAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTa                           <  1:181794/248‑1 (MQ=255)                                                                                                                                                                                                                             tGGCTCTCATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTAcc                         >  2:155513/1‑251 (MQ=255)                                                                                                                                                                                                                                   tcATCCACGCGTATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATACTATTTCCCGACTg                    >  2:302041/1‑250 (MQ=255)                                                                                                                                                                                                                                        cccgcgTATTTCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTACCCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCATATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAg              <  2:283439/249‑1 (MQ=255)                                                                                                                                                                                                                                                    gggAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGgcgc   >  1:136801/1‑250 (MQ=255)                                                                                                                                                                                                                                                    gggAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGGTCAGATGAAACCAATATGGTTGCCTGTCAGGTAGCCTACTCACTTTTCAACACCCCTAATCTCATCTCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGAGCGCACCAGAGCAGTTGGTTATCGATAATAGTTACCGACTGATTGAGTATCCCTGCGCa  <  2:415771/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       CCGAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGTATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTGAGTATCCCGGCGCA  >  NC_000913/3436450‑3436943

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