Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,846,553	G→A	D136D (GAC→GAT)	topB ←	DNA topoisomerase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,846,553	0	G	A	96.7%	183.7 / ‑1.6	61	D136D (GAC→GAT)	topB	DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (32/27);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

CCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGC‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGCAAGATTCCAGCGCGCATAGCGGCT  >  NC_000913/1846309‑1846795                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   ccACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGC‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGttt                                                                                                                                                                                                                                              >  2:383877/1‑251 (MQ=255)       aGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGC‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCa                                                                                                                                                                                                                                          >  2:39111/1‑249 (MQ=255)       aGTTTTCAATCCCTTCATCGCGGCGCACCACCAGCCCAAGC‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCa                                                                                                                                                                                                                                          <  1:39111/249‑1 (MQ=255)           ttCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGC‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAAc                                                                                                                                                                                                                                    >  1:322594/1‑251 (MQ=255)                    ttCATCGCGGCGCACCACCAGCACAAGC‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACct                                                                                                                                                                                                                           >  1:198475/1‑251 (MQ=255)                      cATCGAGGCGCACCACCAGCCCAAGC‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACctgc                                                                                                                                                                                                                         >  2:373577/1‑251 (MQ=255)                               gcACCACCAGCCCAAGC‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCttt                                                                                                                                                                                                                >  1:455718/1‑251 (MQ=255)                                caccacCAGCCCAAGC‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCtttt                                                                                                                                                                                                               >  2:343836/1‑251 (MQ=255)                                   cagccgCCCAAGCCACGGGCGTCTGCACGCGTCCCACGGAAAGTCCGCCCGGATACCCGATCTTGCGACCGAGATGGGTATACGCACGGGTCATATTGATGCCTTACAGCCAGCCGGCACGCGCTCGCGCCAGGGCAGAACCGCACAGCGGTACAAACTCCCGGTTGGAACGAATACGGTCGATCGCCCGCTAAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTcc                                                                                                                                                                                                             <  1:255633/246‑1 (MQ=255)                                        cccaaaGC‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGTGCGAGAATGGTATACGCACGGGTCATATTGGTGCCGTACAGCCAGTCTGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCACCCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGcc                                                                                                                                                                                                       <  1:194030/247‑1 (MQ=255)                                         cccAAGC‑ACGGGCTTCTGCACGCTCCCCACGGAAAGTACGCCCTGATACCCGGCAGTGCGACCGAGAATGGTATCCGCACGGGTCATATTGATGCCTTACAGCCAGTCGGCACGCGCTCGCCCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGAAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCa                                                                                                                                                                                                      <  2:292347/251‑1 (MQ=255)                                          cccacgaaagtgCGTCTGCACGCCCCCCCCGAAAAGTCCCCCCTGATAACCGGCTTTCCAACCGAGAAGTGTTTACGCACTTGTCATATTCATTCCGTCCCGCCGGTCGCCACGCTCCCGCGCCGGCGCAGAATCGCCCAGCGGTTCAAACTTTCTGTTGCACCTAAAACGGTCAAGCCCCCTCTCACCCGCCTGCGGGTTCAGATCGTTTATCAACCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCa                                                                                                                                                                                                      <  1:254808/240‑1 (MQ=255)                                               c‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGACCTGATAACAGGCATTGCGACAGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGACGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGc                                                                                                                                                                                                 >  1:272449/1‑250 (MQ=255)                                                    ggCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATACCCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATa                                                                                                                                                                                            <  1:131462/251‑1 (MQ=255)                                                          tGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTcc                                                                                                                                                                                        >  1:325651/1‑249 (MQ=255)                                                                   cccACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCCCTTCATc                                                                                                                                                                             >  2:87314/1‑251 (MQ=255)                                                                     cACGGAAAGTACGCCCTGATACCCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATcca                                                                                                                                                                           <  1:105214/251‑1 (MQ=255)                                                                        ggAAAGCACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCAAGTCGGCTCGCCCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAGGACGGTAGATCGTCCGCTCAACCGCATGCGGGTTCAGATCGTTTATCAAGAAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATatgctagaacttcatacacc                                                                                                                                                                         >  2:100814/1‑231 (MQ=255)                                                                                     tGATAACCGCCATTGCGACCGAAAATGGTATACGCACGGGTCATATTGATGCCGCACACCCAGTCGGCACGCGCTCGCGCCAGCGCAAAAACGCACGGCGGTCCAAACTCACTTTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCtt                                                                                                                                                            <  1:235012/250‑1 (MQ=255)                                                                                      gATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCtt                                                                                                                                                            <  1:170680/249‑1 (MQ=255)                                                                                                    cGACCGAGAATGGTATACGCACGGGTCATTTTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTcccc                                                                                                                                            <  1:451559/251‑1 (MQ=255)                                                                                                    cGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAACCGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTcccc                                                                                                                                            <  1:343836/251‑1 (MQ=255)                                                                                                     gACCGAGAATGTTATACGCCCGGGTCATATTGTTGCCGTACAGCCAGTCGGCACGCGCTCGCGCAAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCg                                                                                                                                           <  2:147543/251‑1 (MQ=255)                                                                                                      cccGAAAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCgg                                                                                                                                          <  1:413789/250‑1 (MQ=255)                                                                                                        cGAGAATGGTATACGCACGGGTCCTATTGATGCCGTACAGCCAGTCGGCACGAGCTCGCGCCAGCGCAGAAACGCACAGAGGTGCAAACTCACTGTTGGAACGAAGACGGTCGATCGCACGCTCAACCGCCTGCGGGATCAGATCGTTTATCAAACAACGCTGTACCTGCTGGCGCTTTTACGGTGGCAGTTGCAGCTAGTCCAGCAATTCATCCACCATCAATTGCCCCTCACGATCCGGATACCCGGCg                                                                                                                                        >  2:343681/1‑251 (MQ=255)                                                                                                             aTGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAAc                                                                                                                                   >  2:282714/1‑251 (MQ=255)                                                                                                                    aCGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCg                                                                                                                            >  1:27058/1‑251 (MQ=255)                                                                                                                      gcccgggtaaaatTGAGGCCGTACAGCCAGTCGGCCCGCGTTCGCCCCAGCGCAGAAACGCACAGCGGTACAAACTCCCTGTTGGAACGAAGACGGTCGATCGCCCGCTCACCCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTc                                                                                                                             <  2:32582/237‑1 (MQ=255)                                                                                                                      gCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTc                                                                                                                             >  1:32582/1‑248 (MQ=255)                                                                                                                              cATATTGAGGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGACCGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCAt                                                                                                                  <  1:423351/251‑1 (MQ=255)                                                                                                                              aatgttggagCCGTACCGCCGGGCGGCTCGCGCTCGCGCCCGCGCCGAACCGCACAGCGTTCCAACCTCACTTTTGGAACGAAGCCGTTCGTCCGCCAGCGCATCCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCAt                                                                                                                  <  2:9831/242‑1 (MQ=255)                                                                                                                                tatTGATGCCGTACAGCCAGTCGCCACGCCCTCGCGCCAGAACAGAAACGCACAGCGGTCCAAACTCACCGTTGGAACGCAGACGGTCGCTCGCCCGCTCAACCCCCTGCGGGTTCAGATCGTTTATCACGCAACGCTGTACCTGCTCGCGCTTTTCCGGTGCCAGTTGCAGACAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTACCCGGCGTGAACGATTTCGCTGGCTTCATGc                                                                                                                >  2:28850/1‑251 (MQ=255)                                                                                                                                   tGATTCCGCACAGCCAGTCGGCCCGCGCTCGCGCCAGCCAAGAAACGCACAGCGTTACAAACTCTCTGTTGGAACGAAGACGGGCTATCGCCCACTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAg                                                                                                              <  2:429370/250‑1 (MQ=255)                                                                                                                                           tACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGCACCGttt                                                                                                     >  1:261079/1‑251 (MQ=255)                                                                                                                                                ccAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTCCAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTtgatga                                                                                                <  1:106931/251‑1 (MQ=255)                                                                                                                                                      ggtaCGCCCTCGCGCCAGCGCAGCAAAGCACAGCGGGACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCCCAAACGCCTGCGGGTACAGATCGTTTATCAAGCAACGCTGTACATGATGGCGCTTTTCCGGTGCCGGTTGCAGATAGTTCAGCACTTCATGCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGAAGGCACCGTTTGCTGACGTTa                                                                                           >  2:35305/4‑250 (MQ=255)                                                                                                                                                      ggCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTa                                                                                           >  2:120152/1‑250 (MQ=255)                                                                                                                                                      ggCACGCGCTCGCGCCAGCGCAGAAACGAACAGCGGTACAAACTCACTGTTGGAACGAAGACGGGCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGTGTCCCCGGCGTGATCGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTa                                                                                           >  1:429102/1‑250 (MQ=255)                                                                                                                                                      ggCACGCGCTCGCGCCAGCGCAGAAAAACACAGCGGTACAAACCCACTGTTTCAACGACGACGGTCGATCGCCCGCTCAACCACCTGCGGGTTCAGAACGCTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGCGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATACGGGTCCCCGGCGAGAACGCTTTCGCTGGCTTCATGCAGGAACAGTTTGGTGAccntc                                                                                           >  2:410888/1‑246 (MQ=255)                                                                                                                                                         gcgcgCTCGCGCCAGCGCAAAAACTCACAGCGGTACAAACTCGCTGTGGGAACGAACCCGTTCGATCGCCCGCTCAACCGCCTGCGTGTTCAGGTCGTTTATCAAGCAACGCTGGACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCGCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGGTAAgtt                                                                                       <  2:261079/250‑1 (MQ=255)                                                                                                                                                         aCGCGCCCGCGCCAGCGCAGAAGCGCACAGCGGTACAAACTCACTGTTGGACCTAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAgtt                                                                                       <  1:373577/251‑1 (MQ=255)                                                                                                                                                              cTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCCGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGAGCCAGTTGCAGATAGTC‑AGCACTTCATCCACCAGCAATTGCCCTTCACGAACCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAATCGTTTGATGACTTTAAGTTGTTTgg                                                                                 >  1:232113/1‑251 (MQ=255)                                                                                                                                                                   cccagagaaGAAACGCACGGCGGTCCAATCTCTCTTTTGGAACTAAGCCGGTCGATCCCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCAc                                                                             <  2:326382/243‑1 (MQ=255)                                                                                                                                                                        ccccaaaaCGCACAGCGGTACATACACCCTATTGGAACAAATACGGTCGATCGCCCGCTCAACCGCCCTCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCTGTGCCAGTTGCAGATAGTC‑AGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAggg                                                                       <  2:232113/246‑1 (MQ=255)                                                                                                                                                                              aaCGCACAGCGGTACAACCCCACTGTTGGAATGAAGACGGTCGCTCGCCCGCTTATCCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTTCCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCggg                                                                  <  2:327543/251‑1 (MQ=255)                                                                                                                                                                                   acaGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGt                                                             <  1:320095/251‑1 (MQ=255)                                                                                                                                                                                            cAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCAc                                                    >  2:169370/1‑251 (MQ=255)                                                                                                                                                                                                cTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGACCCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACtttt                                                >  2:340376/1‑251 (MQ=255)                                                                                                                                                                                                             cGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCgg                                   <  1:120152/251‑1 (MQ=255)                                                                                                                                                                                                               aagacggacgctagtgctctgctccGCCTGCGGGTTCAGATCGATTATCAAGCAAAGCTGTACCTGCTGGAGCTTTTCCGGTGCAAGTTGCTGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCAAGATCAGGGTCCCCCGCGAGAACGATTTAGCTGGCTGCATGCAGGAACCGGTTGGTGAAGTTAAGCTGTTTGGTCACGGAGGGACGGGGCTGTAATTGCCGCTTTTCCGGGACAATCGGcc                                 >  2:37829/24‑250 (MQ=255)                                                                                                                                                                                                                    ggTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATcc                            >  1:91084/1‑251 (MQ=255)                                                                                                                                                                                                                      tCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGc                          >  2:354954/1‑251 (MQ=255)                                                                                                                                                                                                                      tCGATCGCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGc                          >  2:64509/1‑251 (MQ=255)                                                                                                                                                                                                                        gATCGCCTGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAACCAACGCTGTACCTGCTGGCGCTTTTCCGGAGCCAGTTGCACATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGCCCCCGGCGTGAACGACTTCGCTGGCTTCATGCAGGAGCCGCTTGATGAAGGTAAGATGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGCa                         >  2:46445/1‑250 (MQ=255)                                                                                                                                                                                                                            gCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGCAAGAtt                    >  1:365603/1‑251 (MQ=255)                                                                                                                                                                                                                            gCCCGCTCAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGCAAGAtt                    >  2:132393/1‑251 (MQ=255)                                                                                                                                                                                                                                  ccacccGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGCAAGATTCCAgcg              <  2:455718/247‑1 (MQ=255)                                                                                                                                                                                                                                  ccAACCGCCTGCGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAAACGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGCAAGATTCCAgcg              <  1:330131/250‑1 (MQ=255)                                                                                                                                                                                                                                       cGCCTGCGGGTTCAGATCGTTTATCAATCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGCAAGATTCCAGCGCGCAt         <  1:236276/251‑1 (MQ=255)                                                                                                                                                                                                                                             cGGGTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGCAAGATTCCAGCGCGCATAGCGGc   >  2:154572/1‑251 (MQ=255)                                                                                                                                                                                                                                              gggTTCAGATCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGCAAGATTCCAGCGCGCATAGCGGCt  >  1:74839/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   CCACGAAGTTTTCAATCTCTTCATCGCGGCGCACCACCAGCCCAAGC‑ACGGGCGTCTGCACGCGTCCCACGGAAAGTACGCCCTGATAACCGGCATTGCGACCGAGAATGGTATACGCACGGGTCATATTGATGCCGTACAGCCAGTCGGCACGCGCTCGCGCCAGCGCAGAAACGCACAGCGGTACAAACTCACTGTTGGAACGAAGACGGTCGATCGCCCGCTCAACCGCCTGCGGGTTCAGGTCGTTTATCAAGCAACGCTGTACCTGCTGGCGCTTTTCCGGTGCCAGTTGCAGATAGTCCAGCACTTCATCCACCAGCAATTGCCCTTCACGATCCGGGTCCCCGGCGTGAACGATTTCGCTGGCTTCATGCAGGAACCGTTTGATGACGTTAAGTTGTTTGGTCACGGAGGGTCGGGGCTGTAATTGCCACTTTTCCGGGACAATCGGCAAATCCGCAAGATTCCAGCGCGCATAGCGGCT  >  NC_000913/1846309‑1846795

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