Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,039,623	G→A	H41Y (CAT→TAT)	xerD ←	site‑specific recombinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,039,623	0	G	A	95.8%	141.0 / NA	48	H41Y (CAT→TAT)	xerD	site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/28);  minor base C (2/0);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

AATAAACGTTCGACCTGCGCTTCGCTTAAATCTTTTGGCAAACGCTGGGGCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATGCAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCATTATGCCTTATTGTGCCGTGACTAAAGCGATTCTGATACACTAGCCGCAAAAGCCACAGCAGAATCGAGAAGCTTACGTTATGAATATG  >  NC_000913/3039379‑3039863                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   aataaACGTTCGACCTGCGCTTCGCTTAAATCTTTTGGCAAACGCTGGGGCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCa                                                                                                                                                                                                                                              <  1:87922/251‑1 (MQ=255)    aaaCGTTCGACCTGCGCTTCGCTTAAATCTTTTGGCAAACGCTGGGGCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGGCGCCAGCGTCAACCCGCGGTGATACAACCACTc                                                                                                                                                                                                                                           <  1:262163/251‑1 (MQ=255)      aCGTTCGACCTGCGCTTCGCTTAAATCTTTTGGCATACGCTGGGGCATTTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGAGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAAAGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTATTGCCTGCAAATCGTCACATTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTcc                                                                                                                                                                                                                                          >  1:65521/1‑250 (MQ=255)           cGACCTGCGCTTCGCTTAAATCTTTTGGCAAACGCTGGGGCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCGCGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCTCACTGCACTCAGCAAACGCGCGGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTGTGCGCCGGCGCCAGCTTCAACCCGCGGTGATACAACCACTCCACcatc                                                                                                                                                                                                                                    <  1:454139/251‑1 (MQ=255)           cGACCTGCGCTTCGCTTAAATCTTTTGGCAAACGCTGGGGCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGTCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACcat                                                                                                                                                                                                                                     <  2:386819/250‑1 (MQ=255)            gACCTGCGCTTCGCTTAAATCTTTTGGCAAACGCTGGGGCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACcatca                                                                                                                                                                                                                                   <  1:346060/251‑1 (MQ=255)            gACCTGCGCTTCGCTTAAATCTTTTGGCAAACGCTGGGGCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACcatc                                                                                                                                                                                                                                    >  1:291287/1‑250 (MQ=255)                       cgcttAAAGCTTTGGGCAAACGCTGGGGCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCCCTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATc                                                                                                                                                                                                                        <  2:49061/250‑1 (MQ=255)                            aaaTCTTTTGGCAAACGCTGGGGCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGAc                                                                                                                                                                                                                   <  1:2357/250‑1 (MQ=255)                             aaTCTTTTGGCAAACGCAGGGGCAATTTCCGTCAAGCGAGATGCGCACTGGGACCGTATTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGAAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAACGACTGCAAATAGTAACTTTGCGCCGTAGCCAGCGTCAACCCGAGGCGATACAAACACTCCACCACCATTG‑ACAGATCGCGACgg                                                                                                                                                                                                                 >  2:385268/1‑251 (MQ=255)                                             cTGGGGCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTAt                                                                                                                                                                                                 >  1:331309/1‑251 (MQ=255)                                               ggggCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCTCCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTAtt                                                                                                                                                                                                <  1:31393/250‑1 (MQ=255)                                                  gCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑Ca                                                                                                                                                                                           <  1:222087/251‑1 (MQ=255)                                                    aaTTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGACGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACCGATCGAGCCGGTAAGCGTTCAACGAATTTT‑CAGc                                                                                                                                                                                         >  2:196578/1‑251 (MQ=255)                                                        tCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGcc                                                                                                                                                                                        >  1:200786/1‑248 (MQ=255)                                                         cGGTGAAGCGAGATGCGCACTGGGATGATCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCTCTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCGCCAGCCGTTCGGCCAGTAATGCCTGCAAATCGGCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAACGTTCAACGTATTTT‑CAGTCAGAt                                                                                                                                                                                    <  1:147568/251‑1 (MQ=255)                                                         cGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGAt                                                                                                                                                                                    <  1:19285/251‑1 (MQ=255)                                                          ggTGACGCGCGATGCGCACTGAGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGAGGCTTTATACCCGCCCTCCAGCCGTTCTGCAAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCAAGCGTCTACCCGCGGTGATACCACCACTCCACCATCATTGCACA‑ATCGCGACGGTAAGCGTTCAACGAAttttacgaccaatt                                                                                                                                                                                    >  2:249396/1‑241 (MQ=255)                                                             gAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGAttttt                                                                                                                                                                                <  1:166524/251‑1 (MQ=255)                                                               aGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGAttttttt                                                                                                                                                                              >  1:322816/1‑251 (MQ=255)                                                                                 aTCGTCTTCACTAACCTTTTCGCGATAAAGATACTGGAACAATCGGCGCCCTGCACCCAGCAAACTCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCTTCGCCAGCGTCATCCCGCGGTGATACAACCACTCCCCCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCaa                                                                                                                                                            <  2:436628/251‑1 (MQ=255)                                                                                   cGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGa                                                                                                                                                          <  1:416460/251‑1 (MQ=255)                                                                                    gTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGa                                                                                                                                                          >  2:178382/1‑250 (MQ=255)                                                                                         cACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGc                                                                                                                                                    >  1:386887/1‑251 (MQ=255)                                                                                            gAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGATCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAACCTGCTCg                                                                                                                                                 >  2:122328/1‑251 (MQ=255)                                                                                               aCTTTACGCGATAAAGATACTGTGTCAATCGGCGCTCTGCACTCAGCAATCTCGCTGCGCTGGGTGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAACCGTCACTTGTCGCTGGCGCCAGCTGCAACCCGCGGTTAGACAACCACTCCACCATCATGG‑ATAGTTCTGTACGGTATGCGTTATACGTAGTTT‑CAGCGAGATTTTTTTCCAGCTACAGAGCATCATGAAACTGCTCGATg                                                                                                                                              <  2:268156/251‑1 (MQ=255)                                                                                                    tcgcgaataaGATACTTGAACAATCGGCGCACTGCACCCAGCAAACGCGCTTGGCCGTGGGCTTTATACCCCCCCTCCAGCCGTTCTTCCAGGTATGCCTGCAAACCGTCCCTTTGCGCCGGCGCCAGCGGCAACCCGCGGTGATACAACCACTCCTCCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTg                                                                                                                                          <  2:244964/242‑1 (MQ=255)                                                                                                     cgcgATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGcc                                                                                                                                        >  1:330464/1‑251 (MQ=255)                                                                                                               aTACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGtt                                                                                                                              <  1:264922/251‑1 (MQ=255)                                                                                                                          aTCGGCGCACTGCACTCAGCAAACGCGCTGCGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGACACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATCCAACCACTCCACCATCATTG‑ACAGATCGCGCCGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGcc                                                                                                                   >  1:64588/1‑251 (MQ=255)                                                                                                                          aTCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGcc                                                                                                                   >  1:258883/1‑251 (MQ=255)                                                                                                                             ggCGCTCTGCACCCATCAAACGCTGTGAGCTGGTAGCTTTATCCCCGCCCTCCAGCCGTTCTTCCAGTAAGGCCTGCAAATCGTCACTTTGCGCAGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCTCAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTTCGCCCCt                                                                                                                <  2:336398/251‑1 (MQ=255)                                                                                                                                gcactgcactCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTAt                                                                                                             >  1:361898/1‑251 (MQ=255)                                                                                                                                     tcacccAGCAAACGCGCTGAGTTGGTGGATTATTACCCGCCCTCCAGCCCTTCTGCCAGTAATGCATGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAAACACTCCACCATCATTG‑ACAGATCGCGACGGTCAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGt                                                                                                          <  2:336154/244‑1 (MQ=255)                                                                                                                                               acacgccTGAGCTGGTGGATTTATACCCGCCCTCCAGCCGTTCTGCAAGTAATGCCTGCAAATCGCCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGAt                                                                                              <  2:10044/245‑1 (MQ=255)                                                                                                                                                   cgcTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATAAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCAt                                                                                          >  1:85740/1‑251 (MQ=255)                                                                                                                                                             gtggCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCATTATGCCTTAt                                                                                <  1:15937/251‑1 (MQ=255)                                                                                                                                                                gCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGAGATACAACCCCTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCATTATGCCTTTTtgt                                                                             >  2:109165/1‑251 (MQ=255)                                                                                                                                                                                   gCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGACACAACCACTCCACCATCATTG‑CCAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAACCTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCATTATGCCTTATTGTGCCGTGACTAAAGCGATTc                                                          >  1:258786/1‑251 (MQ=255)                                                                                                                                                                                                 aaggccggCAAAGCGTCCCTTTGCGCCGTCGCCAGCGTCAACCTGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTGCCAGCCACAGATCATCAAGAAACGGCTCGATGCGTGTCAGATCCTGTTTCTCTTGCGCTCCTTATGGTCACTCATTTGTTGCATTATGCCGTATTGTGCCGTGACTAAAGCGATTCTGATACACTAGcc                                             <  1:208997/243‑1 (MQ=255)                                                                                                                                                                                                        gCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTTTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCATTATGCCTTATTGTGCCGTGACTAAAGCGATTCTGATACACTAGCCGCaaaa                                       <  1:318101/249‑1 (MQ=255)                                                                                                                                                                                                                  aCTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCATTATGCCTTATTGTGCCGTGACTAAAGCGATTCTGATACACTAGCCGCAAAAGCCACAGCAGa                            >  1:48155/1‑250 (MQ=255)                                                                                                                                                                                                                       gcgcCGCCGCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCATTATGCCTTATTGGGCCGTGACTAAAGCGAGTCTGATACACTAGCCGCAAAATCCACAGCAGAATCgag                      <  2:171266/251‑1 (MQ=255)                                                                                                                                                                                                                              cGCCAGCGTCAACCCGCGGTGATCCAACCACTCCACCATCATTG‑ACAGATCGCGACGGCAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCCCTCATTTGCTCCATTATGCCTTATTGTGCCGTGACTAAAGCGATTCTGATACACTAGCCGCAAAAGCCACAGCAGAATCGAGAAGCtt                >  2:411878/1‑250 (MQ=255)                                                                                                                                                                                                                              cGCCAGCGTCAACCCGCGGTGATACAACCACTCCCCCATCATGG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTTACTTGCGCCCCTTATGGTCACTCATTTGATCCATTATGCCTTATTGTGCCGTGACTAAAGCGATTCTGATACACTTGCCGCAAAAGCCACAGCAGCATCGAGAAGCtt                <  1:249488/250‑1 (MQ=255)                                                                                                                                                                                                                               gCCAGCGTCAACCCGCGGTGATACAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCATTATGCCTTATTGTGCCGTGACTAAAGCGATTCTGATACACTAGCCGCAAAAGCCACAGCAGAATCGAGAAGCTTAc              <  1:197126/251‑1 (MQ=255)                                                                                                                                                                                                                                       cAACCCGCGGTGTTACAACCACTCCCCCATCATTG‑ACAGAGCGCGTCGTTAAGCGTTCAAC‑TATTTTTCAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCATTATGCCTTATTGTGCCGTGACTAAAGCGATTCTGATATACTAGCCGCAAAGGCCACAGCAGAATCGAGAAGCTTACGTTATGaa      <  2:291287/251‑1 (MQ=255)                                                                                                                                                                                                                                           ccGCGGTGATACAACCACGCCACCATCTTTG‑ACAGATCGCGACGGTAAGCTTTCAACGTAGTTT‑CAGCCAGATTTTTGTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCATTAGGCCTTATTGTGCCGTGACTAAAGCGATGCTGATACACTAGCCGCAAAAGCCACAGCAGAATCGAGAAGCTTACGTTATGAATATg  <  2:168926/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   AATAAACGTTCGACCTGCGCTTCGCTTAAATCTTTTGGCAAACGCTGGGGCAATTTCGGTGAAGCGAGATGCGCACTGGGATCGTCTTCACGAAACTTTTCGCGATAAAGATACTGGAACAATCGGCGCACTGCACTCAGCAAACGCGCTGAGCTGGTGGCTTTATACCCGCCCTCCAGCCGTTCTGCCAGTAATGCCTGCAAATCGTCACTTTGCGCCGTCGCCAGCGTCAACCCGCGGTGATGCAACCACTCCACCATCATTG‑ACAGATCGCGACGGTAAGCGTTCAACGTATTTT‑CAGCCAGATTTTTTTCCAGCCACAGAGCATCAAGAAACTGCTCGATGCGTGCCAGATCCTGTTTCACTTGCGCCCCTTATGGTCACTCATTTGATCCATTATGCCTTATTGTGCCGTGACTAAAGCGATTCTGATACACTAGCCGCAAAAGCCACAGCAGAATCGAGAAGCTTACGTTATGAATATG  >  NC_000913/3039379‑3039863

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