Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,620,145 T→C H181H (CAT→CAC)  deoB → phosphopentomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,620,1450TC95.8% 209.6 / ‑0.5 72H181H (CAT→CAC) deoBphosphopentomutase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (39/29);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01

GGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGA  >  NC_000913/4619903‑4620387
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
ggccACTGGGAAATCGCTGGTGTCCCGTTTTTGTTGGAGTGTGGATATTTCGCCGATCACTAAAACGGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCGCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACgaagg                                                                                                                                                                                                                                                  <  2:236002/245‑2 (MQ=255)
                          ggTTCTGTTTGAGTGGGTATAGTTCTCCGATCACGAAAACAGCTTCCCGCAAGTGCTGCTGGATAAACTGTTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTCCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAAc                                                                                                                                                                                                                     <  1:61272/251‑1 (MQ=255)
                                      ggggggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCGGCCGGGTTACCTGGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCg                                                                                                                                                                                                         <  1:606614/249‑1 (MQ=255)
                                        ggggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACgaag                                                                                                                                                                                                                                                   >  1:394962/1‑207 (MQ=255)
                                           gATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGCTCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTAGATAAACTCTACGAACTGTGAGAAATCGCCCgg                                                                                                                                                                                                    >  2:169618/1‑250 (MQ=255)
                                           gATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGt                                                                                                                                                                                                    >  1:363562/1‑251 (MQ=255)
                                               tttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGCTTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGa                                                                                                                                                                                                  >  2:319516/1‑249 (MQ=255)
                                                         cACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCa                                                                                                                                                                                       >  1:172842/1‑250 (MQ=255)
                                                          aCGAAACCAGTTTCCCCCAAGAGCTGCGGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAc                                                                                                                                                                                     <  1:25639/251‑1 (MQ=255)
                                                          aCGAAAACAGCTTCCCGCAAGGGCTGGTGGATAAACTGTTCGATCGCGCTAATCTGCCGGTTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGTTCATTCTGGATCATCTGGGCGAAGATCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTTGATAAACTCTACGATCTGTGCGAAATCGCCCGTGAAGAGCTGACCAAc                                                                                                                                                                                     <  1:102605/251‑1 (MQ=255)
                                                            gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcgg                                                                                                                                                                                   >  1:21994/1‑251 (MQ=255)
                                                            gAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcgg                                                                                                                                                                                   >  1:224817/1‑251 (MQ=255)
                                                             aaaaCAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAACCCGGCAAGCCGATTTTCTCTACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGCGCTGACCAAaggc                                                                                                                                                                                  >  1:50065/1‑251 (MQ=255)
                                                                  aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCt                                                                                                                                                                              >  1:46657/1‑250 (MQ=255)
                                                                  aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCt                                                                                                                                                                              >  1:388344/1‑250 (MQ=255)
                                                                  aGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCt                                                                                                                                                                              >  1:345111/1‑250 (MQ=255)
                                                                   gCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGCAGGAACCTTTGGCCTGGATAAAATCTACGAACTGTGCGAAATCGCCCGTGCAGAGCTGACCAACGGCGGATAc                                                                                                                                                                            >  2:48938/1‑251 (MQ=255)
                                                                    cTTCCCGCAAGAGCTGCTGGCTAAACTGGCCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGACACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTACGAAATCGCCCGTGAAGCGCTGACCAACGGCGGCtcc                                                                                                                                                                            >  2:651832/1‑248 (MQ=255)
                                                                    cTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTAc                                                                                                                                                                            >  1:106666/1‑250 (MQ=255)
                                                                       cccGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACaacc                                                                                                                                                                        >  1:172989/1‑249 (MQ=255)
                                                                            aaggtttgctgGATAACCGGGCCGACCGCGCTAATGTGCCGG‑TGTCCCCCGTAACCTGCCACTCTTCCGTTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACGCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGt                                                                                                                                                                   <  2:21994/245‑1 (MQ=255)
                                                                                 ctaatagATAAACTGGTCGAACGAGCTAATCTGCCGGGTTACCTCGGTAAA‑TGCCACTCTTCCGGTACGGGCATTCTGGATCAGCTGGGCGCAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGGTCAAGATTGCCTGCCACGAAGGAACCTTGGGCCTGAATAAACTCTACGAACTGTGCGGAGTCGCCCGTGAAGAGCGGACCAACTGCTGCTaacccatcggtagggt                                                                                                                                                              >  2:160310/7‑236 (MQ=255)
                                                                                    gtgGATAAACGGGCCGAACGCGTCAATCTGCCGGTTTGCCTCGGAAAC‑TGACACTCTTCCGGTACGGTCATTCTGGATCAGCGGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAACATCGGTCGTGT‑TAt                                                                                                                                                          <  2:508139/250‑1 (MQ=255)
                                                                                    ctgGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAAcggc                                                                                                                                                                                  >  1:146440/1‑228 (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)
                                                                                                                             cGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAGCTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCACGAAGGAACCTTTGGCCTGGATa                                                                                                                                                                                                                                 <  1:127671/140‑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)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
GGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAAC‑TGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGT‑TATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTG‑TTGAGCCGCCAGCACCGAC‑CGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGA  >  NC_000913/4619903‑4620387

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