Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,726,004 G→A 29.8% intergenic (‑84/‑19) ydhL ← / → nemR DUF1289 family protein/transcriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,726,0040GA29.8% 76.7 / 39.0 55intergenic (‑84/‑19)ydhL/nemRDUF1289 family protein/transcriptional repressor for the nemRA‑gloA operon, quinone‑, glyoxal‑, and HOCl‑activated
Reads supporting (aligned to +/- strand):  ref base G (17/21);  new base A (7/9);  total (24/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 = 9.47e-01

GCTGGCTTATTTGCGCGTAATTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTA  >  NC_000913/1725718‑1726291
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                               
tctggtgtATTTGCGTGTAATTTGCCCATTAATCGCTGACGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTc                                                                                                                                                                                                                                                                                               <  2:376743/282‑1 (MQ=255)
gCTGGCTTATTTGCGCGTAATTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATc                                                                                                                                                                                                                                                                                               <  1:754264/289‑1 (MQ=255)
           tGCGCGTAATTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACTcc                                                                                                                                                                                                                                                                                    <  1:92456/289‑1 (MQ=255)
                gTAATTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCACTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACTATACAACg                                                                                                                                                                                                                                                                                >  2:56331/1‑288 (MQ=255)
                gTAATTTGCGCATTAATCGCGGCCGACAAAGGCGCAGCACCTCTGGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACg                                                                                                                                                                                                                                                                                <  2:417298/288‑1 (MQ=255)
                  aaTTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATGCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCa                                                                                                                                                                                                                                                                              <  1:56331/288‑1 (MQ=255)
                       gcgcATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGaaca                                                                                                                                                                                                                                                                        <  1:505111/289‑1 (MQ=255)
                        cgcATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACTCCAACGCa                                                                                                                                                                                                                                                                              <  1:5161/282‑1 (MQ=255)
                        cgcATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACTCCAACGCa                                                                                                                                                                                                                                                                              >  2:5161/1‑282 (MQ=255)
                           aTTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACTCCAACGCATGaacaaaca                                                                                                                                                                                                                                                                    <  1:195499/289‑1 (MQ=255)
                             taacagCTGCCGACAAAGGCGCAGCACCCCTTGGTTTTCGCCATCGCGAATTTTATTCCAGTTAAAACGCTCAGCACGACGACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAacaca                                                                                                                                                                                                                                                                  <  2:358959/284‑1 (MQ=255)
                             tAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAacaca                                                                                                                                                                                                                                                                  <  2:562895/289‑1 (MQ=255)
                             tAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCGCAGCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAacaca                                                                                                                                                                                                                                                                  <  1:27404/289‑1 (MQ=255)
                                                     caCCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACAt                                                                                                                                                                                                                                           >  1:284890/1‑288 (MQ=255)
                                                          ctTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATTTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACCCCAACGCATGAACAACCACACCGAACATGATACCCGCGAACATCTCCTg                                                                                                                                                                                                                                     >  2:699862/1‑289 (MQ=255)
                                                                tttCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACg                                                                                                                                                                                                                               >  1:62001/1‑289 (MQ=255)
                                                                 ttCGCCAACGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTCGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTAGTCTACTACACTCCAACGCATGAACAAACACACCGTACATGATACTCGCGAACATCTCCTGGCGACgg                                                                                                                                                                                                                              <  2:194737/289‑1 (MQ=255)
                                                                   cGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCTCGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGCAAGAACTCTAATTGATAAGCCACTTCGCCCTCCTCTGATAAGATTATTCCCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTACTTGCAGCAGACCGACTGGTCTACTACATGACAACGAATGAACAAACACCCAGAACATGATACCCGCGAACCTCTCCTGGCTACGGGc                                                                                                                                                                                                                            >  2:61184/1‑289 (MQ=255)
                                                                         cGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCa                                                                                                                                                                                                                       >  1:665349/1‑288 (MQ=255)
                                                                                   tttCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGc                                                                                                                                                                                                            <  1:696819/288‑1 (MQ=255)
                                                                                          ttAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCa                                                                                                                                                                                                                                                                              <  1:348139/216‑1 (MQ=255)
                                                                                          ttAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCa                                                                                                                                                                                                                                                                              >  2:348139/1‑216 (MQ=255)
                                                                                          ttAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACCACCGAACATCTCCTGGCGACGGGCGAACAACTTTGCCTGCAACGTgg                                                                                                                                                                                                     >  2:696819/1‑289 (MQ=255)
                                                                                                      acaCGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACAGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAGCTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGAt                                                                                                                                                                                         <  1:431304/288‑1 (MQ=255)
                                                                                                                 aaaaCAGCCGCTACAAACACCGCGTTCGTCAGACTGGACAATTCCCCGGCACGGGCTCTGTACGTGAAAGAACTCTAATTGCACCGCCACTTCGCCCTCCTCTGATAAGATTATTACCATTATTGAAGCTGTTAATTTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCAAACGCATGAACAAACACACCGAGCATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCg                                                                                                                                                                              <  1:418995/289‑1 (MQ=255)
                                                                                                                            gACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTaaa                                                                                                                                                                    <  1:40830/288‑1 (MQ=255)
                                                                                                                                                   tGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTc                                                                                                                                            <  1:804294/289‑1 (MQ=255)
                                                                                                                                                       aaaTTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCtt                                                                                                                                         <  2:284890/288‑1 (MQ=255)
                                                                                                                                                               cGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCact                                                                                                                                 >  1:128072/1‑288 (MQ=255)
                                                                                                                                                                   aCGGGCTCTGTACGGGAAATAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACttt                                                                                                                            <  1:620892/289‑1 (MQ=255)
                                                                                                                                                                        ctctGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGct                                                                                                                        <  1:328023/288‑1 (MQ=255)
                                                                                                                                                                                  gAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACTCCAACGCATGAATAAACCCACCGAACATGATACTCGCGAACATCCCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTTGTTTCACCGGGATGGGGCTAAGCGAATTACCATAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACGACTTTCGCTCTACAGAAGCg                                                                                                             >  1:780045/1‑289 (MQ=255)
                                                                                                                                                                                                       ccACTTCACCCTCCTCAGATAGGACTATTACCGTTATAGAAGCGATTAATGTCCAAAGTGTCCACGTTGCTGGCACGAGACCGACTGATCTACAACACTCCAACGCATGAACTAACACACCGAACATGATACTCGCGAACATCTCCGGGCAACAGGCGGTCAACGTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAATTGCCGAAAGGTTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAg                                                                                        <  2:780045/289‑1 (MQ=255)
                                                                                                                                                                                                                 ctcctcATATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTAcg                                                                              <  1:112899/289‑1 (MQ=255)
                                                                                                                                                                                                                  tccacAAATAAAGTTATTAACATTATTGAGACTATAAAAGTCCCAAGTAGCAACGTTGCGTTCACTAGACCTACTGGTCTACTACACTCCAACGCATGAACCAACACACCTAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTAcgc                                                                             <  2:437946/289‑1 (MQ=255)
                                                                                                                                                                                                                   cctcAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTAGCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGACCGTCATGAcgcc                                                                            >  2:692053/1‑289 (MQ=255)
                                                                                                                                                                                                                    ctcAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTAcgccg                                                                           <  1:664202/289‑1 (MQ=255)
                                                                                                                                                                                                                    ctcAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTAcgcc                                                                            >  1:437946/1‑288 (MQ=255)
                                                                                                                                                                                                                                        ttattGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCgactg                                                          >  1:449128/1‑286 (MQ=255)
                                                                                                                                                                                                                                        ttattGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCgactg                                                          >  1:349092/1‑286 (MQ=255)
                                                                                                                                                                                                                                        ttattGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCgactg                                                          >  1:688285/1‑286 (MQ=255)
                                                                                                                                                                                                                                          attGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTCGCACTAGACCGACTGATCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGActgactga                                                     >  1:114490/1‑289 (MQ=255)
                                                                                                                                                                                                                                                gCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGtt                                                  >  1:787889/1‑286 (MQ=255)
                                                                                                                                                                                                                                                           ccAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGct                                    <  2:62001/289‑2 (MQ=255)
                                                                                                                                                                                                                                                           ccAAAGTAGCAACTTTGCTTGCACTAGACCGACTGATCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGAGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCg                                    <  1:699862/289‑1 (MQ=255)
                                                                                                                                                                                                                                                              aaGTAGCAACTTTGCTTGCACTAGACCGCCTGGTCTACTACACTCCAACGCATGAACAAAAACACCGAACATGATACTCGCGAACATCTCCTGGCGACAGACGAGCAAATTTGCCTGCAACGTGGATTCACCGGGATGAGGCTAAGCGAATTACTAAAAAACGCAGAAGTGCAGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGACATACTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGACGAAg                                 >  1:542648/1‑289 (MQ=255)
                                                                                                                                                                                                                                                              aaGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCAc                                                                   <  2:152854/255‑1 (MQ=255)
                                                                                                                                                                                                                                                              aaGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCAc                                                                   >  1:152854/1‑255 (MQ=255)
                                                                                                                                                                                                                                                                    cAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCAt              >  1:624251/1‑286 (MQ=255)
                                                                                                                                                                                                                                                                           gcttgcACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACcgcg                     >  1:254161/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                             ttgcACTAGACCGACTGATCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGAc                   <  2:128072/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    aGACCGACTGCTCTACTACACTCCAACGCATAAACAAACACACCGAACATGAGACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGAGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACt           <  1:136168/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          aCTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCATCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCTACTGACTGAGTTGCTGCAATCCGGCCAAGGTAACTACCGCAACCGCGTACTGGCtta     >  2:315302/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                           cTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGTAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTtgctg                                              >  1:59041/1‑247 (MQ=255)
                                                                                                                                                                                                                                                                                           cTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGTAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTtgctg                                              <  2:59041/247‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             ggTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCGTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCttatta  <  2:349092/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                               
GCTGGCTTATTTGCGCGTAATTTGCGCATTAATCGCTGCCGACAAAGGCGCAGCACCTCTTGTTTTTCGCCATCGCTCATTTTATTCCAGTTAAAACGCTCATCACGACTACGAAAACAGCCGCGACAAAACCCGCGTTCGTCAGACTGGCAAATTCCCCGGCACGGGCTCTGGACGGGAAAGAACTCTAATTGCTCCGCCACTTCGCCCTCCTCAGATAAGATTATTACCATTATTGAAGCTGTTAATGTCCAAAGTAGCAACTTTGCTTGCACTAGACCGACTGGTCTACTACACTCCAACGCATGAACAAACACACCGAACATGATACTCGCGAACATCTCCTGGCGACGGGCGAGCAACTTTGCCTGCAACGTGGATTCACCGGGATGGGGCTAAGCGAATTACTAAAAACCGCTGAAGTGCCGAAAGGGTCCTTCTATCACTACTTTCGCTCTAAAGAAGCGTTTGGCGTTGCCATGCTTGAGCGTCATTACGCCGCATATCACCAGCGACTGACTGAGTTGCTGCAATCCGGCGAAGGTAACTACCGCGACCGCATACTGGCTTATTA  >  NC_000913/1725718‑1726291

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: