Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,331,059	2 bp→AT	coding (12‑13/2598 nt)	topA →	DNA topoisomerase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,331,059	0	T	A	97.8%	137.4 / NA	47	A4A (GCT→GCA)	topA	DNA topoisomerase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (24/21);  minor base C (0/1);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
*	NC_000913	1,331,060	0	C	T	97.9%	146.0 / NA	47	L5F (CTT→TTT)	topA	DNA topoisomerase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (24/22);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTTTCGTCAATCTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTAT  >  NC_000913/1330816‑1331296                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                              acgatttcgaccaTCTAGCAATAGACTTGCTTGACATACGACCAAAATTCCGTCGTGCTATAGCTCCGGGAGGCCAAGACCTGTTAACTCAGTCCTCTGAATTTGCTTTAACAGAGGCACTACAAGGTGTTGATAGCCGCAGAGATCGAGTCCATATCGGTAACACGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCCTTTG‑TCa                                                                                                                                                                                                                                        <  2:241647/239‑1 (MQ=255)               cTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCg                                                                                                                                                                                                                          >  1:38854/1‑251 (MQ=255)                   cAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCa                                                                                                                                                                                                                       >  2:460486/1‑250 (MQ=255)                   cAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGATTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGGGAATATGGGGAAAGCATTTG‑TCAGCGTTGAGTCCCCGGCa                                                                                                                                                                                                                       >  2:87535/1‑250 (MQ=255)                    aaTAGATTTGCTTGATATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGGTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGTCaa                                                                                                                                                                                                                      >  2:134078/1‑250 (MQ=255)                          tttGCTTGACATTCGACCAAAATTCCGTCGTGCTTTAGCGCCTGTAGGCCAAGCCCTGTTCACTCAGTCACTTGAATTTTCGTGAACAGAGTCACGACCAGGGGTTGATCTCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGCAGAAACAAATTAGGTAAGGAGAATAGTGGTAAAAC‑CTTGGTCATCGTTGAGTGCCCGGCAAAAGcaa                                                                                                                                                                                                                >  2:186431/1‑248 (MQ=255)                                            aaaaTTCCGTCGTGCTATAGCGCCTGGAGGTCCATACCCGGTAACTCAGGCGCCTGAATGTTCGTGAACAGAGTCACGTCAAGGGGTTGATATTCGCATAGCGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAAGCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAACACGATCAACAAGTATCt                                                                                                                                                                                             <  2:128576/251‑1 (MQ=255)                                                         gCTATAGCGCCTGTAGGCCAAGACCTGTTAACGCAGGCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGAGATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGACAAAACGATCAACAAGTATCTGGGTAGTGACTAc                                                                                                                                                                                <  2:38854/251‑1 (MQ=255)                                                                         gCCAAGACCTGTTAACTAAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGGTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGAAAGGTTTATCACAGGGAGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGAAAAGCATTTG‑TCATCGTTGAGTCCCAGGCAAAAGCCACAACGATCAACAAGTATCTGGGTAGTGACTCCGTGGTGAAATCAAGCg                                                                                                                                                                >  2:280591/1‑251 (MQ=255)                                                                            aaGACCTGTTAACTCAGTCGACTGAATTTTCGTGAACAGAGTCACGACCAGGGGTTGATATACGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTTGTGAAATCCAGCGTCg                                                                                                                                                             <  2:357860/251‑1 (MQ=255)                                                                               aCCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTc                                                                                                                                                          >  1:202365/1‑251 (MQ=255)                                                                                      aaCTCAGTCACCTGAATTGTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACGACGTGGTGAAATCCAGCGTCGGTCACATcc                                                                                                                                                    <  1:118894/250‑1 (MQ=255)                                                                                                 cTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTgg                                                                                                                                                                                           >  2:409235/1‑200 (MQ=255)                                                                                                 cTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTgg                                                                                                                                                                                           <  1:409235/200‑1 (MQ=255)                                                                                                           gTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTc                                                                                                                              <  2:425306/251‑1 (MQ=255)                                                                                                               aCAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTTGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGAGCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCt                                                                                                                          <  1:60815/251‑1 (MQ=255)                                                                                                                     tCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaa                                                                                                                     >  1:475443/1‑250 (MQ=255)                                                                                                                      cACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCaaa                                                                                                                    >  2:437592/1‑250 (MQ=255)                                                                                                                           cAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTTGGTAGTGACTACGTGGGGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGt                                                                                                              <  1:328091/251‑1 (MQ=255)                                                                                                                             aGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGc                                                                                                            >  1:139114/1‑251 (MQ=255)                                                                                                                                 gTTGATATCCGCAGAGAGCGAGTCCATATCGGAAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCAAATCCGCGATTTGCCGGCCAGTGGCTCAGCTGCCAAAAAGAATGCCgcc                                                                                                        >  2:431437/1‑249 (MQ=255)                                                                                                                                            cAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGGAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTccac                                                                                             <  2:202365/251‑1 (MQ=255)                                                                                                                                                gagCGAGTCCATATCGGTAACTCGTTGCCAGTGGCAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATAAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGAGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAACTGCCCAAAAGAGTGCCGACACTACCTCCACCaa                                                                                          >  1:296685/1‑250 (MQ=255)                                                                                                                                                      gTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCt                                                                                   <  1:437592/251‑1 (MQ=255)                                                                                                                                                      gTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGAGCCGACTCTACCTCCACCAAGACGGcc                                                                                   >  2:349211/1‑250 (MQ=255)                                                                                                                                                        ccATATCGGTAACTCGGTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTaa                                                                                 <  2:109569/251‑1 (MQ=255)                                                                                                                                                         cATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTaaa                                                                                <  1:374566/251‑1 (MQ=255)                                                                                                                                                                 tAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTa                                                                         >  1:366536/1‑250 (MQ=255)                                                                                                                                                                          gCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATAGGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGa                                                               <  1:354422/251‑1 (MQ=255)                                                                                                                                                                            cAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAAc                                                             <  2:475443/251‑1 (MQ=255)                                                                                                                                                                              gTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGt                                                           <  1:87625/251‑1 (MQ=255)                                                                                                                                                                                 gAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGgc                                                        >  1:449733/1‑251 (MQ=255)                                                                                                                                                                                 gAAGGGTTATCAACGTGCGACGCATTCCTGGAATAATCAAAGTAGGTAAGGTGAATATGGGTAAAGCATGTG‑TCAGCGTTGAGTCCCCGGCAAAAGCCGAAACGAGCAACAAGTATCTGGGTAGTGACTACGTGGTGTAATCCAGCTTCGGTCACATCCGCGATTTGCCGACCGGTGGCTCATCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGgc                                                        <  1:280591/251‑1 (MQ=255)                                                                                                                                                                                         aTCAACGTGCGACGCATTACTGGAAGAATCAAATTAGGTAAGGTGCATATGGGTAAAGCATTTG‑GCATAGTTGAGTCCCCGGCAAAAGCCAAAAGGATAAACAAGTCTCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGACACATCCGCGATTTGCCGACCAGTGGCTAAGCTGCCAAAAGGAGGGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGCTGAACGAGGCTCTCTCGt                                                >  2:436739/1‑251 (MQ=255)                                                                                                                                                                                                tGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCGCGTCAACCgg                                         >  2:102526/1‑250 (MQ=255)                                                                                                                                                                                                             tGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGAcc                            >  1:273085/1‑251 (MQ=255)                                                                                                                                                                                                                   aaTCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGc                      >  1:43290/1‑251 (MQ=255)                                                                                                                                                                                                                    aTCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGCCGGCTAAAAAGCCTAAAAAGGATGCACGTGGCGCTCTCGACAACCGTATGGGGGTTGACCCGTGGcc                     >  2:278090/1‑250 (MQ=255)                                                                                                                                                                                                                    aTCAAATTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGca                     <  2:196157/251‑1 (MQ=255)                                                                                                                                                                                                                    aTCAAATTAGGTAAGGGGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAGGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGca                     <  2:106818/251‑1 (MQ=255)                                                                                                                                                                                                                       aaaTTAGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACaa                  <  2:190036/251‑1 (MQ=255)                                                                                                                                                                                                                           taggacaggaggctatGGGTAAAGCATTTG‑TCATAGTTGAGTCCCCGGCAAGAGCCAAGACGCTCAAAAAGCATCTGGGTAGTGACTACGTGGTGAAATCCCGCGTCGGTCACATCCTAGATTTGCCGACCCGCGGCTAATCTGCCAAAAAGAGAGCCGACTCTACCGCCACCACGACGGCTAAACAGCCTAAAGAGGATGAACGTAGCGCCCTCGTCAACAGTATGGGTGATTACCCGTG‑CACaacttgg             >  2:74081/14‑246 (MQ=255)                                                                                                                                                                                                                            aGGTAAGGTGAATATGGTTAAAGCATGTG‑TCATCGTTGAGACCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGGCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGATTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCGCACGTCAACCGTATGGGGGTTGACCCGTGGCACAATTggg             <  2:159280/251‑1 (MQ=255)                                                                                                                                                                                                                            aGGTAAGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTggg             <  1:160111/251‑1 (MQ=255)                                                                                                                                                                                                                             ggTAAGGTGAATCTGGGTAAAGCATTTG‑TCATCGTTGGGTCCCCGGCAAAAGCAAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAACGAGTGCCGACTCTACCTCCCCCAAGACGGCTAAAAAGCCTAACAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCCCCATTgggg            >  2:465695/1‑250 (MQ=255)                                                                                                                                                                                                                                 aGGTGAATATGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAACAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGGTGACCCGTGGCACAATTGGGAGgcg        >  2:457152/1‑251 (MQ=255)                                                                                                                                                                                                                                       atatGGGTAAAGCATTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGCCCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAGGACGGCTAAAAAGCCTAACAAGGAGGAACGTGACGCTCTAGTCAACCGTATGGGGGTGGACCCGGGGCACAATTGGGAGGCGCACTAt  >  2:354359/1‑251 (MQ=255)                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                              ACGCTTTCGTCAATCTGGCAATAGATTTGCTTGACATTCGACCAAAATTCCGTCGTGCTATAGCGCCTGTAGGCCAAGACCTGTTAACTCAGTCACCTGAATTTTCGTGAACAGAGTCACGACAAGGGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTG‑TCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCAGCGTCGGTCACATCCGCGATTTGCCGACCAGTGGCTCAGCTGCCAAAAAGAGTGCCGACTCTACCTCCACCAAGACGGCTAAAAAGCCTAAAAAGGATGAACGTGGCGCTCTCGTCAACCGTATGGGGGTTGACCCGTGGCACAATTGGGAGGCGCACTAT  >  NC_000913/1330816‑1331296

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