Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,524,401 C→T G204S (GGC→AGC)  yjhH ← putative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,524,4010CT97.9% 134.5 / NA 48G204S (GGC→AGC) yjhHputative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/30);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTG  >  NC_000913/4524155‑4524644
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
ttGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTcag                                                                                                                                                                                                                                                  <  2:147948/250‑1 (MQ=255)
           gATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTAtt                                                                                                                                                                                                                                      <  1:127674/251‑1 (MQ=255)
                   aaGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGGTATCGCAGCGTCGCCGCTCAGCAGCATCGTATTCATCAAAt                                                                                                                                                                                                                              <  1:372977/251‑1 (MQ=255)
                      aaTCGGTGGTAAGCAATATGTCTCTACAGGCATCCCGACGCACTGCATGCTGTATTTGATCAGTGAGACAGACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTAGTCAGCGTCGCAGCGGTCGCCTGATCGCGTTGACGCCAGGCGCGTGAGAGGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGat                                                                                                                                                                                                                           <  1:453111/251‑1 (MQ=255)
                            tggtAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCCCCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGt                                                                                                                                                                                                                     >  1:465538/1‑251 (MQ=255)
                             ggtAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATGCAGCAGATGATCATCGt                                                                                                                                                                                                                     <  2:465538/250‑1 (MQ=255)
                                       atGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAAt                                                                                                                                                                                                          <  1:49684/251‑1 (MQ=255)
                                       atGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAAt                                                                                                                                                                                                          <  1:102521/251‑1 (MQ=255)
                                                     aGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCCCCGTCGCCGCTCCGCAGCCTCGTATTCCGCCAAAGATCCTCGTAACCGCAGAATACCGAAAACGACgg                                                                                                                                                                                            >  1:199794/1‑251 (MQ=255)
                                                            cacacTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTACATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCCCCGTCGCCCCTCCCCACCATCGTATTCAGGAAAAGATCATCGTTACCGCCGAATAcagcaaacaaagggcggccc                                                                                                                                                                                     >  2:428895/1‑232 (MQ=255)
                                                                        tGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACtgtg                                                                                                                                                                         >  1:142173/1‑251 (MQ=255)
                                                                                   aGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTAATCAGCAAATGATCATCGTAACCGCAGAATACCGAAAAATACGGGCGTACCGACTTACATGTGTTGATCATCGt                                                                                                                                                              >  1:200103/1‑251 (MQ=255)
                                                                                   aGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCGTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGt                                                                                                                                                              <  1:44288/251‑1 (MQ=255)
                                                                                       agaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGc                                                                                                                                                          <  1:1670/251‑1 (MQ=255)
                                                                                               ggTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTTGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTAGTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTGCGCAAGTGAcc                                                                                                                                                  <  1:250555/251‑1 (MQ=255)
                                                                                               ggTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATAGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTCTTCAGCAAATGATCATCGTAAACGCCGAATACCGAAAACAGCGGGCGGAACGACTTCACTGTGATGATCCTCGTACGCAAGGGAcc                                                                                                                                                  >  2:344107/1‑251 (MQ=255)
                                                                                                        aGGGCGTAAATAGCGGGCATTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCTCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCCAGTGACCAACGCTGTc                                                                                                                                         <  1:1292/251‑1 (MQ=255)
                                                                                                        aGGGCGCAAATAGCGGGCCGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCCGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGGGACCAACGCTGTc                                                                                                                                         >  1:349446/1‑251 (MQ=255)
                                                                                                         gtgCGTAAAGAGCGTGCAGTTGTATGAGTTTTTTATGCAGCGTCGCAGCGGTCGTCAGATCGCCTGCATTCCTGGCGCGGTAGACGCCGTCGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTGATCGCTCCTTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATGCCGAAATCGACTGGCGGTCCGACTTAACGGTGTTGATCTCCGTACGGTAGGGACCACCGCTGTCg                                                                                                                                        <  2:441121/249‑1 (MQ=255)
                                                                                                                   ggCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGATTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATTCCGAAAACGACGGGCTTACCGACTTATGTGTGTTGACCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCtt                                                                                                                              <  1:51254/250‑1 (MQ=255)
                                                                                                                      gggCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGa                                                                                                                           <  1:428895/251‑1 (MQ=255)
                                                                                                                       ggCAGTTGGAGTAGGTTTTTAGTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCTACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGAt                                                                                                                          <  2:96990/251‑1 (MQ=255)
                                                                                                                              gtagtaGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGGACGCAAGGGACGAACGCTGTCTATGGTGTCTTTGATGCCAAc                                                                                                                    <  2:319501/250‑1 (MQ=255)
                                                                                                                                      tttttATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGGACGCAAGTGACCAACGCTGTCGATGGGGTCTTTGATTCCAACGATATtctc                                                                                                           >  1:44281/1‑251 (MQ=255)
                                                                                                                                       tgttATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATCCCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCGTTGATGCCAACGATATtctc                                                                                                           <  2:253728/248‑1 (MQ=255)
                                                                                                                                         ttATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGAGATTCTCGt                                                                                                         <  1:344107/250‑1 (MQ=255)
                                                                                                                                          tATTCCGCGTCGCCGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCCCAACGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGGAAACGACGGGCGTACCGACTTAACTGTGGTTCTCATCGTACGCAAGGGACCAACGGTTGCGGTGGGGTCTTTGATGCCAACGATTTTCTCGttt                                                                                                       >  1:372053/1‑251 (MQ=255)
                                                                                                                                                gCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGAGGTTATCGCACCGTCGCCGCTCAGCAGCATCGCATTCAGCAAATGATCATCGTa                                                                                                                                                                                                                    <  2:428945/136‑1 (MQ=255)
                                                                                                                                                                     cGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTcc                                                                            <  1:290461/251‑1 (MQ=255)
                                                                                                                                                                      gCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGGGTTGATCCTCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGGTTCCg                                                                           >  2:390215/1‑251 (MQ=255)
                                                                                                                                                                                 aGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCGCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTcc                                                                <  1:381783/251‑1 (MQ=255)
                                                                                                                                                                                    cgcgGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCCTCGCCGCTCAGCAGCATCGTATTCCGCAAATGATCATCGTAACCGCAGAATACCGAAAAAGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGGTTCCGGGGTTAAGTCCTg                                                              >  2:345975/1‑250 (MQ=255)
                                                                                                                                                                                           gATGCCGACGGAGAGTTCCGGAGCAAAGGTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCggca                                                      <  2:66439/251‑4 (MQ=255)
                                                                                                                                                                                                cGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGGCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGAGTCACGGTTTCCGGGGTTAAGTCCGGACCCGTCAGAttc                                                 <  2:15629/251‑3 (MQ=255)
                                                                                                                                                                                                      agagTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCCGCAGCATCGTATTCCGCCAATGATCATCGCAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGGTCATCGTACGCCAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCCACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGaa                                             >  1:218249/1‑249 (MQ=255)
                                                                                                                                                                                                        agTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGtt                                         <  2:352354/251‑1 (MQ=255)
                                                                                                                                                                                                        agTTCCGGAGCAAAGTTAGCGCTGGCGGTTATAGCAACGTCGCCGCTCAGCAGCATCGTATTCAGCAAAAGATCATCGTAACCGCAGAATACCGAACACGGCGGGCGTAACGACTTAACTGTGTTGGTCATCGTACGCAAGTGACCAACGCTGTCGATGGGGTCTTTGATGCCAACGGTATTCTCGTTTTTCAGAGCCAGACGCGTCCCGGTTTCCGGGGTTAAGTCCTTGCCCGTCAAATTCGGAAAGtt                                         >  2:314551/1‑251 (MQ=255)
                                                                                                                                                                                                         gTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTg                                        <  1:390215/251‑1 (MQ=255)
                                                                                                                                                                                                                 gCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAgggt                                >  1:340984/1‑249 (MQ=255)
                                                                                                                                                                                                                         aGCGCTGGCGTTTAGCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACTACGGGCGTACCGACTTAACTGTGATGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTGTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACTGTTTCCGGGGTTAAGTCCTGACCCGTCAGTTTCGGAAAGTTGTACAGGATCACCGGt                         <  2:199794/250‑1 (MQ=255)
                                                                                                                                                                                                                         aGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTGTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTa                        <  2:150405/251‑1 (MQ=255)
                                                                                                                                                                                                                                gCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGAc                 <  1:314551/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGc               >  1:37444/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGCTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGc               >  2:252836/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cgctccGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTGACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTTTGTTGATCATCGTACGCAAGTGACCAACGCTGGCGATGGTGTCTTTGATGCCAACGATATTCTCGGTTTGCATAGCCAGAGGCGGCACGGGTGCCGGGGTTAATTCCTGACCCGTTTGATGCGGAAAGGTGTACTGTATCTCCGGTATGGttgcgttacgg          <  2:206141/246‑11 (MQ=255)
                                                                                                                                                                                                                                            ccGTCGCCGCTCAGCAGCATCGTATTCAGCACATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGAt     >  1:130455/1‑251 (MQ=255)
                                                                                                                                                                                                                                               tCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGACTACCGCAAACGACCGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCACCGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATctg  >  2:29015/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cgccgcTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATctg  <  1:182519/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
TTGCTTCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTG  >  NC_000913/4524155‑4524644

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