Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,669,782	A→T	E64V (GAG→GTG)	treF →	cytoplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,669,782	0	A	T	93.6%	127.8 / NA	47	E64V (GAG→GTG)	treF	cytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (24/20);  minor base G (0/3);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGACGAACTGATGATCGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGAGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGCTA  >  NC_000913/3669545‑3670017                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            tAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGACGAACTGATGATCGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTc                                                                                                                                                                                                                                 <  1:7762/250‑1 (MQ=255)            cGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGACGAACTGCTGATCGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAAtttt                                                                                                                                                                                                                     >  1:403876/1‑251 (MQ=255)             gATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGACGAACTGATGATCGAAGTCGATCTCTGCTATGAGCTGGGCCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTc                                                                                                                                                                                                                    <  1:106012/251‑1 (MQ=255)                 aTTAGCGACGAATTTCGGCGGTCGGATCCTTATGCTAAATCAGAAAATTCCAAACCCTAATCCAGACGAACTGATGATCGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGCGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCATGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCg                                                                                                                                                                                                                >  2:222001/1‑251 (MQ=255)                 aTTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGACGAACTGATGATCGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCg                                                                                                                                                                                                                <  1:95326/251‑1 (MQ=255)                           aaTTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGACGAACTGATGATCGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAAc                                                                                                                                                                                                      >  1:14496/1‑251 (MQ=255)                                  gAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGACGAACTGATGATCGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTccc                                                                                                                                                                                               <  2:421600/251‑1 (MQ=255)                                                                 tCAAAACCCTAATCCAGACGAACTGATGATCGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGat                                                                                                                                                                 >  1:169967/1‑250 (MQ=255)                                                                   aaaaCCCTAATCCAGACGAACTGATGATCGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATc                                                                                                                                                              <  1:301171/251‑1 (MQ=255)                                                                         cTAATCCAGACGAACTGATGATCGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATc                                                                                                                                                        >  1:370633/1‑251 (MQ=255)                                                                                               cGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGt                                                                                                                                   <  1:379298/250‑1 (MQ=255)                                                                                                        tctcTGCTATGAGCTGGACCCGTATGAATTCAAACTGGATGAGATGCTCGCGGCAGAACCGGAACCCGAGATGATTGCAGGGCTGCCTGCCTCTGATGCGCTGACCCCTGCCGATCGCTATCTCGACCTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGAGCACCTAAAATGAACCCGCTGGATCTCTTAATCCGCTCCCGTAAAGTGCGCCGTCATCGTGAtt                                                                                                                         >  2:78621/1‑251 (MQ=255)                                                                                                             gctaggggCTGGACCCGAATAAATTAAAACTGGATGAGATGTTCAAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGa                                                                                                                     <  2:397692/243‑1 (MQ=255)                                                                                                             gCTATGAGCTGGCCCCGTATGAATTAAAACTGGATGAGATGATCAAAGCAGACCCGGAACCCGAGAGGATTGAAGGGCGGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTTTCCTTGCATGTCCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCTTCTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGAc                                                                                                                    <  2:463441/251‑1 (MQ=255)                                                                                                               tATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGCGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAgaaacgccgtgctcgagcaaaatttc                                                                                                                  >  2:41089/1‑226 (MQ=255)                                                                                                               tATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACtt                                                                                                                  >  1:144161/1‑251 (MQ=255)                                                                                                                       ggACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAATTGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTTCGCCGTCATCGTGATTTTGACTTGCGCAAGt                                                                                                          >  2:109951/1‑251 (MQ=255)                                                                                                                                tGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGaa                                                                                                  >  2:119350/1‑250 (MQ=255)                                                                                                                                   aTTAAAACTGGATGAGATGATCGAGGCACACCCGGAACCCGGGATGATTGAAGGGCTGCCTGCTTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAAcc                                                                                              <  2:457574/251‑1 (MQ=255)                                                                                                                                     tAAAACTGGCTGAGATGATCGAGGCAGAACCAGAACCCGAGATGATTGAAGGGATGCATGCCTCTGATGCGCTGACGCCTGCCGCTCGCTATCTAGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGAAATCTTAATCCGCTACCGTAAAGTGAGCCGGCCTCGTGCTTTTGACTTGCGCACGTTTGTTGACAACCAc                                                                                            >  2:180258/1‑251 (MQ=255)                                                                                                                                       aaaCTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACtt                                                                                          >  2:149742/1‑251 (MQ=255)                                                                                                                                           tGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACCGTAAAACCTTTCCCGACTGCGCACCTAAACTGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGCTTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTgg                                                                                      >  1:461651/1‑251 (MQ=255)                                                                                                                                              atgagatgaTCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACATCTGGCTg                                                                                   >  2:222803/1‑251 (MQ=255)                                                                                                                                                 agatgaTCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCg                                                                                >  2:211425/1‑251 (MQ=255)                                                                                                                                                  gatgatCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCACGACTGAGCACCTAAAATGGAACCGCTGGATATCTTAATACGCTACCGTAAAGTGCGCCGTAATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCgg                                                                               >  1:496937/1‑251 (MQ=255)                                                                                                                                                      atCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGAAGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGAGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACATAAAATGGACCCGCTGGATATCTTACTCCGCTCCCGTAAATGGCGCAGTCTTCGTGATTTTCACTTGCGCAAGTTTGTTGAAAACAACTTCTGGCTGCCgacggt                                                                           >  2:136266/1‑246 (MQ=255)                                                                                                                                                             cAGAACGGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTTGCCCGAACTGTTCGTGCATGTGCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTcc                                                                    <  2:181644/251‑1 (MQ=255)                                                                                                                                                              aaaaCCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTTATGCGCTGACGCCTGCCGATCGCTATCTCGAGCTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTATCGTAAAGTGCGCCGTCATCGTGATTTTTACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCa                                                                   <  1:328860/249‑1 (MQ=255)                                                                                                                                                              aGAACCGGAACCCGACATGATTTAAGGGCGGCCTGCCTTCGAGGCGCTGACGCCTGCCGATCGATATCTCGAACTATTCGTGCATGTACAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACGACGCACCTAAAATGGACCCTCTGTATATCTTAATCCGCGACCGTTAATTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCa                                                                   <  1:494889/251‑1 (MQ=255)                                                                                                                                                               gAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTcc                                                                    >  1:364233/1‑249 (MQ=255)                                                                                                                                                                aaCCGGAACCCGAGCTGATTGAAGGGCTGCCTGCCTCTTATGCGCGGCCGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGc                                                                 <  2:14496/251‑1 (MQ=255)                                                                                                                                                                    ggAACCCGGGAGGATGGAAGGGCTGCCTGCCGCTGTGGGGCTGACGCCGGCCGATCGCTATCTCGACCTTTTCGGGCAGGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCTGCATCGTGATTGTGACTTGCGCAATTTTGTTGAAAACCACTTCTGGCTGCCGGAGGGCTACTCCAGCGAgt                                                             <  1:52351/251‑1 (MQ=255)                                                                                                                                                                           gagaTGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATc                                                      <  2:370633/251‑1 (MQ=255)                                                                                                                                                                                gatTGAAGGGCAGCCTGCCTATGATGCGCTGACGCCTGCCGATCGATATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTCACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGTGGTCTCATCACGCGCGTATGTATCGGAcc                                                 >  2:498235/1‑251 (MQ=255)                                                                                                                                                                                            gcctgccttggttgggcTGAGGCCTGCCGATCGCTATCTCGATCTGTTCGTGCTCGTCCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCTCCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTTGCTGCCGGAGGTCTTCTTCAGCGATTATGTAGCGGACCCGCAAAATTcc                                      <  2:260390/235‑1 (MQ=255)                                                                                                                                                                                            gcctgcctCTGATGCGCTGGCGCCTGCGGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTccc                                     <  1:198941/251‑1 (MQ=255)                                                                                                                                                                                             cctgcctCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCt                                    >  1:289855/1‑251 (MQ=255)                                                                                                                                                                                                  ctctGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGaa                                 >  1:449188/1‑249 (MQ=255)                                                                                                                                                                                                         gctcTGACGCCTGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCata                         <  1:4116/247‑1 (MQ=255)                                                                                                                                                                                                         gagctaaCGCCTGCCGATCGCTATCGCGAAATGTTCGGGCATGTTCAGTGGGCAAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGTACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCatat                        <  2:375979/245‑1 (MQ=255)                                                                                                                                                                                                              taCGCCTGCCGATCGCTATCTCGAACTGTTCGTGCAGGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTGTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTTAAAGAGCATATCGAc                    <  1:259586/249‑1 (MQ=255)                                                                                                                                                                                                                    tGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTg                                                                                       >  1:487538/1‑177 (MQ=255)                                                                                                                                                                                                                    tGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTg                                                                                       <  2:487538/177‑1 (MQ=255)                                                                                                                                                                                                                    tGCCGATCGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCtgt             >  1:320032/1‑251 (MQ=255)                                                                                                                                                                                                                     gCCGATCGATATCTCGAACTGTTCGGGCATGTTCAGTCGGCGAAAATTTCCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGTTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCtgt             <  2:106901/250‑1 (MQ=255)                                                                                                                                                                                                                           cGCTATCTCGAACTGTTCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGACGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGt      >  1:211030/1‑251 (MQ=255)                                                                                                                                                                                                                               agctcGAACTGTGCGTGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTGATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGCTa  <  1:210238/249‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            TAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGACGAACTGATGATCGAAGTCGATCTCTGCTATGAGCTGGACCCGTATGAATTAAAACTGGATGAGATGATCGAGGCAGAACCGGAACCCGAGATGATTGAAGGGCTGCCTGCCTCTGATGCGCTGACGCCTGCCGATCGCTATCTCGAACTGTTCGAGCATGTTCAGTCGGCGAAAATTTTCCCCGACAGTAAAACCTTTCCCGACTGCGCACCTAAAATGGACCCGCTGGATATCTTAATCCGCTACCGTAAAGTGCGCCGTCATCGTGATTTTGACTTGCGCAAGTTTGTTGAAAACCACTTCTGGCTGCCGGAGGTCTACTCCAGCGAGTATGTATCGGACCCGCAAAATTCCCTGAAAGAGCATATCGACCAGCTGTGGCCGGTGCTA  >  NC_000913/3669545‑3670017

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