Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	254,249	A→T	intergenic (+1088/+10)	yafP → / ← pepD	putative N‑acetyltransferase YafP/peptidase D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	254,249	0	A	T	98.6%	270.2 / NA	69	intergenic (+1088/+10)	yafP/pepD	putative N‑acetyltransferase YafP/peptidase D
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (28/40);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAC  >  NC_000913/254004‑254486                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              tATTAGCGTGAAGGTTCAGTCAGAGCGTAGTAAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAACTAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTACATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGAATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcatatca                                                                                                                                                                                                                                           >  1:402521/1‑250 (MQ=255)  aTTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcatatcaaa                                                                                                                                                                                                                                         <  1:38840/251‑1 (MQ=255)         tGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACACGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATGGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAAAAATTAc                                                                                                                                                                                                                                  >  1:495302/1‑251 (MQ=255)            aGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTc                                                                                                                                                                                                                               <  1:489463/251‑1 (MQ=255)             ggTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCg                                                                                                                                                                                                                              >  1:22302/1‑251 (MQ=255)             ggTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCg                                                                                                                                                                                                                              <  2:554056/251‑1 (MQ=255)                tcagtcagAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTAc                                                                                                                                                                                                                                  <  1:416070/244‑1 (MQ=255)                tcagtcagAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTAc                                                                                                                                                                                                                                  >  2:416070/1‑244 (MQ=255)                          cGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtttcttt                                                                                                                                                                                                                 <  1:16943/251‑1 (MQ=255)                               tCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagc                                                                                                                                                                                                             >  1:214661/1‑250 (MQ=255)                                 aGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGt                                                                                                                                                                                                          <  1:363667/251‑1 (MQ=255)                                    aTGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtc                                                                                                                                                                                                        >  1:146222/1‑250 (MQ=255)                                    aTGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtc                                                                                                                                                                                                        >  1:54854/1‑250 (MQ=255)                                        tAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtc                                                                                                                                                                                                    >  1:114451/1‑250 (MQ=255)                                          aCAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagc                                                                                                                                                                                                 <  2:305532/251‑1 (MQ=255)                                            aaGCGGGTGGGACGATTGCGGATTGCATGGAAGCTGGAGGAACAGCAACCGGAAAATAGCGCGGCGCTGAAATCGCAGCAGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCCACGGACATTTACA‑GGGATGGCTTTTATCGAATGATAATGAAGGATGAAACTGCCGGAAGGCGAGTAAACGCCATCCGGCAGCATATCAAATAATTACATCGCCGGAATTTCTTTCAGGAGTTCAGTcagcag                                                                                                                                                                                               <  2:505100/251‑1 (MQ=255)                                                 gCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTcc                                                                                                                                                                                          <  1:70539/251‑1 (MQ=255)                                                    ggCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt                                                                                                                                                                                       <  2:47506/251‑1 (MQ=255)                                                            tGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCTCCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGGCAGCAGTGTCCAGTAATGAcc                                                                                                                                                                                <  2:618632/250‑1 (MQ=255)                                                              cTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTAc                                                                                                                                                                             >  1:516053/1‑251 (MQ=255)                                                                gATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTAc                                                                                                                                                                             <  2:189539/249‑1 (MQ=255)                                                                    gCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTc                                                                                                                                                                       <  1:336206/251‑1 (MQ=255)                                                                        ggAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCTGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGa                                                                                                                                                                     <  2:45792/249‑1 (MQ=255)                                                                        ggAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCGCCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCGTTCGAt                                                                                                                                                                    <  1:335622/250‑1 (MQ=255)                                                                        ggAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACCTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGAtg                                                                                                                                                                   >  2:580679/1‑251 (MQ=255)                                                                               ggAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACtt                                                                                                                                                            <  2:359846/251‑1 (MQ=255)                                                                                 agcaacagcaacAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagca                                                                                                                                                                                                            >  1:217584/1‑201 (MQ=255)                                                                                 agcaacagcaacAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcagca                                                                                                                                                                                                            <  2:217584/201‑1 (MQ=255)                                                                                       agcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATcc                                                                                                                                                    >  2:107583/1‑251 (MQ=255)                                                                                        gcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTTGATGTGAACTTGCTCATCCg                                                                                                                                                   <  1:481105/251‑1 (MQ=255)                                                                                                 aaaTAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGtgtg                                                                                                                                          >  1:489476/1‑251 (MQ=255)                                                                                                 aaaTAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGtgtg                                                                                                                                          <  1:214710/251‑1 (MQ=255)                                                                                                          ggcgCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGAAACCCGCGACGGACATTTACA‑GGGATGGCTTTTAACGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTAAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGACCTTGCTCATCAGGAGAGTGTGGCCCGGTg                                                                                                                                  >  2:6764/1‑250 (MQ=255)                                                                                                            cgcTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATa                                                                                                                               <  1:562979/251‑1 (MQ=255)                                                                                                               ggAAATCGCAGCGGAGAAGTTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGtt                                                                                                                            <  2:131012/250‑1 (MQ=255)                                                                                                                  aaTCGCAGCGATGACTGTTCCGTCACCAGATTGAACGTAGCATCCCGCGCCGACCCTTTACA‑GGGATAGATTTGATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGGGATTAAACGCCGTCCGGCAGCAACTAAAATAATTAATTCGCAGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGCTGTGAACTTGCTCATACGGTGAGTGTGGCCCGGTGATAGTTGGc                                                                                                                         >  2:208869/1‑251 (MQ=255)                                                                                                                     cGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGccc                                                                                                                       <  1:572392/250‑1 (MQ=255)                                                                                                                           ggcgAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa                                                                                                                 >  2:301061/1‑250 (MQ=255)                                                                                                                              gAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGaaccc                                                                                                             <  1:544037/251‑4 (MQ=255)                                                                                                                              gAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAcc                                                                                                             >  1:504595/1‑251 (MQ=255)                                                                                                                                aTGTCCCATCACCAGATTGAACGTGGCAAACCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt                                                                                                           >  2:97947/1‑251 (MQ=255)                                                                                                                                  gTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGt                                                                                                         <  1:386673/251‑1 (MQ=255)                                                                                                                                       aTCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGAGAATGAAGGATGAAACTGCCGGAAGGGGATTAAACGCCATCCGGCAGCATATCAAATAATTGCTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTTCGCTTTCGATGTGAACTTGCTCATCCGGTGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTACAtt                                                                                                    <  2:402521/251‑1 (MQ=255)                                                                                                                                         cACCAGATTGAACGTGGCAACCCGGGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCAGCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAGCCATGTCCAttt                                                                                                   <  2:599953/250‑1 (MQ=255)                                                                                                                                           ccAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCg                                                                                                <  1:107583/251‑1 (MQ=255)                                                                                                                                              gATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTCCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCACTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGAt                                                                                             <  1:250192/251‑1 (MQ=255)                                                                                                                                               aTTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATa                                                                                            <  2:22302/251‑1 (MQ=255)                                                                                                                                                  gAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACg                                                                                          <  1:38902/250‑1 (MQ=255)                                                                                                                                                         gCAACCCGCGACGGACATTT‑CACGGGATGGCTTTTATCGAAGGATACTGCAGGATGAACCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATCCGGTTTTTTg                                                                                  >  2:161616/1‑251 (MQ=255)                                                                                                                                                            aCCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAAc                                                                               <  1:294961/251‑1 (MQ=255)                                                                                                                                                               cgcgACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAAc                                                                               >  1:594310/1‑248 (MQ=255)                                                                                                                                                                  gACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGAcc                                                                          >  2:60285/1‑250 (MQ=255)                                                                                                                                                                   aCGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACcac                                                                        <  1:180365/251‑1 (MQ=255)                                                                                                                                                                    cGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATCCGGTTTTTTGAACAGACcac                                                                        >  2:267128/1‑250 (MQ=255)                                                                                                                                                                       aCATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCa                                                                                                                                                                                         >  1:36977/1‑134 (MQ=255)                                                                                                                                                                       aCATTTACA‑GGGATGGCTTTTATCGAAGGATAACGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCa                                                                                                                                                                                         <  2:36977/134‑1 (MQ=255)                                                                                                                                                                         aTTTACA‑GGGATGGCGTTTAGCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCTGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGCGAGTTTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTGTTTTGAACAGTCCACATTCCa                                                                  <  2:106632/251‑1 (MQ=255)                                                                                                                                                                           ttACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCTATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACGTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGGTTTTTGAACAGACCACATTCCAg                                                                 <  2:579972/250‑1 (MQ=255)                                                                                                                                                                               a‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCcg                                                            >  1:263268/1‑251 (MQ=255)                                                                                                                                                                                           ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAAt                                                  <  2:516053/250‑1 (MQ=255)                                                                                                                                                                                              tCGAAGGATAATGAAGGATGACACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATAAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTgg                                              >  2:65111/1‑251 (MQ=255)                                                                                                                                                                                                    gATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACCTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTc                                        <  1:580679/251‑1 (MQ=255)                                                                                                                                                                                                       aaTGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGc                                     <  1:1479/251‑1 (MQ=255)                                                                                                                                                                                                          gAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTc                                  >  1:322874/1‑251 (MQ=255)                                                                                                                                                                                                          gAAGGATCAAACTGCCTGAAGGCGATTAAACGCAATCCGGCCGCATATCAAATAATTACTTCGCCGTAATTTCTTTCAGCAGTTCAGACAGCAGTGTCCAGTAATGACCTACGCTTTCGATGGGAAGTGGCTCATCCGCAGAGTGTGGACCGGTGATAGTTGTCCCGATAGAAATCATGTCCATTTCCGGATACGGTTTTTTGAACAGTCCACAGACCAGACCCGCGTGGATTATCTGAATGtccgaggc                                   <  2:588942/250‑8 (MQ=255)                                                                                                                                                                                                                  aaaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGCAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATATGGATGTTCGGCGTCTTGTTGa                           >  2:338774/1‑250 (MQ=255)                                                                                                                                                                                                                       gCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAg                       <  1:109216/249‑1 (MQ=255)                                                                                                                                                                                                                            aaGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTg                 >  1:294861/1‑250 (MQ=255)                                                                                                                                                                                                                               gCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAg             <  1:89457/251‑1 (MQ=255)                                                                                                                                                                                                                                         gCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGt    >  1:304028/1‑250 (MQ=255)                                                                                                                                                                                                                                          ccATCCGGCAGTATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAc  >  1:84455/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              TATTAGCGTGAAGGTTCAGTCAGAGCGTAGTCAGCATGCTAACAAGCGGCTGGCACGATTGCTGATTGCCTGGAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACA‑GGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTGGATGTTCGGCGTCTTGTTGAACAGGCGCTGATAGGTTTCACGTAC  >  NC_000913/254004‑254486

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