Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,223	T→C	N207N (AAT→AAC)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,223	0	T	C	95.6%	190.5 / ‑0.0	70	N207N (AAT→AAC)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (30/35);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

AAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAG  >  NC_000913/4619979‑4620464                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     aaggtttgctgGATAACCGGGCCGACCGCGCTAATGTGCCGG‑TGTCCCCCGTAACCTGCCACTCTTCCGTTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACGCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGt                                                                                                                                                                                                                                                <  2:21994/245‑1 (MQ=255)         gtgGATAAACGGGCCGAACGCGTCAATCTGCCGGTTTGCCTCGGAAAC‑TGACACTCTTCCGGTACGGTCATTCTGGATCAGCGGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TAt                                                                                                                                                                                                                                       <  2:508139/250‑1 (MQ=255)            gATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCCAAAACATCGGTCGGGT‑TATCGc                                                                                                                                                                                                                                    >  2:294656/1‑251 (MQ=255)                  cgggTCGAACGCGCTAATCGGCCGGGTTACCTCGGTAAC‑TGCAACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTcc                                                                                                                                                                                                                              <  1:508257/249‑1 (MQ=255)                   tGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCg                                                                                                                                                                                                                             >  1:564585/1‑251 (MQ=255)                       caaccgcgcagttcttccGGGTTCCCCCTGTAAC‑GGCCCCTCTTCCGGGACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCGATACCTCCGCTGCCCCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTGTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCTTCCGTTTa                                                                                                                                                                                                                         <  2:350011/235‑1 (MQ=255)                                tAATCTACCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGCGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAACCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCCACACAATCGGTCGGGTGTA‑CGCCCATCCGTTTATCGGTGAc                                                                                                                                                                                                                 >  2:587606/1‑250 (MQ=255)                                 aaTCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGT‑TATAGCTCGTCCGTTTATCGGTGACa                                                                                                                                                                                                                >  1:78078/1‑250 (MQ=255)                                  atcggCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATGCGGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACaaa                                                                                                                                                                                                              <  1:607517/247‑1 (MQ=255)                                          ggTTACCTCAGTAAC‑GGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTc                                                                                                                                                                                                                                  <  1:27740/223‑1 (MQ=255)                                                      aaccggCCACTCTGCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTa                                                                                                                                                                                           <  2:316872/246‑1 (MQ=255)                                                       aC‑TGCCACTCTTCCGGTACGGTCATTCTGGAACAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTGCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCTGTAACTTCCAGCGTAcc                                                                                                                                                                                         <  2:526310/251‑1 (MQ=255)                                                           gCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCgg                                                                                                                                                                                       <  1:319516/250‑1 (MQ=255)                                                            cccgacgctcttccgaactaTTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                                                                                                                                                                                     <  2:29512/232‑1 (MQ=255)                                                             cACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                                                                                                                                                                                     >  1:301911/1‑250 (MQ=255)                                                                       gTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGGCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGAc                                                                                                                                                                          <  1:26666/251‑1 (MQ=255)                                                                         aCGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCt                                                                                                                                                                        <  1:301163/251‑1 (MQ=255)                                                                           ggTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAACTCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACAACCGTCGTGC‑TATCGCTCGTCCGGTTATCGGTGACCAAGCCGGTAACATCCCGCGTACCGGTCACCGTCACGAACTgg                                                                                                                                                                      >  2:251500/1‑251 (MQ=255)                                                                              cATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTg                                                                                                                                                                   <  1:130413/251‑1 (MQ=255)                                                                               aTTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTa                                                                                                                                                                                     >  1:29512/1‑232 (MQ=255)                                                                                  cTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTgg                                                                                                                                                              >  1:94773/1‑250 (MQ=255)                                                                                      aTCAGCTGGGCGAAGAGCACATGAAAAACGGCAAGCCGATTTTCTATACCTCCGCTGACTACGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAAATCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCCGTCGGTT‑TATCGCTCGTCCCTTATTCGGTGACAAAACCGGTCACTTCCCACGTACCGGTAACCGtaccgaccgggaagtggcgacg                                                                                                                                                           >  2:675156/1‑230 (MQ=255)                                                                                             gggCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAACCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCGG‑TTGAGCCGCCCGCac                                                                                                                                                   >  1:465348/1‑251 (MQ=255)                                                                                                  aaGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATAGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCaccgac‑c                                                                                                                                             >  1:478383/1‑251 (MQ=255)                                                                                                      gCACATGAAAACCGGCAAGACGATTTTCTATACCTCCGCTGACTCCGTGATCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGGGCGAAATCGCCCGAGAAGAGCTGACAAACGGCGGCTACAAACTCGGTCGGGG‑TAACGCACGAACGTTTAACGGGGACCAAGCAGGTAAATTCCAGCGTACCGGTAAACGTGACGACCTGGCTG‑GGGGGCCTCCAGCaccgac‑cgggc                                                                                                                                         >  2:9619/1‑248 (MQ=255)                                                                                                       cacaTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGtgct                                                                                                                                        >  1:539849/1‑251 (MQ=255)                                                                                                         caTGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAAACGTCACGACCTGGCTG‑TTGGGCCGCCAGCCACGAC‑CGtgctgg                                                                                                                                      >  1:133902/1‑250 (MQ=255)                                                                                                          aTGAAACCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCa                                                                                                                                     <  1:57819/251‑1 (MQ=255)                                                                                                             aaaaCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGGGACAAAGACGGTAACTTCCAACGTACCGGGAACCGTCACGACCGGGCTG‑GTGGGCCCCCAGGacccgccctggtgcaga                                                                                                                                    >  2:558897/1‑236 (MQ=255)                                                                                                                aCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTTACAAAGCCGGTAACTTCCAGCGTACCGGGAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACt                                                                                                                               >  1:537280/1‑251 (MQ=255)                                                                                                                aCCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACt                                                                                                                               <  1:651832/251‑1 (MQ=255)                                                                                                                aCCGGAAAGCCTATTTTCGATCCCTCTGCTGATTCCGTGTTCCAGATTGCCTTCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTTCGAAATCGCCCTTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCGCGTCCGTTTACCGGTGACAAAGCCGGTAACTTCCAGCGTTCCGTTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACt                                                                                                                               <  2:106666/251‑1 (MQ=255)                                                                                                                 ccGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGGCACGCCCTGGCTG‑TTGGGCCGCCAGCACCGAA‑CGGGCTGCAGAAACTg                                                                                                                              >  1:557871/1‑251 (MQ=255)                                                                                                                    gCAAGCCGATTTTCTATTCCTCCGCGGACGCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGAAAAAGCCGTTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGtt                                                                                                                           <  1:47343/251‑1 (MQ=255)                                                                                                                          cGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTTTCGGGGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCCGCCCCCAC‑CGCGCTGCAGAAACTGGTtgatta                                                                                                                      >  2:235200/1‑248 (MQ=255)                                                                                                                             ttttCTATACCTCCGCTGACTCCGTGTTCCCGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGaaaa                                                                                                                   >  1:520621/1‑250 (MQ=255)                                                                                                                               ttCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAAATGGTTGATGAAAAac                                                                                                                 >  1:90776/1‑250 (MQ=255)                                                                                                                                         ccGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTGCAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtg                                                                                                      <  1:235200/251‑1 (MQ=255)                                                                                                                                          cGCTGACTCCGTGTTCCAGTTTGCTGGACACGAAGGAACCGTTGGCCTGGATAAACACTACGAAATGTGCGAAATCGCCCGTGAAGAGCGGCCCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTGGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAggtgg                                                                                                     <  2:172989/251‑1 (MQ=255)                                                                                                                                            cgtactccgtttgtcAGATTGCCTGCCCAGAAGGAACCTTTGGCCTGGATAAACTCTACGTACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGtt                                                                                                   <  2:564585/237‑1 (MQ=255)                                                                                                                                             tGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGA‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGACACTTCCAGCGTACCGGGAAACGGCACGACCTGGCTGGTTG‑GCCGCCAGCCCCGACACGGCCTAC‑GAAACCTGTTGATGAAAAACACGGGCAggagggtt                                                                                                  >  2:332042/1‑245 (MQ=255)                                                                                                                                                cTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAac                                                                                                                 <  1:235364/233‑1 (MQ=255)                                                                                                                                                cTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGAGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGAACTGGCTG‑TTGGGCCGCCAGCACCGAC‑CGGGCTGCAGAAACTGGGTGGTGTAAAacagaacggaagagaggcgg                                                                                               >  2:235364/1‑234 (MQ=255)                                                                                                                                                  ccGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGTTTCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTc                                                                                             <  1:587606/251‑1 (MQ=255)                                                                                                                                                       ttCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGGGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGGGCCGCCAGCAACGAC‑CGTGCTGCAGAAACCGGTTGATGAAAAAAACGGTAAGGGGGGTTCTGTCGGTa                                                                                         >  1:371230/1‑250 (MQ=255)                                                                                                                                                         ccAGATGGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTTTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTaaa                                                                                       <  2:94773/250‑1 (MQ=255)                                                                                                                                                           aGATTGCCTGCCACGAAGGAACCTTGGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACGTCCAGCGTACCGTTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAAtt                                                                                    <  1:319198/251‑1 (MQ=255)                                                                                                                                                              ttGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTtgatga                                                                                                                      >  1:360305/1‑214 (MQ=255)                                                                                                                                                                gcctgccACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGcaga                                                                               <  1:626593/251‑4 (MQ=255)                                                                                                                                                                              aacctttggccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGTCGGCTACAACATCGGTCAAGT‑TATCGCTCGACCGTTTATCGGAGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGATGCATAAACTGGTTGATGAAAAAAAAGGACAGGGGGTTTCTGTCGGTAAAATTGAAGACATATACGCCAACt                                                                 >  1:397843/11‑251 (MQ=255)                                                                                                                                                                                  tttggccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGATACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCAAGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGGGCTGCCGAAACTGGTTGATGAAAAACACGGTCAGGGGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCACACTGCgg                                                             >  1:208218/7‑251 (MQ=255)                                                                                                                                                                                    tggccTGGATAAACCCTACGAACTGGGCGAAATCGCCCGTGTAGAGCTGACCAACGGCGGCTACAACAGCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGt                                                            <  2:345062/246‑1 (MQ=255)                                                                                                                                                                                      gccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGACGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCACCACCGCC‑CGTGCTGCAGAAACTGGTTGATGGAAAACACGGTCAGGTGGTTTCTGTAGGTAAAATTTCAGACATCTAAGACAACTCCGGtgtc                                                         >  1:105480/3‑248 (MQ=255)                                                                                                                                                                                       ccTGGATTAACTCTCCGAACTGTGCGAAATCGCGCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGATAAAGCCGGTAATGTCCAGCGTACCGGTAACCGTCACGACCTGGTTG‑TTGAGCCGCCATCACCGAC‑CGTGCTTCAGAAACTGGTTGATGAAAAACACTGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCa                                                        <  1:617375/250‑1 (MQ=255)                                                                                                                                                                                       ccTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCa                                                        <  1:446750/250‑1 (MQ=255)                                                                                                                                                                                           gtaaaaccctccGAACTGTTCGAAACTGCCAGTGAAGTGCTGCCCAACGGCGGCTACAACTTCGTCCGTTT‑GATCGCTCGTCCGTTTATCGGTTCCAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACAACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaa                                                    <  2:172842/240‑1 (MQ=255)                                                                                                                                                                                           gATAAACTCTACGAACTGGGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaa                                                    <  1:501122/251‑1 (MQ=255)                                                                                                                                                                                              aaaCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCaaga                                                  >  1:513447/1‑248 (MQ=255)                                                                                                                                                                                                            gtgCGAAATCGCCCGTGAAGAGCTGTCCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTg                                   <  2:363562/251‑1 (MQ=255)                                                                                                                                                                                                            gtgCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACt                                    <  1:231169/250‑1 (MQ=255)                                                                                                                                                                                                                     cGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGGGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGAc                           <  1:168044/250‑1 (MQ=255)                                                                                                                                                                                                                      gCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGAcgc                         >  1:156971/1‑251 (MQ=255)                                                                                                                                                                                                                           tGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTaa                                                                                        >  1:235136/1‑183 (MQ=255)                                                                                                                                                                                                                           ggAAGAGCTGACCAACGGCGGCTACAACATCTGTCGTGT‑TATCGCTCGTCCGTTTACCGGTGACAAAGCCGGTACCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTaa                                                                                        <  2:235136/182‑1 (MQ=255)                                                                                                                                                                                                                             aaGAGCTGACCAACGGCGGCTACAACATCGGTCGGGT‑TATCGCTCGTCTGTGTATCGGTGAAAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGGTGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGttt                  <  2:388344/251‑1 (MQ=255)                                                                                                                                                                                                                                   tGACCAACGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGcc            >  1:386476/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaCGGCGGCTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCAt       <  1:294656/251‑1 (MQ=255)                                                                                                                                                                                                                                             cggcTACAACATCGGCCGTGT‑TATCGCTCGCCCGTTTATCGGCGACAAAGCCGGTAAATTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCCCCGAC‑AGCGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGCGGTTTCTGTCGGTAAAATTGCAGACCTCTACGCCAACTGCGGTACCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGAAGCCACCAGCaa    >  2:293571/1‑249 (MQ=255)                                                                                                                                                                                                                                               gcTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAag  >  1:315824/1‑249 (MQ=255)                                                                                                                                                                                                                                                cTACAACATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGTGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAag  >  1:34968/1‑248 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     AAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAG  >  NC_000913/4619979‑4620464

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: ‑