Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,728,754	C→T	A136A (GCC→GCT)	rnt →	ribonuclease T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,728,754	0	C	T	98.4%	204.2 / NA	63	A136A (GCC→GCT)	rnt	ribonuclease T
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (42/20);  minor base A (1/0);  total (43/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGAAGAGGTGTAATCGA  >  NC_000913/1728514‑1728998                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     cacaTTACATTGCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAAc                                                                                                                                                                                                                                            <  1:237212/251‑1 (MQ=255)  acattacattTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACGGAAACg                                                                                                                                                                                                                                           <  1:126218/251‑1 (MQ=255)                 ggcGAACCATTTGACGGCCCAAATTTGCACCCAGAAGCCCCCGCCTTCAACGTCAGTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCTAGCGGCGGTATCCGCGCCATTATGGTGGCGCAAAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCAccc                                                                                                                                                                                                                            <  2:17642/249‑1 (MQ=255)                    gAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCttt                                                                                                                                                                                                                         >  1:476794/1‑251 (MQ=255)                               cGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTA‑CCCGTTCCACCCTTTCGCCACTTTTg                                                                                                                                                                                                              >  1:77746/1‑250 (MQ=255)                                    cAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACAct                                                                                                                                                                                                         >  2:120619/1‑251 (MQ=255)                                         ttccAACCAGGAGCCCTCGCCTTCAACGACATTGCCCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAATTTGGACGTAAAGGTATTAAAGCGAGCGTCTGTATCCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCTACCCTTTCGCCACTTTTGACActgctg                                                                                                                                                                                                     <  2:501801/247‑1 (MQ=255)                                            cAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACt                                                                                                                                                                                                 >  2:171914/1‑251 (MQ=255)                                             aaCCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTg                                                                                                                                                                                                >  1:411531/1‑251 (MQ=255)                                             aaCCAGAAGCACTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTg                                                                                                                                                                                                >  1:126539/1‑251 (MQ=255)                                              aCCAGAAGACCTAGCCTTCAACGGCATTGACCCGAACGATACCGATCGCGGCGCGGTCAGCGAAAACGAGGCGCTGCACGAAATTTTTACAGTTGTCCGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACCATGCCAATTTTGATCACCGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTACACTgg                                                                                                                                                                                               >  2:474390/1‑251 (MQ=255)                                                  gAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCgg                                                                                                                                                                                           >  1:153774/1‑251 (MQ=255)                                                         tCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGc                                                                                                                                                                                     >  1:19861/1‑250 (MQ=255)                                                           cccTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTACCCGCGCCATTAGGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTA‑CCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTc                                                                                                                                                                                 <  2:77746/250‑1 (MQ=255)                                                                cAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGAc                                                                                                                                                                             >  1:36174/1‑251 (MQ=255)                                                                 aaCGGCATTGACCCGCCCGATCCCGATCGCGGCGCGGTCAACGAAGACGAGACGCTGCACGAAATTTTTAAAGTTGTCCGTAACGGTATTAAAGCGAGCGGCTGTAACCGCGCCCTTATGGTGACGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTCACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACCGGCCGGGCTGGCACTCGGAc                                                                                                                                                                             >  2:515597/1‑250 (MQ=255)                                                                  aCGGCCTTGACCAGAACGATCCCGATCGCGGCGCGGTATGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGa                                                                                                                                                                              <  2:468154/248‑1 (MQ=255)                                                                         tGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTAt                                                                                                                                                                    >  2:143401/1‑251 (MQ=255)                                                                              cGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGCCAAACCGTATTGTc                                                                                                                                                                >  2:426847/1‑250 (MQ=255)                                                                                  cGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAgg                                                                                                                                                           >  1:231280/1‑251 (MQ=255)                                                                                         gATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCtttcc                                                                                                                                                     >  2:442817/1‑247 (MQ=255)                                                                                          aTCGCGGCGCGGTCAGCGAATCCGGGGCGCGGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTAGGGTGGCGCACAATGCCAATTTTGATAACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCGGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAg                                                                                                                                                   <  2:116445/251‑1 (MQ=255)                                                                                                         gCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTTTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTAGGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGAc                                                                                                                                    <  2:124003/251‑1 (MQ=255)                                                                                                           gAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACtt                                                                                                                                  >  1:313178/1‑251 (MQ=255)                                                                                                                  aGGCGCGGCACGAAATTTTTAAAGTTGTACGTTAAGGTATTAAAGCGAGCGGATTTACCCGCGCCGTTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACGTTTGACACTGCTGCACTGGCCGTGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGc                                                                                                                           <  2:2688/251‑1 (MQ=255)                                                                                                                    gcgcTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCAc                                                                                                                         <  1:79114/251‑1 (MQ=255)                                                                                                                        tGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCAccc                                                                                                                       >  1:458560/1‑249 (MQ=255)                                                                                                                           aCGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGgcg                                                                                                                  >  1:355982/1‑251 (MQ=255)                                                                                                                            cGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGgcgc                                                                                                                 >  2:157381/1‑251 (MQ=255)                                                                                                                                aTTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCATAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTc                                                                                                             >  1:72030/1‑251 (MQ=255)                                                                                                                                aTTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGTGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCCGACCGCTGGCATGGACTTCGCAAGCACCCAGGCGCACTc                                                                                                             >  2:121281/1‑251 (MQ=255)                                                                                                                                   tttAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGAATGGACTTCGACAGCCCCCAGTCGCACTCcgc                                                                                                          >  2:79796/1‑251 (MQ=255)                                                                                                                                           ggTACGTAAAGGTAGTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTAc                                                                                                   <  2:411531/249‑1 (MQ=255)                                                                                                                                                gTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCGCCCTTTCGCCTCTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATGGTCAAAGGCTTGCCAGTCCGCTGGCAGGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACAcc                                                                                             <  1:220391/251‑1 (MQ=255)                                                                                                                                                gTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACAcc                                                                                             <  1:410063/251‑1 (MQ=255)                                                                                                                                                  aaaGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCg                                                                                            >  1:505774/1‑250 (MQ=255)                                                                                                                                                    aGGTATTAAAGCGCGACGCTGTAAACGCGACCTTATGGTCGCGCACAATGCAACTATTGATCACAGCTCTCTGATGGCCGCCGACGAACGCGAATCACTGCAACGTAACCCGTTCCGCCCTTTCGCCAATTCTGACCATGCTGCCCTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCACAGGATTGCCAGACCGCTGGCCTGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGaa                                                                                          >  2:506417/1‑250 (MQ=255)                                                                                                                                                    aGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTCTGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCAACAGCACCAAGGCGCACTCCGCGCTGAACGAAACCGAAc                                                                                         >  2:111064/1‑251 (MQ=255)                                                                                                                                                       tATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGCCACCGTACGc                                                                                       >  1:182935/1‑250 (MQ=255)                                                                                                                                                        aTTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCAc                                                                                     >  1:272537/1‑251 (MQ=255)                                                                                                                                                         aaaaggtgagcgGTTGTAACGGCGTCAGTATGGTGACGCCCTCGGCCAATTTGGATCACATCTGTATGAGGGCCGCCTCAGAACGCGCTTCACTGAAACGTAACCCGTGCCACCCTTTCGCCACTGTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTGGTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCCCCCAGGCTCACTCCGCGCTGTACGACACCGAACGCAct                                                                                    <  2:79706/244‑1 (MQ=255)                                                                                                                                                          tAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACACTGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACATCACCCAGGCGAACTCCGCGCTGTACGACACCGAACGCActg                                                                                   >  2:367456/1‑251 (MQ=255)                                                                                                                                                          aaaaGCGAGCGGCGGTACCCGCGCCATTATGGTGGCGCAAAATGCCAATTTTGATCACTGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCActg                                                                                   <  2:395152/250‑1 (MQ=255)                                                                                                                                                           aaaGCGAGCGGCTGTAACCGCGCCATTATGGTTGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGACACGTAACCCGTTCCGCCCTTTCGCCACTTTTGACACCGCTGCACTTGCCGGGCTGGCACTCGGACAAATCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCActgc                                                                                  <  1:67305/251‑1 (MQ=255)                                                                                                                                                           aaaGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCActg                                                                                   <  1:83854/250‑1 (MQ=255)                                                                                                                                                                        gTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTtgt                                                                     >  2:138262/1‑251 (MQ=255)                                                                                                                                                                               cgcCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGATGTGCTGTTTTGTGAAATCg                                                              >  2:126504/1‑251 (MQ=255)                                                                                                                                                                                           ggCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGa                                                  <  1:444111/251‑1 (MQ=255)                                                                                                                                                                                                         aTTTTGATCACAGCTTTATGATGGCCGCCGCCGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGc                                    >  1:169857/1‑251 (MQ=255)                                                                                                                                                                                                         aTTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAAc                                                         >  2:136566/1‑230 (MQ=255)                                                                                                                                                                                                         aTTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAAc                                                         <  1:136566/230‑1 (MQ=255)                                                                                                                                                                                                             ggATCACAGCTTGATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTGGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGc                                <  2:187196/250‑1 (MQ=255)                                                                                                                                                                                                               aTCACAGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGTGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCGGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTCCGACACCGAACGCACGGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGGCTGGGAGGcgtgccg                              >  2:120457/1‑245 (MQ=255)                                                                                                                                                                                                                  acaGCTTTATGATGGCCGCCGCAGAACGCGCTTCAATGATACGTAACCCGTTCCACACTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTACCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCAGCGATGTACGACACCGAACGCACTGCTGTGCTGTTTTTTGAAATCGTAAACCGCTGGCAACGACTGGGCGGCTGGCCGcgc                           >  2:432029/1‑249 (MQ=255)                                                                                                                                                                                                                    aGCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATc                         >  1:306609/1‑251 (MQ=255)                                                                                                                                                                                                                     gCTTTATGATGGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTAACAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGACATCGTCACCCGCTGGAAACGTCTGGGAGGCTGCTCGCTATCt                        >  2:234834/1‑251 (MQ=255)                                                                                                                                                                                                                             atgGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGccgccg                 >  1:500347/1‑250 (MQ=255)                                                                                                                                                                                                                             atgGCCGCCGCAGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGccgccg                 >  1:275678/1‑250 (MQ=255)                                                                                                                                                                                                                               ggCCGCCGCAGAATGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGa                >  2:138847/1‑249 (MQ=255)                                                                                                                                                                                                                                        aGAACGCGCTTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGAAGAGGTGTAAt     >  1:90383/1‑251 (MQ=255)                                                                                                                                                                                                                                         gAACGCGCTTCACTGAAACGTAACCCGTTCCACCATTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGAATTCGACAGCACCCAGGCGAACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAACCGTCTGGGAGGCTGGACGCTATCTGCCGCCGAAGGGGTGTCATc    >  2:388919/1‑251 (MQ=255)                                                                                                                                                                                                                                          caCGCGCTTCACTGAAACGTAAACCGGTACACACTTTCGCCACTTTTGACACTGCTGCACTGACCGGGCTGGCACTAGGACCAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAC‑CGTCTGGGAGGCTGGCCGCTATCTGCCGCCTAAGAGGTGTAATCg   >  2:98666/2‑250 (MQ=255)                                                                                                                                                                                                                                           aCGCGCTTCACTGAACCGTAACCCGTTCCACCCTTTCGCCACTGTTGCCAACACTGAACCGACCGCACTGGAACACGGACAAACCGTATTGTCAAAGGATTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGAACTCCGCGCCGTACGCCCCCGAACGCACTGATGAGCGGAATTGTGAAATCGACAACCGCTGGAAACGTCTGGGAGGACGGCAGCTATCTGCCGCCGGAGCGGAGTAATctc  >  2:222155/1‑249 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     CACATTACATTTCCACGTCGAACCATTTGTCGGCGCAAATTTGCAACCAGAAGCCCTCGCCTTCAACGGCATTGACCCGAACGATCCCGATCGCGGCGCGGTCAGCGAATACGAGGCGCTGCACGAAATTTTTAAAGTTGTACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGTCAACCGCTGGAAACGTCTGGGAGGCTGGCCGCTATCTGCCGCCGAAGAGGTGTAATCGA  >  NC_000913/1728514‑1728998

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