Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	23,161	A→T	P257P (CCA→CCT)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	23,161	0	A	T	83.3%	113.8 / 15.9	54	P257P (CCA→CCT)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (3/6);  new base T (23/22);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCGGCACGGTAAAAGGTG‑CGGAG‑CTTGAGCTGCTGCGC‑TTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACG  >  NC_000913/22946‑23405                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         cACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCCCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCTCTGGTGCAGATCGACGGTCa                                                                                                                                                                                                                       <  1:176412/251‑1 (MQ=255) cACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTc                                                                                                                                                                                                                        <  1:381459/250‑1 (MQ=255)  aCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAg                                                                                                                                                                                                                      <  1:242473/251‑1 (MQ=255)       tGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGt                                                                                                                                                                                                                 >  1:123931/1‑251 (MQ=255)        gcggtGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTg                                                                                                                                                                                                                <  1:28281/247‑1 (MQ=255)         cGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGACCGTGa                                                                                                                                                                                                               <  2:80197/251‑1 (MQ=255)           ttGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCGGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTCTTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCCGGCCGTGAtt                                                                                                                                                                                                             >  1:318516/1‑251 (MQ=255)           ttGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGAtt                                                                                                                                                                                                             >  1:283002/1‑251 (MQ=255)               ctgacGTTCTGCGCTGGCGGAAGCGGCAGTTGAGTATTACGACAAAACGTCTCCGTCCATCGACGGTGCTTTTCAGGCGGTCGATCCGGATGCACTGACAGCAACATTTGCCGTAAGCAACGTTAACGGCCCTATCTCGCTGGTGATCTGGACAACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTTTTGCACCATATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGcagctcttct                                                                                                                                                                                                           >  2:270170/5‑240 (MQ=255)                 gCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTAAGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCACGCAATCTCTATTTCACCTGATTTAGACTATGCGCTGGTGCAGAACGACGGGCAGGCCGTGGTTCTGGCg                                                                                                                                                                                                       >  1:193842/1‑251 (MQ=255)                         gcgcTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGACTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCACTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGa                                                                                                                                                                                                  <  1:78572/248‑1 (MQ=255)                           ccTGGAGGAAGCGGAATTTGAGTTTTACGACAAAACTTCTCCTGCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTTGCGAAAGAGCTgg                                                                                                                                                                                             <  2:300165/250‑1 (MQ=255)                                cggaagcggcagTTGAGTATTACGACAAAACTTCTCCGTCCATAGACGTTGCTTTTCAGGCGGTCGATCAGGATGCAATGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGAACCTGATTTCGACTATGCGCTGGTGAAGATCGACGTCCAGGCCGTGATTCTGGCGAAAGATCTGGTTGa                                                                                                                                                                                         >  2:362690/1‑250 (MQ=255)                                    agcggaagTTGAGTATTACGACAAAACTTCTCCGTGCATCTACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCCCGCCGTTGACTCTGCCTGCCAACCGCGCAATCTCTATTGCAGCAGATTTCGACTATGCGCTGGTGCAGAGCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCg                                                                                                                                                                                    <  2:81508/251‑1 (MQ=255)                                       ggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGTTTCTGGCGAAAGATCTGGTTGAAAGCGTa                                                                                                                                                                                  >  2:193650/1‑250 (MQ=255)                                          agTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGc                                                                                                                                                                              <  1:144319/251‑1 (MQ=255)                                           gTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCa                                                                                                                                                                             <  1:113558/251‑1 (MQ=255)                                             tGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGAGGTCGATCAGGCTGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGAGATCTGGACAACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATATCTATTGCAACTGATTTCGACTATGCGCTGGTGCAGATCGCCGGACAGCCCGTGATTCTGGCGAGAGATCTGGTTGGAAGCGGAGTGCAGc                                                                                                                                                                           >  2:850/1‑251 (MQ=255)                                                           aaaCTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAc‑c                                                                                                                                                            >  1:230837/1‑251 (MQ=255)                                                           aaaCTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAc‑c                                                                                                                                                            <  1:70087/251‑1 (MQ=255)                                                                    cGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACAcc                                                                                                                                                   <  1:73616/251‑1 (MQ=255)                                                                    cGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACAcc                                                                                                                                                   >  1:137541/1‑251 (MQ=255)                                                                               gTTGATTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGAAACGTTAACGGCCCACTCTCGCTGGGGATCTGGACCACTACGACGTGGCCTCTGCCTGCGAACCTCGCAATCTCTATTGCACCTTATTTCGACTATGCGCTGGCGCAGAGCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTTAAAGAGTAATGCAGCTTCTCGGCGTAGC‑CGATTaccacaatttcgacac                                                                                                                                        >  2:25393/1‑237 (MQ=255)                                                                                 tGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTg                                                                                                                                                                <  1:293802/226‑1 (MQ=255)                                                                                 tGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGagatcggaagagcgacgggtcggg                                                                                                                                        >  2:293802/1‑227 (MQ=255)                                                                                     tttcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTaa                                                                                                                                   <  1:211986/247‑1 (MQ=255)                                                                                       tcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTa                                                                                                                                                       >  1:207447/2‑228 (MQ=255)                                                                                       tcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCGCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTa                                                                                                                                                       <  2:207447/227‑1 (MQ=255)                                                                                         aGGCGGTCGATCAGGGTGCACTGAAAGCAAAATTTGCCGTACGCAACGTTAACCGCCCAAGCTCGCTGGTGCTCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTGTGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTaaaatgt                                                                                                                              <  1:346754/251‑4 (MQ=255)                                                                                             ggTCGCTCAGGATGCACTGAAAGCAAAATTTGCCGTACGCAACGTTAACGGCCCAATCTCGCTGTTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTACACCTGATTTCGACTCTGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCGTGCGACATATCTGTTTGAAGAAGGAATACAGAGGATCGGCGTAAC‑CGATTCCACCATtctcaacacgtaaaaaggagcgg                                                                                                                          >  2:225540/2‑232 (MQ=255)                                                                                              gTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑Cggg                                                                                                                        >  1:6508/1‑250 (MQ=255)                                                                                                cGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑c                                                                                                                     <  1:126546/251‑1 (MQ=255)                                                                                                  aTCAGGATGCACTGAAAGCAAAATTTGCCGGAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAATATCTTGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑Ctt                                                                                                                   <  1:170099/251‑1 (MQ=255)                                                                                                     aGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTg                                                                                                                  >  1:359574/1‑249 (MQ=255)                                                                                                      ggATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGCCGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGAAATGCAGCGGATCGGCGTGAC‑CGATTACACCATTTTCAGCACGGTAACAGGTG‑CGGAG‑CTTGAgt                                                                                                               >  2:62424/1‑250 (MQ=255)                                                                                                           cACTGAAAGCAAGATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGGTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAACACGGTAAAAGGTG‑CGGAG‑GTTGAGgtgctg                                                                                                          >  1:259725/1‑251 (MQ=255)                                                                                                              tGAAAGCAAAATTTCCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTgcgc                                                                                                       >  1:311499/1‑251 (MQ=255)                                                                                                               gAAAGCAAAATTTGCCGGAAGCAACGGTAACGGCTCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACGCTGCCTGCGAACCGCGCAATCTCTATTGCACTTGATTTCGACGATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGGAATGCAGCGTATCGGCGTGAC‑CTATTACACCAGTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGGTGCTgcg                                                                                                        <  2:378096/249‑1 (MQ=255)                                                                                                                    caaaattttccgatatcaactgtAACGGCCCAATCTCGCTGGTGATCTGGACCTCTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGAGTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGG‑TTTAca                                                                                                <  2:157882/229‑2 (MQ=255)                                                                                                                            tccGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGC‑TTTACCCATCCGt                                                                                         <  2:168669/249‑1 (MQ=255)                                                                                                                               gTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGC‑TTTACCCATCCGTTTAt                                                                                     <  1:103580/251‑1 (MQ=255)                                                                                                                                     aaCGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑AGATTAAACCATTCTCCACACGGAAAAAGGTGACGG‑G‑CTTGAGTTGCTTCGC‑TTTACCCATCCGTTTATGGGCtt                                                                               >  2:273366/1‑251 (MQ=255)                                                                                                                                      aCGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGC‑TTTACCCATCCGTTTATGGGCTTc                                                                              <  1:368873/251‑1 (MQ=255)                                                                                                                                       cGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAACGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTACAAGGGG‑CGGAG‑ATTGAGTTGCTGCGC‑TTTACCCATCCGTTTATGGACTTCg                                                                             >  1:42938/1‑251 (MQ=255)                                                                                                                                                   aaTCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGAGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCCTTCTCAGAACGGTAAAAGGAG‑CGGGG‑TTTGAGATGCCGCGC‑AGTACCCATCCGTTTATGGGCCTCGACGTACCGGcc                                                                  >  2:145304/1‑249 (MQ=255)                                                                                                                                                      ctcGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGATAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCGGCACGGGAAAAGGTG‑CGGAG‑CTTGAGCTGCTGCGC‑TTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATcc                                                              >  1:84614/1‑251 (MQ=255)                                                                                                                                                         gCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGAGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGTCGCAACGCG‑TTACAACATTCTCAGCACGGAAAAAGG‑G‑ACGAGATTTGTGTTGGTGGGC‑TTTACCCATCCGTTTATGAGCtttggcggtctggaactacacg                                                           >  2:303397/1‑231 (MQ=255)                                                                                                                                                                              gCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGC‑TTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATg                                      <  1:93284/251‑1 (MQ=255)                                                                                                                                                                                              gcgaaCCGCGCAATCTCTATTGCATCTGATTTCGACTATGCGCTTGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGC‑TTTACCCATCCGTTTATGGGCTGCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCg                       <  1:148750/247‑1 (MQ=255)                                                                                                                                                                                              gccAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCGGCACGGTAAAAGGTG‑CGGAG‑CTTGAGCTGCTGCGC‑TTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGt                      >  1:189423/1‑251 (MQ=255)                                                                                                                                                                                                     gcgcAATCGCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGC‑TTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACAcc               <  1:160371/251‑1 (MQ=255)                                                                                                                                                                                                        cactctcTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGC‑TTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACcgcg            <  1:235491/248‑1 (MQ=255)                                                                                                                                                                                                           tctcTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCCGGCCGTGTTTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGAGTGAC‑AGATAACACCAATCTCCGCAAGGGAAAAGGGG‑TGGAG‑CGTGAGTTGCT‑CGCCTTCACCCATCCGTTTAAGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGATACCAGAGATGCCTGTACCgatcccattcccacccgcacc         >  2:310830/1‑231 (MQ=255)                                                                                                                                                                                                                   gCACCTGACTTCGAGTATGAGTTGGTGCAGATCGACGGTCAGGTCGTGATTCTGGCGAAAGAGCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACTCCATTCCCAGCCCGGTAAAAGGTG‑CGGAG‑CTTGAGTGGCTGCGC‑TTTGCCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTGCCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCAACg  <  2:283002/250‑1 (MQ=255)                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                         CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCGGCACGGTAAAAGGTG‑CGGAG‑CTTGAGCTGCTGCGC‑TTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACG  >  NC_000913/22946‑23405

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