Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,384	A→C	31.3%	V289G (GTG→GGG)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,384	0	A	C	31.3%	102.3 / 22.1	77	V289G (GTG→GGG)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base A (41/11);  new base C (4/21);  total (45/32)
Fisher's exact test for biased strand distribution p-value = 2.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAA  >  NC_000913/3760159‑3760620                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               aTGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACCCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCGGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGAcgcc                                                                                                                                                                                                                                  <  1:66783/239‑1 (MQ=255)     ccTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGAcgcccgc                                                                                                                                                                                                                               >  1:1027718/1‑238 (MQ=255)         tgcgctgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTCGCTTTGCAGGCGGT‑CGGTCCTCTCTTTTTTACCAACCATTATTACTGCCTGCCGCCCGCCTTACCACTT‑TTCGCCGGCAAAAACCACCAG‑ACTGACGCCCCCCCCa                                                                                                                                                                                                                           >  2:1084389/1‑238 (MQ=255)              tgcgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAAAAATAATTCCCGCCTGGGGCCCGCCTAACAACTT‑TTCGCCCGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCac                                                                                                                                                                                                                      >  2:941188/1‑237 (MQ=255)                cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTTCTGCCTGCGGCCCGCCTCACAACTT‑GTCGCCGGTGAAACTCACCAG‑ACTGACGCCCCCCGCaacctcat                                                                                                                                                                                                                    >  1:234477/1‑232 (MQ=255)                cgGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑AGATCATCTCTTTTTTACCAACAAAAATTTCTGCCTGCCGCCCCCCTAACAACTT‑GTCTCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCCATCAACt                                                                                                                                                                                                                    >  1:70131/1‑238 (MQ=255)                    tAAGCGGGCGCAGGGTCGTTAATTGTTCACTCAGGGCTTCAGGGTGTAAATAAAGACGAAACGTGCCTCCCGGCGCGGCGAGGGTCTTTTTACCCGCGCGTAATGCTCGCTCCAGCGGGTGGCTTGGCGGGCGGG‑CGATCATCTCTTTTTTCCCAACAATATTCCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCGGGAGAACCTCCCCAG‑ACTGACGCCCGCCGCAATCAACTCCTg                                                                                                                                                                                                                <  2:817594/238‑1 (MQ=255)                            cgcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGGCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCCCCGCAATCAACTCCTGCGGCAttt                                                                                                                                                                                                        >  1:960411/1‑237 (MQ=255)                            cgcAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCCATCAACTCCTGCGGCATTg                                                                                                                                                                                                        >  1:444207/1‑238 (MQ=255)                                 ggTCGGTAATTTTTCACTCAGAGCTTCAGGGCGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGTG‑CGATTATCTCTTTTTTACTAACAATAATTCCTGCCTCCGGCCCCCCTTACATCTT‑TTCGCCCGGGAAACTCCACAG‑CATCACCCCCCGCGCAACACACTCCTCCGGCATTGGgtcct                                                                                                                                                                                                  >  2:900750/1‑237 (MQ=255)                                  gTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGGTCATCTCTTTTTTTCCACAAATAATTCCTGCCTGCGGCCCGCCTTACAACTT‑CTTGCCCGCGCAACTCCCCAC‑GCTTGAGCCCCCCGCCATCCACTCCTCCGGCCTTGGCTCttt                                                                                                                                                                                                 >  2:1174002/1‑239 (MQ=255)                                     ggTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCCACAATAATTCCTTCCTGCGGCCCGCCCACCAACTT‑GTGCCCGGCGCAACACACCCG‑TCTGAGGCCCGCCGCAATCAACTCCTTCGGTAGTGGctttctct                                                                                                                                                                                               >  2:1057581/1‑232 (MQ=255)                                         tttttcCACCAAGGGCTTGGGGGTTTAAATAAGGCCGCAAGGCGGCGTCCAGCGCTGGGGGGTTCTTTTTTCCCGCGCGTAATCCACGCTTCAGCGGGTGGCTTT‑CCGGCGGGCCGCTCATCTTTTTTTTCCCAACAATAATTCCTGCCTGCGCCCCGCCTAACAACTT‑GCCGCCGGAGAACCTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCaaa                                                                                                                                                                                         <  2:305989/238‑1 (MQ=255)                                         aTTTTTCCCCCGGGGCTTCGGGGTTTAAAAAAAGCCGCAACGTGGCTTCCGGCGCCGCGGGGTTCATTTTCCCCGCGCTTAATGCCCGCTTCAGCGGGTGGCTTTGCGGGCGGG‑CGATCATCTTTTTTTTCCAAAAAATACTCCCTGCCTGCGGCCCGCCTAACACCTT‑GCCGCCGGAAAACCTCCCCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCaaa                                                                                                                                                                                         <  2:502082/239‑1 (MQ=255)                                             ttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGGTCATCTCTTTTTTTACAACAAAAATTCCTGTCTGCGGGCCGCCTAACACATT‑TTCTCCGGGCAAACTCCCCAG‑ACTGCCGCCCGCCGCCAACAAATCCCGCGGCATTGTCTCTTT‑CGGCAAAACg                                                                                                                                                                                      >  2:905426/1‑238 (MQ=255)                                             ttCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGCTCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGGCCGCCTAACCACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCGTTGGCTCTTT‑CGGCAAACCGt                                                                                                                                                                                     >  1:21135/1‑239 (MQ=255)                                              tCACTCAGAGCTTCAGGGTATAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCTTTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCGGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAACCTCCCCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTa                                                                                                                                                                                    <  1:764069/239‑1 (MQ=255)                                              acccccAGAGTTTAGGGGTGAAAAAAAAGACGCAACGGGGCTCCCAGCGCGGGGAGGGTCATTTTATCCGCGGGTAATGCACGCTTCAGCGGGGGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAAAAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCCCCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTa                                                                                                                                                                                    <  2:1029858/234‑1 (MQ=255)                                               cACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGGTTATCTCTTTTTTACCACAAATAATTTCTGCCTGGCGCACGGCTAACCACTT‑TTTCCCGGAGAAACCCACCCG‑GCTGAAGCCCGCCGCAATCAACTTCAGCGGCATCGGCTCTTTCCGCCAAAccaa                                                                                                                                                                                     >  2:327966/1‑236 (MQ=255)                                                           cAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCCTCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCCCTTT‑CGGCAAACCGTACTGGCTAAGATc                                                                                                                                                                        >  1:364016/1‑238 (MQ=255)                                                              ggTATAAAAAAACACGCAACGGTGCCTCCGGCGCGGGGGGGGATATTTTACCCGCGCGTAATGCACGCCGCAGCGGGTGGCTTTGCGGGGGGG‑CGATCATCTCTTTTTCACCAAAAATAATCCCTGCCTGCGCCCCGCCTAACAACTT‑GTCGCGGGAGAACCCCCCCAG‑CCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGAc                                                                                                                                                                     <  2:278515/238‑1 (MQ=255)                                                                tgtAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAATCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCCTT‑CGGCAAACCGTACTGGCTAAGATCGcccc                                                                                                                                                                   >  1:272990/1‑237 (MQ=255)                                                                 gtAAATAAGGCCGCATCGTGCCTTCCTGGGCCGCGGGGTTATTTTTACCCGCGCGTAATGCGCGCTGCAGCGGGTGGTTTTGCGGGCGGG‑CTTTCATCTTTTTTTTACAAACAATACTCCCCCCCTGCGCCCCCCCAAACAACTT‑GTCGCCGGAGAACCCCCCCAG‑CCGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAg                                                                                                                                                                  <  2:1027718/238‑1 (MQ=255)                                                                   aaataaaGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGGTCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCCGAGAAACTCACCAG‑ACTGACCCCCGCCCCCATCCACTCCCGCGGCATTGGCTCTTT‑TCGCAAACCGCACTGGCTAAGATccaccaacca                                                                                                                                                               >  2:1098570/1‑235 (MQ=255)                                                                            cGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCa                                                                                                                                                                                                                           >  1:425278/1‑171 (MQ=255)                                                                              cAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCAGCGAGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑AGATCATCTCTTTTTTACCCACAATAATTCCTGCCTGCGGACCGCCTAACAACTT‑GTCTCCGGAGAAACTCACAAG‑ACTGACGCCCGCCGCCATCAACTCCTGCGGCCTTTGCGCTTT‑AGGCATACCACTATGGCTAACAACGACCAGCGAGCCCATGCCtt                                                                                                                                                    >  2:160737/1‑238 (MQ=255)                                                                                                               ctCCGCGTGTAACTCTCGCTTCGGGGGGTGCCTTGCCGGGCGGC‑CACTCACCTCTTTTGTACCAACAAGAATCCCCGCCCGCGCCCCCCACAACACCTT‑CTCGCCGGAAAACCCCCCCAG‑ACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGc                                                                                                                   <  2:1150263/238‑1 (MQ=255)                                                                                                               aTCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCCCTACGGGAACATCCAGCTCTTTg                                                                                                                    >  1:333412/1‑238 (MQ=255)                                                                                                                    cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCccgccg                                                                                                                                                                                                                             >  2:883705/1‑129 (MQ=255)                                                                                                                    cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCCGAGAAACTTACCAC‑ACTGACGCCCGCCGCAAACAACTCCTTCGGCATTTGCTCTTT‑CGGCAAACCCTACTGGCTCAAATCCACCAGCGAGCCACCGCCTAAATCCATCACTAAGGGCACATCCAGCTCTTTGCCgcg                                                                                                               >  2:154336/1‑236 (MQ=255)                                                                                                                    cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GCCGCCGGAGAAACTCACCAG‑ACTGACGCccgccg                                                                                                                                                                                                                             <  1:883705/129‑1 (MQ=255)                                                                                                                    cgcgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCCCTTTTTTTCCAACAATAATTCCTGCCTTCGGCCCGCCTAACCACTT‑GTCCCCGGAGAAATTCCCCAG‑ACTGACGCCCCCCCCCATCAACCCCTCCGGCAACGGCCCTTT‑CGGCAAAACCCACTTGCTAAGATTGACCCGCGTGCGCCTCCCTACAACAGGCCCCACGGGGACATCCCGCCCTTTGCCCAg                                                                                                               >  2:644311/1‑238 (MQ=255)                                                                                                                      cgTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGGTCATCTCTTTTTTACCAACAATAATTCCTTCCTGCGGCCCGCCTAACAACTT‑TTTGCCGGAGAAAATCACCAG‑CCTGACGCCCGCCCCCATCAACTCCTTCGGCCATGGCCCTTT‑CGGCAAACCGGACTGGCTAACATCGACCAGCGAGCCACTTCCCAAATCACTCCACACGGGGACATCCACCTCCTTGccccccc                                                                                                             >  2:410863/1‑233 (MQ=255)                                                                                                                                  ttCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAAt                                                                                                                                                 >  1:233750/1‑190 (MQ=255)                                                                                                                                   acAGCGGGTGGTTGTGCAGGCGGG‑CGATCATCTCTTTTTAACAAAAAAAAATCCCTGCCGGCGGCCCGCCTAACAACTT‑GCCGCCGGAGAAACTCCCCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCt                                                                                               <  1:405956/238‑1 (MQ=255)                                                                                                                                       cGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑TTCGCCGGAGAAACTCACCAG‑ACTGACGCCCCCCGCAATCAACTCCTGCGGCCTTGGCGCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACCGCCCTAATAAGTCCCTACGGGGACATCCCGCTCTTTTCCGACCGCCCCCAGTTCCGCtttc                                                                                            >  2:1174111/1‑236 (MQ=255)                                                                                                                                          gTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAAATCACCAG‑ACTGACGCCCGCCGAAATCAACTCCTGCGGCATTGGCTCTTT‑CCGCAAACCGTACTGGCTAAAATCCACCAGCGAGCCACTGCCCTAATCAGTCACTCCGGGAACATCCAGCTCCTTGCCCGCCGCCCCCCGTTCCCCTTTATCTa                                                                                        >  2:788547/1‑239 (MQ=255)                                                                                                                                            ggCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCCGCGAGCCACTGCCTAAATCAGTTACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCCGTTCCCCTTCATTTAt                                                                                       >  1:681652/1‑238 (MQ=255)                                                                                                                                            ggCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCCCCAGTTCCGCTTCATCTAt                                                                                       >  1:667153/1‑238 (MQ=255)                                                                                                                                            ggCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATc                                                                                      >  1:974549/1‑239 (MQ=255)                                                                                                                                                   cAGGCGGG‑CGATCATCTCTTTTTTCCCAACAATAATTCCTGCCCGCGCCCCGCCTAACAACTT‑GCCGCCGGAAAAACTCCCCAG‑ACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAg                                                                                                                                                                           <  2:131538/147‑1 (MQ=255)                                                                                                                                                      gcggg‑cgATCATCTTTTTTTACCCAACAATAATTCCTGCCTGCGGCCCGCCTACCAACTT‑GTCGCGGGAGAACCTCCCCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                             <  1:1174111/238‑1 (MQ=255)                                                                                                                                                      ccggg‑cgATCTTCTTTTTTTTACCACCAATAATTCCTCCCTGCGCCCCCCCTACCACCTT‑GTCGCCGGAGAACCTCCCCAG‑CCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTg                                                                             <  1:674227/237‑1 (MQ=255)                                                                                                                                                        ggg‑cgATCATCTCTTTTTTACAAACAATAATCCCTGCCTGCGGCCCGCCTAACAACTT‑GCCGGCGGAGAAACTCCCCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                           <  1:941188/238‑1 (MQ=255)                                                                                                                                                            cgATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGGACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAAAATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAAccc                                                                        >  1:202393/1‑238 (MQ=255)                                                                                                                                                              ctcatcTCTTTTTTACAAACAATAATCCCTGCCCGCGCCCCGCCTAACAACTT‑GTCGCCGGAGAACCTCACCAG‑ACGGACCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCt                                                                      <  1:644311/237‑1 (MQ=255)                                                                                                                                                              atcatcTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGc                                                                                                                                                                 >  1:131538/1‑147 (MQ=255)                                                                                                                                                                  tctcTTTTTTACCAACAAAAATCCCTGCCTGCGGCCCGTCTAACAATTT‑GCCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATg                                                                 <  1:410863/239‑1 (MQ=255)                                                                                                                                                                    tgttttttACCAACAATAGTTCCTGCGTGCGCCCCGCTAAACAACTT‑GCCGCCGGAGAACCCCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                               <  1:1084389/237‑1 (MQ=255)                                                                                                                                                                    tctTTTTTACCAACAATAATTCCCGCCTGCGCCCCGCCTAACACCTG‑GTCGCCGGAAAAACTCCCCAG‑ACGGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                               <  1:152371/239‑1 (MQ=255)                                                                                                                                                                    ctttttttACAAACATTACTCCCTCCCCGCGCCCCGCCTAACAACTT‑GCCGCCGGAAAACCCCCCCAG‑ACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                               <  1:444212/237‑1 (MQ=255)                                                                                                                                                                    cctgtagtACAAACAATACTCCCACCCTGCGGCCCGCCT‑ACAACCT‑CGGCCCGGGAAAACCCCCCAGTACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGCCCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCt                                                               <  2:272990/232‑1 (MQ=255)                                                                                                                                                                     ctTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTCAATCAGCCACTACGGTAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAAt                                                                  >  2:1083922/1‑235 (MQ=255)                                                                                                                                                                     ctTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAAt                                                                  >  2:1009469/1‑235 (MQ=255)                                                                                                                                                                         tttACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAACGCTGTAGt                                                          >  1:378630/1‑239 (MQ=255)                                                                                                                                                                           tcccAACAATATTCCCCCCCTGCGGCCCGCCAAACAACTT‑GGCGCGGAAGAATCTCCCCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTa                                                        <  1:469717/237‑1 (MQ=255)                                                                                                                                                                                acaaTAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAACATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCACCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGTTGAACCCCTGACTGCTGTAGTTACTGGt                                                   >  2:637232/1‑239 (MQ=255)                                                                                                                                                                                  aaaaatCCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGCCCCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTAt                                                 <  1:245273/236‑1 (MQ=255)                                                                                                                                                                                     actccctgcctcCGCCCCCCCAAACACCTG‑TTCCCCGGAGAAACCCCCCAG‑CTGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                              <  1:618753/235‑1 (MQ=255)                                                                                                                                                                                     acccccgccctgc‑GCCCCCCCAACCACTTGGTCCCCGGAAAATCTCCCCAC‑GCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                              <  2:97368/231‑1 (MQ=255)                                                                                                                                                                                     aatCCCTGCCGGCGCCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACGGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                              <  1:1174002/239‑1 (MQ=255)                                                                                                                                                                                     aatCCCCGCCTGCGGCCCGCCTAACAACTT‑GCCGCGGGAAAACCCCCCCAG‑ACTACCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTa                                              <  1:154336/239‑1 (MQ=255)                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                           >  2:987507/1‑203 (MQ=255)                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGaa                                                                           <  1:987507/203‑1 (MQ=255)                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                        >  2:290151/1‑238 (MQ=255)                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                        >  1:1016477/1‑238 (MQ=255)                                                                                                                                                                                            gcctgcGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTtca                                        >  1:952470/1‑238 (MQ=255)                                                                                                                                                                                               tgcGCCCCCCCAAACAACTT‑GCCGCGGGACAACCTCCCCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCaa                                    <  1:30437/239‑1 (MQ=255)                                                                                                                                                                                                     ccGCCTAAAAACTT‑GTCCCGGGAAAAACGCACAGC‑CCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTC‑TT‑CGGCAAACCGTACTGGCTAAGACCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCg                             <  1:584192/239‑1 (MQ=255)                                                                                                                                                                                                                    gTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCGAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGAATGTACTTTCATCaa                                    >  1:605714/1‑219 (MQ=255)                                                                                                                                                                                                                    gTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc                 >  2:60411/1‑238 (MQ=255)                                                                                                                                                                                                                    gTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc                 >  1:838445/1‑238 (MQ=255)                                                                                                                                                                                                                    gTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc                 >  1:88549/1‑238 (MQ=255)                                                                                                                                                                                                                    gTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTttcattc                 >  2:1075598/1‑238 (MQ=255)                                                                                                                                                                                                                         cGGAGAACCTCACCAG‑CCTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGcc           <  1:1098570/239‑1 (MQ=255)                                                                                                                                                                                                                             gaAACTCCCCAG‑ACTGCCGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGa        <  2:944497/238‑1 (MQ=255)                                                                                                                                                                                                                                   cACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATaa  >  2:1036102/1‑238 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               ATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTTCAGGGTGTAAATAAAGACGCAACGTGGCTTCCAGCGCCGCGAGGGTCATTTTATCCGCGCGTAATGCACGCTTCAGCGGGTGGCTTTGCAGGCGGG‑CGATCATCTCTTTTTTACCAACAATAATTCCTGCCTGCGGCCCGCCTAACAACTT‑GTCGCCGGAGAAACTCACCAG‑ACTGACGCCCGCCGCAATCAACTCCTGCGGCATTGGCTCTTT‑CGGCAAACCGTACTGGCTAAGATCGACCAGCGAGCCACTGCCTAAATCAGTCACTACGGGAACATCCAGCTCTTTGCCGAGCGCCACCAGTTCCGCTTCATCTATCGCTTTGGTGAACCCCTGAATGCTGTAGTTACTGGTATGTACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAA  >  NC_000913/3760159‑3760620

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