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,4010CT90.9% 149.6 / ‑2.3 68G204S (GGC→AGC) yjhHputative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase
Reads supporting (aligned to +/- strand):  ref base C (6/0);  new base T (19/41);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 1.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTG  >  NC_000913/4524160‑4524644
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
tcttcAGATGCTTCAAGAATCGGTTGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTAGTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCTGTTGTAGTAGTTTTTTATTCAGCGTCTCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGTCGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAgcagca                                                                                                                                                                                                                                               <  1:547635/250‑1 (MQ=255)
   tcAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCg                                                                                                                                                                                                                                            <  1:218246/250‑1 (MQ=255)
    cAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGCGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCgg                                                                                                                                                                                                                                           >  1:482324/1‑249 (MQ=255)
    cAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGt                                                                                                                                                                                                                                           <  1:341205/250‑1 (MQ=255)
       aTGCTTCAAGAATCGGTGGTAAGCAATAAGTCTCGACAGGCAGCCCGACACACTGCATGCTGGATTTGATCAGTGAGACAAAAGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGGAGTTTTTTAGTCAGCGAAGCAGCGGTAGCCAGATCGCCTTCCCGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTAGCGCACCGTCGCCGCTCAGCAGCATCGGAtta                                                                                                                                                                                                                                       <  2:118153/251‑2 (MQ=255)
        tGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCCCTCAACAGCATCGTATTc                                                                                                                                                                                                                                       >  1:209405/1‑250 (MQ=255)
              aaGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCCTCGCCGCTCAGACGCATCGTATTccgccaat                                                                                                                                                                                                                                >  1:205381/1‑244 (MQ=255)
                 aaTCTGTGGTAAGCAATACGTCTCTCCAGGCAGCCCGACACACTGCATGTTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGat                                                                                                                                                                                                                             <  2:129373/251‑1 (MQ=255)
                 aaTCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGat                                                                                                                                                                                                                             <  1:134682/251‑1 (MQ=255)
                      gtggtAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCCTCGTATTCAGCAAAAGATCATCg                                                                                                                                                                                                                        >  1:55543/1‑251 (MQ=255)
                      gtggtAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAACCAAAGATCATCg                                                                                                                                                                                                                        >  1:506716/1‑251 (MQ=255)
                         gtAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGGCGCCGCTCAGCCGCCTCGTATTCAGCCAATGATCATCGTa                                                                                                                                                                                                                      >  1:82352/1‑250 (MQ=255)
                                atatGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTCAGCGCTGGCGGTTATCGCCCCGTCGCCGCTCAGCAGCATCGTATTCAACAAATTATCATCGTAACCGccg                                                                                                                                                                                                               >  1:172984/1‑248 (MQ=255)
                                        ctACAGGCTGCACGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTGTCGAGGGCGGAAATAGCGGGCAGGTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGAGCGTCTTCACGCCAGGCGCGGTATATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCTTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATAcc                                                                                                                                                                                                         <  2:154238/248‑1 (MQ=255)
                                        ctACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATAcc                                                                                                                                                                                                         <  1:45377/248‑1 (MQ=255)
                                          aCAGGCTGCCCGACACACGGCAGGCTGTATTAGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTAGACAGCAAATGATCATCGTAACCGCAGAATTCCGaaa                                                                                                                                                                                                     <  2:157753/250‑1 (MQ=255)
                                          aCAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGAAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGttaacgacgaataccgcaa                                                                                                                                                                                                     >  1:361970/1‑232 (MQ=255)
                                                            tacaacctttATTTGAGCAGTCCGACAAGCGGTGGTTTGAGTGCGCAAATAGCGGGCGGTTGTAGTATATTGTTAGTCAGCGGCGCAGCGGTATACAGATCGCCTTCGCTCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAGGTTAGCGCTGGCGGTTACCGCATCGTCGCCGCTCAGCAGCATCGTATGCAGCAAATGATCATCGTAACCGCAGAATATCGAAAACGACGGGCGTACCgta                                                                                                                                                                                    <  2:311144/240‑3 (MQ=255)
                                                              cATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACTACGGGCGTACCGACTTaa                                                                                                                                                                                <  1:169843/251‑1 (MQ=255)
                                                                     tATTTGATCAGTGAGCCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGCCGCCGCTCCGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAAATGTGtt                                                                                                                                                                        >  1:443515/1‑251 (MQ=255)
                                                                      aTTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGtt                                                                                                                                                                        <  1:197239/250‑1 (MQ=255)
                                                                                 gagaCAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATAGCACCCTCGCCGCCCTGCCGCATCCTATTCCGCAAATAATCATCATAACCCCAGAATACCGAAAACGACGGGCGGAACGACTTTACTGTGTTGGTCATCgaacc                                                                                                                                                            >  1:193595/1‑247 (MQ=255)
                                                                                      aaaCGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTGATAGCACCCACGCCGCACAGCCACATCCTATTATGCAAAAGATCCTCGGAACCCCACAAGACCGAAAAAGAAGGGGGTACAGACTTAAcatagttgatactcggcaccaagg                                                                                                                                                        >  2:299846/1‑228 (MQ=255)
                                                                                             gtTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAAc                                                                                                                                               <  1:236098/251‑1 (MQ=255)
                                                                                             gtTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTGAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTATGCAAGTGACCAAc                                                                                                                                               <  1:35092/251‑1 (MQ=255)
                                                                                                  gAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTTTCGCACCGGCCCCCCCCCGCCGCATCCCATTCAGCAAATGCTCCTCGTCACCGCAGAATACCGCAAACGACGGGCGGCCCCACTTAACTGTGTTGatcatccttcccacgtggccaacgcctc                                                                                                                                           >  2:52938/1‑229 (MQ=255)
                                                                                                      gCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGAt                                                                                                                                      <  1:384153/251‑1 (MQ=255)
                                                                                                       cGTAAATAGCGGGCAGTTGTAGTAGGTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATg                                                                                                                                     <  2:567663/251‑1 (MQ=255)
                                                                                                           aaTAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGgtgt                                                                                                                                 <  1:460280/251‑1 (MQ=255)
                                                                                                            aTAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGatc                                                                                                                                                                    <  1:352751/216‑1 (MQ=255)
                                                                                                              aGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCtt                                                                                                                              <  1:72700/251‑1 (MQ=255)
                                                                                                                 gggCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCCCTCCACAGCATCGTATTCCGCAAATTATCATCCTAACCGCCGAAAACCGGAAACGACCGGCGTACCGACTTAACTGTTTTGGTCCTCGTACGCAAGTTACCACCGCTGTCTATGTTGTCTTTgt                                                                                                                           >  1:302128/1‑250 (MQ=255)
                                                                                                                       ttGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGGAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGGTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGGTCATCGTAACCGCAGAATACCGAAAACGACGGGCGGACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCggaggggtattttgtgcca                                                                                                                      >  1:275736/1‑232 (MQ=255)
                                                                                                                         gtagtaGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTGCGCGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAAc                                                                                                                    <  2:412874/250‑1 (MQ=255)
                                                                                                                           agtagtTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCCGAATACCGAAAACGACGGGCGGACCGACTTAACTGTGTTGATCATCGTAAGCAAGTGACCAACGCTGTCGATGGGGTCTTTGATACCAACGat                                                                                                                 >  1:248518/1‑251 (MQ=255)
                                                                                                                            gtagtTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGGTAATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGGCTTAACTTAGGTGATTATTGTACGCAAGTGCACAACGCCGTCGTTGGTGTTTTTTATGCCAACGaa                                                                                                                 >  1:445909/1‑249 (MQ=255)
                                                                                                                              agtTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATAtt                                                                                                              <  1:11211/251‑1 (MQ=255)
                                                                                                                                  ttttATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTACGGAGCAGAGTTAGCGCAGGCGGTTATCGGACCGTCGCCGCGCAGCAGCCTCGTATTCAGCAAAAGATCATCGTAACCGCAGAAAACCGGAAACAGCGGGGGTAACGTATTAACATGGGGGATCATCGGACGAACGCTACCAACGCTGTCGGTGGGGTCTTTGAAGCCAAcaagaatctcg                                                                                                          >  1:517076/1‑241 (MQ=255)
                                                                                                                                  tgttATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATtctc                                                                                                           <  2:504626/248‑1 (MQ=255)
                                                                                                                                     tATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGGCACCTCGTATTTAGCAAATGATCATCGTAACCGCCGAATACCGAAAACGACGGGCGGACCGACTTAAATGGTTTGATCATCGTACGCAAGTGACCAAAGCTGTCGCTGGGTTCTTTGTTGTCAACGGTATTCTCGttt                                                                                                       >  1:380416/1‑251 (MQ=255)
                                                                                                                                      aTTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTCTCGCCCCGTCGCCGCCCAGCAGCATCGTATTCCGCCAATGATCATCGAAACCCCCGAATACCGGAAAAGACGGGCGTACCCACTTAAATGTGTTGATCATCGTACGCAAGTGGCCAACGCTGGCGGTGGGGGCTTTTATGCCAACGATATTCTCGtttt                                                                                                      >  1:477648/1‑251 (MQ=255)
                                                                                                                                      aTTCAGCCGCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATAGTCTCGttt                                                                                                       <  2:254307/250‑1 (MQ=255)
                                                                                                                                         cAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCCCCGTCGCCGCTCAGCACCATCGTATTCCGCAAATGATCATCGGAACCGCAGAATACCGAAAACGACGGGCGTAACGACTTAACTGTGTTGATCATCGTACGCCAGTGACCACCGCTGTCGCTGGGGGTTTTTTTGCCCACGGTATTCTCATTTTGcc                                                                                                   >  2:155726/1‑250 (MQ=255)
                                                                                                                                            cGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCTCCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCTCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGTCCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCagag                                                                                                <  1:469691/251‑1 (MQ=255)
                                                                                                                                              gcgcctcgggcGCCAGATCGCCTTCACGACAGGCGCGGTAGATTCCGACGGAAAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCAACGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGcc                                                                                              <  1:380301/241‑1 (MQ=255)
                                                                                                                                                ggaGCGGGCGCCAGAGCGCCGTCACGCCAGGCGCGGTAGAGGCCGACGGAGAGTTCCGGAGCAAAGGTAGCGCGGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAg                                                                                            <  2:302128/249‑1 (MQ=255)
                                                                                                                                                     ggTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGAcgcg                                                                                       <  2:506716/251‑1 (MQ=255)
                                                                                                                                                     ggTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGAcgcg                                                                                       <  1:155726/251‑1 (MQ=255)
                                                                                                                                                      gTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCCCCGGCGCCGCTCAACAGCATCGTATTCCGCAAATGATCATCGTCAACGCAGAATACCGAAAACGACGGGCGGCCCGACTTAACTGTGTTGATCATCGTACGACAGTGACCAACGCTGTAGATGGGTTCTTTGATTACAACGCTATTTTCGTTTTGCCGAGACCGAcgcgg                                                                                      >  2:508455/1‑250 (MQ=255)
                                                                                                                                                        cGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCTCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCa                                                                                    <  2:522392/251‑1 (MQ=255)
                                                                                                                                                         gCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCCGCAGCATCGGATTTAGCAAATAATCATCGTAACCGCAGAATACCGAAAACGACCGGCGTACCGACTTTACTGTGTTGATCATCGGAACCAAGATGACAACGCTGTCGATGGTGGTTTTGATACCCACGATATTTTCGATTTTCAGAGCCAGAAGCGTCa                                                                                    >  2:322705/1‑250 (MQ=255)
                                                                                                                                                            aGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGt                                                                                <  1:52938/251‑1 (MQ=255)
                                                                                                                                                            aGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGt                                                                                <  1:370406/251‑1 (MQ=255)
                                                                                                                                                              aTCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGGCGCCGCTCAGCAGCCTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACAGGCGTACCGAATTAACTGTGTTGATCAACGTAAGCAAGTGAACAACGCTGTCGATGGGGTCTTTGATGCCAACGATATTCTCGGTTTGCAGAACCAGACGCGTCACgggtt                                                                              >  1:314960/1‑248 (MQ=255)
                                                                                                                                                                                gcgGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGGCGCCGCTCAGGAGCATCGTATTCAGCAAATGATAATCGGAAACGCAGAATAACGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGGTGGTGGCTTTGATTCCAACGATATTATCGTTTTGCAGAGCCAGACGCGGCACGGGTTCCGGGGTTAAGTCCTGAc                                                            >  1:165875/1‑251 (MQ=255)
                                                                                                                                                                                          ccGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAATTCCTGACCCGTCAGAtt                                                  <  2:321783/251‑2 (MQ=255)
                                                                                                                                                                                            gACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCCGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTAAGCCAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGAttcg                                                >  1:513785/1‑248 (MQ=255)
                                                                                                                                                                                            gACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGAttcg                                                <  1:299846/251‑4 (MQ=255)
                                                                                                                                                                                                   agTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATAGTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGtt                                         <  1:254907/251‑1 (MQ=255)
                                                                                                                                                                                                       ccGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCATCGCTGTCGATGGTGTCTTTGATGCCAACg                                                                                                                   <  1:388008/173‑1 (MQ=255)
                                                                                                                                                                                                        cGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACa                                    <  2:209405/251‑1 (MQ=255)
                                                                                                                                                                                                        cGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTTCGCAAATGATCATCGGAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGGCCAACGCTGTCGGTGGGGTCTTTGATGCCAACGAAATTCTCGTTTTGCAGAGCCCGACGCCTCCCGGTTTTCGGGGTTAAGTCCTGGCCCGGCAGaattggaaaggtgtta                                     >  1:309323/1‑235 (MQ=255)
                                                                                                                                                                                                          gAGCAAAGTTAGCGCTGGCGGTTATCGCAACGGCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCCTCCTAACCGCAGAATACCCAAAACGACCGGCGTACCGCCTTAACTGTGTTTATCATCGTACGCAAATTACCACCGCTGTCGATGGTGTCTTTGATTCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCCCGGTTTCCGGGGTTTAGTCCTGAcccctcctattctgcaagttttcccgg                                  >  1:124213/1‑227 (MQ=255)
                                                                                                                                                                                                               aaggtAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATAATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTTCAGGATCAc                             <  2:127818/247‑1 (MQ=255)
                                                                                                                                                                                                                            cGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACg                <  1:508455/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cACCGTCTCCGCTCAGCAGCATCGTATTCATCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTTATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGGTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGGCCTGTCCCGGCAGATTCGGAAAGTTGTACAGGATCACCGGTGGGGTGACGCTACGGGCg       <  2:368330/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cccctCGCCGCTCAGCAGCATCGTATTCAGAAAATGATCATCGTAACCGCAGAATACCTAAAACGTCGGTCGTACCTACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGGCGATGGTGTCTTTGATGCCAACGTTATTCTCGTGTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCGTACCCGTCAGATTCGGAAAGGTGTACAGGAACACCGGTAGGGTGACTCTACGGGCGa      <  2:443515/247‑1 (MQ=255)
                                                                                                                                                                                                                                           cgccgcTCAGCAGCATTGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCCGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATctg  <  1:532767/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
TCTTCAGATGCTTCAAGAATCGGTGGTAAGCAATATGTCTCTACAGGCAGCCCGACACACTGCATGCTGTATTTGATCAGTGAGACAAACGGTGTTTCGAGGGCGTAAATAGCGGGCAGTTGTAGTAGTTTTTTATTCAGCGTCGCAGCGGTCGCCAGATCGCCTTCACGCCAGGCGCGGTAGATGCCGACGGAGAGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTG  >  NC_000913/4524160‑4524644

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