Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,878,947	0	C	T	26.0%	99.0 / 16.1	51	E276K (GAA→AAA)	cas1	multifunctional nuclease Cas1
Reads supporting (aligned to +/- strand):  ref base C (5/32);  new base T (13/0);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 2.41e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCC  >  NC_000913/2878703‑2879194                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                               ccAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCATTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGc                                                                                                                                                                                                                                                         <  1:445387/251‑1 (MQ=255)        aGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGc                                                                                                                                                                                                                                                  <  1:421252/251‑1 (MQ=255)         gCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTTTTTTTC‑CAGCGGCaaac                                                                                                                                                                                                                                                  >  1:19569/1‑248 (MQ=255)             tttGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGGATTTTCC‑CAGCgggaaagacgcc                                                                                                                                                                                                                                             >  1:134723/1‑241 (MQ=255)                gTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCctc                                                                                                                                                                                                                                          <  2:123863/251‑1 (MQ=255)                  tCCCAGATCATTTCACGAATTTTTTCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCctct                                                                                                                                                                                                                                         <  2:427739/250‑1 (MQ=255)                   cccAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCtat                                                                                                                                                                                                                                       <  2:19569/251‑1 (MQ=255)                      aGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGAGGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTGCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAg                                                                                                                                                                                                                                    <  2:508758/251‑1 (MQ=255)                       gATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATCCTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGc                                                                                                                                                                                                                                   <  1:56991/251‑1 (MQ=255)                          cATTTCACGAATTTTTGCGGATACCTCACCTACATATACCCCTGCACGTCCCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTGGTGTTTTTC‑CACCGGCCAGccccgcctctgttcacggg                                                                                                                                                                                                                                >  2:41581/1‑233 (MQ=255)                          cATTTCACGAATTTTTGCGGAGACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGAGGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACGCCGATGGCCTGCATCGCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGGTGTATTTCTC‑CATCGGCAAGCACGTCCTCTATAAGCgg                                                                                                                                                                                                                                 <  2:253514/250‑1 (MQ=255)                               cACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCaa                                                                                                                                                                                                                          <  2:360151/251‑1 (MQ=255)                                        ttGCCGATACATCACCTACATATACCCCTGCACGTACCTCTCACAACCACATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGATACTCCGCTGGCCTGCATCTCCCAGCGAAACAGGACGCGGACTGGCAACAGGCTGAGCATCTTCAGGTGGGGCCGGCGGTTGGATTT‑TCTCAGCGGCCAGCAAGGCCACTATAAGCGG‑GATTAAttttgcca                                                                                                                                                                                                                  >  1:323288/1‑245 (MQ=255)                                          gCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAAt                                                                                                                                                                                                                <  2:104482/250‑1 (MQ=255)                                            ggATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTC‑CCGCGGCAAGAACGTCCTCTATAAGCGG‑GATCAATTTTGCTAATG‑tt                                                                                                                                                                                                            >  1:164466/1‑251 (MQ=255)                                            ggATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCGCTATAAGCGG‑AATCAATTTGGCTAATG‑tt                                                                                                                                                                                                            <  1:383568/251‑1 (MQ=255)                                                    cACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTactact                                                                                                                                                                                                    <  2:164466/251‑1 (MQ=255)                                                         aCATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTT‑CAGCGGGAAGCACGGCCCCTATAAACGG‑AATCAATTTGGCTAATT‑TTTTTCTACTGCG‑a                                                                                                                                                                                               >  1:166441/1‑250 (MQ=255)                                                         aCATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCg                                                                                                                                                                                                 <  1:387003/249‑1 (MQ=255)                                                             ataCCCCTGCACGTACCTCCAACAACCAGATGGCGAATCTGCCTCGTAAGCGCGGAGTTACATTTTCAGTGACCACGACCAACATACTCATGTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGTCAACAGGCTGTGCGTCTTCAGGTGGGGGCGGCGGTTGTAGTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTGGGCTAATG‑TTTTACTTCTGCG‑GAAAat                                                                                                                                                                                          <  2:44768/251‑1 (MQ=255)                                                                ccccTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTC‑CCACGGCAAGCCCGGCCTCTATAAACGG‑AATCAATTTGGCTAATG‑TTTTTCTACTGCG‑GAAAATATc                                                                                                                                                                                       >  2:147985/1‑251 (MQ=255)                                                                        cGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTCTTTTTC‑CAGCCGCAAAGACGGCCTCCATAAGCGG‑AATCAATTTGGGTCATGGTTTTA‑TACTGCG‑GAAAAAATACCTGCAcc                                                                                                                                                                               >  1:228303/1‑250 (MQ=255)                                                                               ccAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGCACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCCGCATCTCCCAGTGCAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CACCGGCAAGCACGGCCTCTATAAGCGG‑AATCAATTTTGCTAATG‑TTTTACTACTGGG‑AAAAATAACCCTGCACGCCCAAcc                                                                                                                                                                        >  1:497518/1‑250 (MQ=255)                                                                                caacaaCCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACgg                                                                                                                                                                       <  2:154402/251‑1 (MQ=255)                                                                                caacaaCCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGGGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCGCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACgg                                                                                                                                                                       <  1:462060/251‑1 (MQ=255)                                                                                 aacaacCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCGCTATAAGCGG‑AATCAAt                                                                                                                                                                                                                         <  1:377119/202‑1 (MQ=255)                                                                                             gCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCa                                                                                                                                                                                                                                                 <  1:223259/167‑1 (MQ=255)                                                                                                aaTCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATATTCAGGTGGGGCCGGCGGTTGTATTTTTC‑CAGCGGCAAGAAAGACCTTTATAAGCGG‑AATCAAATTTGCTAATG‑TTTTACTACTGCG‑GAAAATATCCCCTCACGCCCAAAGGACTTCCCGGGCCGGcc                                                                                                                                                       >  1:382738/1‑250 (MQ=255)                                                                                                   cTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCCCGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTC‑CCCCGGCCAGCACGTCCCCTCTTAGCGGGAACCAATTTG‑CTTATG‑TTTTTCTACTTCG‑CAAAATATTCCTCCACCCCCAACAGACTTccccggcccgctccc                                                                                                                                                    >  2:192404/1‑239 (MQ=255)                                                                                                      ccTCGGAAGCGCGGAGGTACATTGTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACATGAAGCGGAATGTCAACAGGCTGTGCATCTTCGGGGGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGTGG‑AATCAATTTGGCTAATG‑TTTTACTACTTCG‑AAAAATAGCCCTGCACGCTAAACGGACGTCCCGGTCCGGCTCACC‑A‑g                                                                                                                                               <  1:355369/251‑1 (MQ=255)                                                                                                                          aTTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTC‑CAGCGGCAAGCAAGTCCTCTATAAGCGG‑AATCAATTTGGCTCATG‑TTTTACTACTTCG‑AAAAATATCCCTGCACGCCCAACGGACTTCCCGGTCCCGCTCACC‑A‑GGGTTACGACGCGCTAtctcc                                                                                                                           >  1:102922/1‑250 (MQ=255)                                                                                                                            tttCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCaa                                                                                                                          <  1:237568/250‑1 (MQ=255)                                                                                                                                       aCGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGGCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTa                                                                                                              <  2:318208/251‑1 (MQ=255)                                                                                                                                         gACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTTTC‑CAGCGGCAAGAACGACCTCTATAAACGG‑AATCAATTTTGCTTATG‑TTTTTTTACTGCG‑AAAAATATACCCGCACGCCCAACGGACTTCCCGGTCCGGCTCACC‑ACGGGTT‑CGACGCGCTATCTCAAAAACTTTCGGttc                                                                                                             >  1:514205/1‑248 (MQ=255)                                                                                                                                                           aGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTGACGACGCGCTATCGCAAAAGCTTTCTGTGCAACAGTGT‑CAAATTTAatg                                                                                         <  2:269782/251‑1 (MQ=255)                                                                                                                                                            gCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CCGCGGGAAGCACGTCCCCTATAAGCGG‑GATCAATTTTGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTTCCCGCCCAACGGACTTCCCCGTCCGGCTCACCCA‑GGGTT‑CGAAGCGCTATCTCAAAAGCTTTCGGTACAACAggggccaattttatg                                                                                          >  1:259247/1‑236 (MQ=255)                                                                                                                                                               aCTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAatgat                                                                                       <  1:50589/249‑1 (MQ=255)                                                                                                                                                                      gggCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAat                                                                              <  2:166441/250‑1 (MQ=255)                                                                                                                                                                                  cccAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGAGGTCGGCAATATCGTAAACaa                                                                   <  2:497518/250‑1 (MQ=255)                                                                                                                                                                                   ccAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGAATTTTTC‑CAGCGGCACGCACGTCCTCTATAAGCGG‑AAACAATTTGGCTAATG‑TTTTACTACTGCG‑GAAAATATCCCTGCACGCCCAACAGACATCCCGGTCCTGCGCACC‑A‑GGGCTACGACGCCCTCTCTCAAAAGCTTTCGGTACAACAGTGT‑CAAACTTAATGATGTCTGCAATATCGTAAACaaa                                                                  >  1:244277/1‑250 (MQ=255)                                                                                                                                                                                         aaaCAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGATGTATTTCTC‑CCGCGGCCAACACGTCCCCTTTAAGCGG‑AATCAATTTTGCTAATG‑TTTTACTACTGCG‑AAAAATATACCCGCACGCCCAACGGACCTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCCAAAGCTTTTGGTACAACAGGGT‑CAAATTTAATGATGTCTTCAATATCGTAAACAAAAGaaaa                                                            >  1:102013/1‑249 (MQ=255)                                                                                                                                                                                           aCAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTTTATTTTTC‑CAGCGGCAAGCACCTCCTCTATAAGCGG‑AATCAAATTTGCTAAAG‑TTTTTCTACTTCGGAAAAAATTCCCT‑CCCGCCCAACGGACTTCCCGGGCCGGCTCCCC‑A‑GGGGTACAACGCGCTATTTCCAAAGCTTTCCGTTCAACAGAGTACAAATTAAATG‑TTTCTGCTATAACGTAAAAAAaagaaaga                                                           >  2:83709/1‑249 (MQ=255)                                                                                                                                                                                                    ggAATGGCAACAGGCTGTTCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGt                                                <  2:382738/251‑1 (MQ=255)                                                                                                                                                                                                          gCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGcac                                          <  2:514205/251‑1 (MQ=255)                                                                                                                                                                                                                    gtgCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATa                                <  1:192404/251‑1 (MQ=255)                                                                                                                                                                                                                       cATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑GTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCt                             <  2:37456/251‑1 (MQ=255)                                                                                                                                                                                                                                  ggggCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAAc                   <  1:418489/250‑1 (MQ=255)                                                                                                                                                                                                                                    ggCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGAGG‑AATCAATTTGGCTCATG‑ATTTCCTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAAcc                  >  1:492678/1‑249 (MQ=255)                                                                                                                                                                                                                                        ggcggTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTgg            >  2:3102/1‑250 (MQ=255)                                                                                                                                                                                                                                            gTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGt        <  1:200228/251‑1 (MQ=255)                                                                                                                                                                                                                                                   ttCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATcgcc  <  1:41581/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                               CCAGTCCAGCTATTTGTTCCCAGATCATTTCACGAATTTTTGCGGATACATCACCTACATATACCCCTGCACGTACCTCCAACAACCAGATGGCTAATCTGCCTCGTAAGCGCGGAGGTACATTTTCAGTGACCACGACCAACATACTCATTTCAGCTACTCCGATGGCCTGCATCTCCCAGTGAAACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTC‑CAGCGGCAAGCACGTCCTCTATAAGCGG‑AATCAATTTGGCTAATG‑TTTTACTACTGCG‑AAAAATATCCCTGCACGCCAAACGGACTTCCCGGTCCGGCTCACC‑A‑GGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGT‑CAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCC  >  NC_000913/2878703‑2879194

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