Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,331,947	C→T	G300G (GGC→GGT)	topA →	DNA topoisomerase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,331,947	0	C	T	97.6%	136.7 / ‑7.0	42	G300G (GGC→GGT)	topA	DNA topoisomerase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/22);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

CTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGTCAATGTGATGGCGGAATCGCTGAA  >  NC_000913/1331704‑1332189                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   cTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGaaa                                                                                                                                                                                                                                              >  1:102407/1‑250 (MQ=255)   gggAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGGACCCGTCTTGGATGTGGTGTGaaaa                                                                                                                                                                                                                                             <  2:54344/249‑1 (MQ=255)    ggAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGaaaaaaa                                                                                                                                                                                                                                          <  1:157810/251‑1 (MQ=255)              gCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCatgatga                                                                                                                                                                                                                                 <  1:272662/250‑1 (MQ=255)               ccAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGGGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTTCGCGGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCatgatgatg                                                                                                                                                                                                                               <  1:209448/251‑1 (MQ=255)                    aCGACCACGCCAGCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTGCCGTCCGGGCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAGCCTGGCGCTCCTTGTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTGTTGGATTTGGTGTGAAAAAAGCCATGAGGATGGCGCa                                                                                                                                                                                                                          <  2:15363/251‑1 (MQ=255)                           cGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTg                                                                                                                                                                                                                   >  2:257101/1‑251 (MQ=255)                            gCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGt                                                                                                                                                                                                                  <  1:209129/251‑1 (MQ=255)                            gCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGGCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTTCCTCTACGCTGCAACAAGCTGCCAGCAGCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGt                                                                                                                                                                                                                  <  2:139444/251‑1 (MQ=255)                               aTCTGGTGAAGCGTTGGCGTTACAGGTGTCTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATg                                                                                                                                                                                                               >  1:89583/1‑251 (MQ=255)                                   ggTGAAGCGGTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATg                                                                                                                                                                                                               <  2:108361/247‑1 (MQ=255)                                           gTTGGCGTTCCAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCCGGCTGCGGAAACTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCAGGCGCTCCATTTATTACCTCTACCCTGCAACCAGCTGCCAGCACCCGTCTTGGCTTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGATTGTATGAAGCAGGCtat                                                                                                                                                                                                    >  2:323585/1‑250 (MQ=255)                                                    aCAGGTGACTCATCAGAACGAAAAAGCGTTCCGTCCGTTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTAAAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCAGCAACAAGCTGCCAGCACCCGTCTTGGCTTTGGTGTGAACAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTAcc                                                                                                                                                                                          >  2:153251/1‑250 (MQ=255)                                                      aGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATAGCACTTACATg                                                                                                                                                                                        <  2:238333/251‑1 (MQ=255)                                                                           aaCCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGaaa                                                                                                                                                                                                                                              <  1:263534/176‑1 (MQ=255)                                                                           aaCCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGaaa                                                                                                                                                                                                                                              >  2:263534/1‑176 (MQ=255)                                                                           aaCCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCGGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTg                                                                                                                                                                   <  2:119968/251‑1 (MQ=255)                                                                              cGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAGCGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGt                                                                                                                                                                >  2:185003/1‑251 (MQ=255)                                                                                 tccgtccgGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCa                                                                                                                                                              <  1:94932/250‑1 (MQ=255)                                                                                       cgGTCAATAAACAAAAAAATCCGGTTGCGGTAAGTCTGCTGGAAAAAGCGCGCCACAGCGTGCTGGAACGTGAAGCCAAACCGACAACCAGTAAACCTGGCGCTCCTCTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGATATATCACTTACATGCGTATCGACTCCACTAACCTGAGTCAGGAcgcg                                                                                                                                                       >  2:149808/1‑251 (MQ=255)                                                                                       cgGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGAcgcg                                                                                                                                                       >  2:149826/1‑251 (MQ=255)                                                                                        ggTCAACAAAGCACAAACCGAGGATGCGGTCCGTCTGCTGGAAAAAGCGCGCTACAGCCCCCTGGAACGTGAAGACAAACAGACAACCAGTAAACCAGCCGCTCCTTTTATTACCTCTACGCTGCAACAAGCAGACAGCACCCGTCTTGGATTTGGTGTGCAAAAAACCATGATGATGGCGAAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGAcgagg                                                                                                                                                      >  2:225914/1‑248 (MQ=255)                                                                                                     cAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTcg                                                                                                                                         >  2:116922/1‑251 (MQ=255)                                                                                                              gCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTatat                                                                                                                                >  2:316468/1‑251 (MQ=255)                                                                                                                   ggTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCAGTTATATCAGCg                                                                                                                           <  1:185003/251‑1 (MQ=255)                                                                                                                                                        gAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGa                                                                                       >  2:168563/1‑250 (MQ=255)                                                                                                                                                                                                  tgtattaCCTCTACGCTGCAACAAGATGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATAGCACATACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGaa                                              <  2:89583/247‑1 (MQ=255)                                                                                                                                                                                                                   gCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGAc                            <  2:286100/250‑1 (MQ=255)                                                                                                                                                                                                                   gCAACAAGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGGCCTTCTGACg                           <  2:9402/251‑1 (MQ=255)                                                                                                                                                                                                                    caacaaGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAACAACCAATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAAACTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATAAGTACGCCAGAAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGt                          >  1:225411/1‑251 (MQ=255)                                                                                                                                                                                                                       caaGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGCTTc                                       >  2:328278/1‑235 (MQ=255)                                                                                                                                                                                                                       caaGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTc                                       <  1:328278/235‑1 (MQ=255)                                                                                                                                                                                                                         aGCTGCCAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGTCAAtg                     >  1:288656/1‑251 (MQ=255)                                                                                                                                                                                                                               cAGCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGTCAATGTGATgg               >  1:182448/1‑251 (MQ=255)                                                                                                                                                                                                                                 gCACCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGTCAATGTGATGGCg             <  1:26025/251‑1 (MQ=255)                                                                                                                                                                                                                                   aCCCGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACGAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGAGTCGTCCTTCTGACGTCAATGTGATGGCGGa           <  2:323631/251‑1 (MQ=255)                                                                                                                                                                                                                                    cccGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGTCAATGTGATGGCGGa           >  2:87238/1‑250 (MQ=255)                                                                                                                                                                                                                                     ccgttgttgATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCTTTTGTATGAAGCAGGCTATATCGCTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGTCAATGTGATGGCGGAAt         <  2:219686/243‑1 (MQ=255)                                                                                                                                                                                                                                      cGTCTTGGATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGTCAATGTGATGGCGGAATc        <  2:306756/251‑1 (MQ=255)                                                                                                                                                                                                                                            ggATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGTCAATGTGATGGCGGAATCGCTGaa  <  1:46504/251‑1 (MQ=255)                                                                                                                                                                                                                                             gATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATCTGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGTCAATGTGATGGCGGAATCGCTGa   >  2:349217/1‑249 (MQ=255)                                                                                                                                                                                                                                             gATTTGGTGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGTCAATGTGATGGCGGAATCGCTGa   >  1:284459/1‑249 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   CTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGCGATTCGTCCTTCTGACGTCAATGTGATGGCGGAATCGCTGAA  >  NC_000913/1331704‑1332189

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