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	98.0%	154.3 / NA	49	E60E (GAG→GAA)	trkA	NAD‑binding component of Trk potassium transporters
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (27/21);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTG  >  NC_000913/3436452‑3436929                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         gAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAgg                                                                                                                                                                                                                                     >  1:325777/1‑251 (MQ=255)  aaGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGt                                                                                                                                                                                                                                    >  2:315034/1‑251 (MQ=255)                 tcttttacaatatgatcaccaacaaGTGATGAGCTACCCGGTCGCGACTGATGAAAATCATCATTCTGGGTGCCGGCCAGGATGGCCGCACACTGGAGGAAAACCTGGTTGGCGAGCACAGCCATATTACTGTTGTCAATACCAACGGAGAGCGTCTCCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGAAGGGGCCTGGCTCTCCTCCACTCGTCTTTGGGGGAG‑CAGGGTCCGACGACGCCGat                                                                                                                                                                                                                     >  2:433781/24‑251 (MQ=255)                    tttACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATg                                                                                                                                                                                                                  >  1:102150/1‑251 (MQ=255)                             agctaaTCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGTCCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCt                                                                                                                                                                                                         >  1:15000/1‑251 (MQ=255)                                taaTCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAACATTATCATTATGGGTGCCGGACAGGTTGGCGGCACACTGGCGGCAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCATACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGAAGCCGATATGCTGGTTGCTgg                                                                                                                                                                                                       >  1:435098/1‑249 (MQ=255)                                    cAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTAGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACCGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCCGGTGCCGACGACGCCGATATGCTGGTTGCTGTAAcc                                                                                                                                                                                                   >  1:246452/1‑250 (MQ=255)                                     aaaaaGTGATGAGATAACGGGTCGCGACTGATGAAAATTATTCTTCTGGGTGCCGGCCAGGTTGGCGCCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGCTATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACATTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCCAGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAAcc                                                                                                                                                                                                   >  1:183395/1‑249 (MQ=255)                                            gatgaGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGAt                                                                                                                                                                                          <  1:69612/251‑1 (MQ=255)                                                  aTAACGGGTCGCGACTGATTAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACAATGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGAGATGCTGGTTGCTGTAACCAGTTCAGATGAAAcc                                                                                                                                                                                    <  2:85843/251‑1 (MQ=255)                                                    aaCGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCa                                                                                                                                                                                   >  1:160398/1‑250 (MQ=255)                                                         gTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATgg                                                                                                                                                                             >  2:146005/1‑251 (MQ=255)                                                           cgcgACTGATGAAAATTATCATTCTGTGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGt                                                                                                                                                                            <  2:348921/250‑1 (MQ=255)                                                            gcgACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACGGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCATGGGCATGGCGCTCATCCGCGCGT‑AGTGCTGGAAGCAGGTGCCGACGACGCCGATATGGTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTtg                                                                                                                                                                          <  1:274528/251‑1 (MQ=255)                                                              gACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTtgcc                                                                                                                                                                        <  2:174593/251‑1 (MQ=255)                                                                  gatgaAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTtgcctg                                                                                                                                                                      <  2:325777/249‑1 (MQ=255)                                                                   atgaAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCa                                                                                                                                                                   <  2:344498/251‑1 (MQ=255)                                                                            aTCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCt                                                                                                                                                           <  1:363484/250‑1 (MQ=255)                                                                                               aGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCt                                                                                                                                       >  1:92152/1‑251 (MQ=255)                                                                                                 gTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGCACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACCGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTa                                                                                                                                      >  1:167856/1‑250 (MQ=255)                                                                                                                gCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTAt                                                                                                                      >  1:217728/1‑251 (MQ=255)                                                                                                                                  gagaACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGCCTACGt                                                                                                    >  1:153919/1‑251 (MQ=255)                                                                                                                                    gaACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTgc                                                                                                  <  1:106888/251‑1 (MQ=255)                                                                                                                                        aacGATATTACTGTTGTCGATACCAACGGTGAGCGTATGCGGCCCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGcgcg                                                                                               >  2:285404/1‑250 (MQ=255)                                                                                                                                              atTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCg                                                                                         <  1:146005/250‑1 (MQ=255)                                                                                                                                                tACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGAt                                                                                       >  1:1891/1‑250 (MQ=255)                                                                                                                                                   tgttgtCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAACTTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGc                                                                                   >  2:189060/1‑251 (MQ=255)                                                                                                                                                       gtCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTAt                                                                                <  2:256150/250‑1 (MQ=255)                                                                                                                                                         cGATACCAACGGTGAGCGTCTGCGGACCTTACAGGCTAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATttc                                                                             >  2:73864/1‑251 (MQ=255)                                                                                                                                                           aTACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATttcat                                                                           >  2:213789/1‑251 (MQ=255)                                                                                                                                                             aacAACGGGGAGGGGCTGCGGACCGTACAGGATAAAATTGACCTTCGGGGCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACCCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATttcattc                                                                         <  2:160398/249‑1 (MQ=255)                                                                                                                                                              ccAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAGCCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTTATCGCATCGCTCGTATCCGCTCAGCAGACTACGTGCGCGATGCCGATAAGCTATttcattca                                                                        <  1:196457/251‑1 (MQ=255)                                                                                                                                                                   ggTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGAAGGGGCATGGCTCTCATACACGCGT‑ATTGCGGGAAGCAGGTGACGACGACGCCGATATGATGGTTGCTGTAACCAGTTAAGATGAAACCAATATGGTTGCCTGACAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGc                                                                   >  1:53760/1‑251 (MQ=255)                                                                                                                                                                    gTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCt                                                                  <  2:102150/251‑1 (MQ=255)                                                                                                                                                                         ggTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCTCCAGACTACGTGCGCGATGGCGATAAGCTAtgttaggcaggtgctgtgcc                                                             <  2:206500/250‑20 (MQ=255)                                                                                                                                                                               cGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGCTATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTg                                                        >  1:290284/1‑250 (MQ=255)                                                                                                                                                                                ggACCTTACAGGATAAATTTGACATGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGACGACGACGCCGATATGCTGGTTGATGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCC‑TAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGatc                                                     >  1:407845/1‑251 (MQ=255)                                                                                                                                                                                      tACAGGATAAATTTGACCTGCGGGTCGTGCAGGAGCATGGCTCTCATCCACGCGT‑ATTTCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTGGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGatcatc                                                  <  2:260229/249‑1 (MQ=255)                                                                                                                                                                                       acaggatgagtttcacgagCGGGTGGTGCAGGCGCATGGCTCTCATCCAAGCGT‑ATTGCGCGAAGCAGGCGCCGACGACACCGCCATGCTGGTGGCTGTTGCCAGTTCGGATGCAACCAACATGGGTGCCTGCCAGGTAGCATGCTCACTTTTCAACACCGCTAATCGGCTCGCTGGTATACGCTCCCCAGACTAGGTCGGAGACGGCGCTAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTg                                                >  2:264946/19‑250 (MQ=255)                                                                                                                                                                                                  ttGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGATGCTGTAACCAGTTCAGATGCAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGATCGTATCAGCTCACCAGACTACGTGCGCGATGCCGACACGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCag                                      >  1:75141/1‑249 (MQ=255)                                                                                                                                                                                                    gACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTGCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGc                                   <  1:73864/250‑1 (MQ=255)                                                                                                                                                                                                     aCCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCAGCGCTCGTAGCCGGTCACCAGACTACGTTCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAg                                 <  1:435437/251‑1 (MQ=255)                                                                                                                                                                                                         gCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTg                              <  1:213789/250‑1 (MQ=255)                                                                                                                                                                                                           gggTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGtt                           >  1:22950/1‑251 (MQ=255)                                                                                                                                                                                                                gTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCg                       >  1:247097/1‑250 (MQ=255)                                                                                                                                                                                                                  gCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGat                     >  2:33396/1‑250 (MQ=255)                                                                                                                                                                                                                                  cATCCACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTGGCGGTATCCAGGTCAGATGAATGCAATATGGTTGCCTGCCAGGGAGCCTACTCGCTTTTCAACACGCCTAATCGCATCGCTCGTAGCCGCTCACCAGACGACGTGCGCGATGCCGGTAAGCTAGTTCATTCAGATGCTGTGCCGATTGATCATCGTAGCGCTCGAGATCAGTTGGTTTTCGATAAGATTTACCGACTg      <  2:286961/249‑1 (MQ=255)                                                                                                                                                                                                                                     ccACGCGT‑ATTGCGGGAAGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTg  <  2:217728/250‑1 (MQ=255)                                                                                                                                                                                                                                      cACGCGT‑ATTGCGGGAAGCAGGCGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATCTTTACAGACTGATTg  >  2:129510/1‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         GAAGAGGGCGACGGCTTCTTTTACGCTAAGCTAATCAAAAAGTGATGAGATAACGGGTCGCGACTGATGAAAATTATCATTCTGGGTGCCGGCCAGGTTGGCGGCACACTGGCGGAAAACCTGGTTGGCGAGAACAACGATATTACTGTTGTCGATACCAACGGTGAGCGTCTGCGGACCTTACAGGATAAATTTGACCTGCGGGTCGTGCAGGGGCATGGCTCTCATCCACGCGT‑ATTGCGGGAGGCAGGTGCCGACGACGCCGATATGCTGGTTGCTGTAACCAGTTCAGATGAAACCAATATGGTTGCCTGCCAGGTAGCCTACTCACTTTTCAACACCCCTAATCGCATCGCTCGTATCCGCTCACCAGACTACGTGCGCGATGCCGATAAGCTATTTCATTCAGATGCTGTGCCGATTGATCATCTGATCGCACCAGAGCAGTTGGTTATCGATAATATTTACCGACTGATTG  >  NC_000913/3436452‑3436929

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: ‑