Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,524,596 C→T D139N (GAT→AAT)  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,5960CT95.7% 204.0 / ‑2.5 69D139N (GAT→AAT) yjhHputative 2‑dehydro‑3‑deoxy‑D‑pentonate aldolase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (26/40);  total (27/42)
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.99e-01

AGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGC  >  NC_000913/4524355‑4524838
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
agTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATAGTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGtt                                                                                                                                                                                                                                            <  1:254907/251‑1 (MQ=255)
     cGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCTCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACa                                                                                                                                                                                                                                       <  2:209405/251‑1 (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)
                                           cgctcagcagCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATc                                                                                                                                                                                                                                  <  1:96152/214‑1 (MQ=255)
                                           cgctcagcagCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATc                                                                                                                                                                                                                                  >  2:96152/5‑218 (MQ=255)
                                            gctccgcagcaTCGTATTCATCAAAATATCATCGGAACCGCAGAATACCGAAAACGACTGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAATTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGTATCACCGGTAGGGTGACGCTTCTGGCGATCTGCTGGt                                                                                                                                                                                               <  2:517076/246‑1 (MQ=255)
                                            gctcagcagCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGt                                                                                                                                                                                               <  1:394271/248‑1 (MQ=255)
                                               cagcagAATCGTATTCAGCAAAGGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGGGACGCTACGGGCGATCTGCTGGTAAt                                                                                                                                                                                            <  2:82352/251‑1 (MQ=255)
                                                agcagcATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAAt                                                                                                                                                                                            >  1:61579/1‑250 (MQ=255)
                                                    gcaTCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAgtcg                                                                                                                                                                                       <  1:129306/251‑1 (MQ=255)
                                                                 cAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCgtg                                                                                                                                                                          >  1:364404/1‑251 (MQ=255)
                                                                  aaaTGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACACGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCgtgg                                                                                                                                                                         >  1:304889/1‑251 (MQ=255)
                                                                          agcGTAACCGCAGAATACCGAAAACGACGGGCGTACCTACTTAACTGTGGTGATCATCGGACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACtt                                                                                                                                                                 <  2:481487/249‑1 (MQ=255)
                                                                           tcGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCg                                                                                                                                                                                                          >  1:301174/1‑210 (MQ=255)
                                                                           tcGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCg                                                                                                                                                                                                          <  2:301174/210‑1 (MQ=255)
                                                                             gtaatcgtaGAATATCGAAAACGACGGGCGTACAGACTGAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGGCGATGGTGGCTTTGATGCCAACGATATTCTCGTGTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTcc                                                                                                                                                              <  2:11705/243‑1 (MQ=255)
                                                                                  cGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCagtgg                                                                                                                                                         >  2:514391/1‑249 (MQ=255)
                                                                                   gCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCagtagt                                                                                                                                                        >  1:524501/1‑251 (MQ=255)
                                                                                      gAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAgg                                                                                                                                                     >  1:74088/1‑251 (MQ=255)
                                                                                        aTACCGAAAACGACGGGCGTACCGCCTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCCAGGTTTCGTGGTGCGACTTTCCAGTAGTAgggg                                                                                                                                                   >  2:155008/1‑251 (MQ=255)
                                                                                          gccGAAAACGACGGGGGTACCGACTTAACGGTGTGGATCATCGGACGCAAGTGACCAACGCTGTCGATGGGGTCTGTGATGCCTACGATATTCCCGTTTTTCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGt                                                                                                                                                  <  2:482324/249‑1 (MQ=255)
                                                                                           ccGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTg                                                                                                                                                <  1:322705/251‑1 (MQ=255)
                                                                                                  cgacgGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGtagtag                                                                                                                                                      >  1:432791/1‑238 (MQ=255)
                                                                                                  cgacgGGCGTACCGACGTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGtagtag                                                                                                                                                      <  2:432791/238‑1 (MQ=255)
                                                                                                   gacgGGCGTACAGACTTAACTGTGTTGATCATCGTACGAAAGTTACCAACGCTGTCGATGGAGTCTTTGATGCCCACGATATACTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGACAGCTTCGGAAAGTTGTACAGGATCACCGGTAGGTTGACGCTTAGGGCGTTCTGCTGGTAATAGTCGACAAGATATCGTGGTGCGACTTTCCATTAGTAGGGGTTGATGGcgac                                                                                                                                        >  2:181487/1‑251 (MQ=255)
                                                                                                    acgGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGcgacg                                                                                                                                       <  1:481721/251‑1 (MQ=255)
                                                                                                            aCCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCagta                                                                                                                                                          <  1:54865/224‑1 (MQ=255)
                                                                                                            aCCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCagta                                                                                                                                                          >  2:54865/1‑224 (MQ=255)
                                                                                                                 cTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTACGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCAGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGGTGATGGCGACGATCCCATCAgcgc                                                                                                                          >  1:63025/1‑251 (MQ=255)
                                                                                                                     aCTGTGTTTATCATCTCACGCACGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCTTTGTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGTGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTa                                                                                                                      <  2:165875/251‑1 (MQ=255)
                                                                                                                              atcatcGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGGTGGCGATGATACCATCAGCGACGTAGGCTTgcg                                                                                                              >  2:393045/1‑250 (MQ=255)
                                                                                                                                     tACGCAAGTGACCAACGCTGTCGATGATGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAACCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAgc                                                                                                                            >  1:140482/1‑229 (MQ=255)
                                                                                                                                     tACGCAAGTGACCAACGCTGGCGAGGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGTCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCTGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAgc                                                                                                                            <  2:140482/229‑1 (MQ=255)
                                                                                                                                           aGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGt                                                                                                >  1:34533/1‑251 (MQ=255)
                                                                                                                                                ccAACGCGGGCGAGGGTGTCTTGGATGCAAACGATATTCGCGTTGTGCAGAGCCAGGTGCTGAACGGTTTCCGGTGTTAAGACGAGAGCCGTCAGGATCGGAAAGGTTTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGGCAAGATGTCGTGGGGCGACTTTCCAGTAGTAGGCGTGGATGGCGACGATGCCATCAGCGCCGTAGGGTTGCGCATGCGGCGCCCGTTCGAc                                                                                           <  2:193595/251‑1 (MQ=255)
                                                                                                                                                  cccctattccgatctatcggtGATGCCAACGATATTCTCGTTGTGCAGAGCCAGACTCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGGCAGATTCGGAAAGTTGTACAGGAGCACCGGTAGGGTGACGCTACGGGCGATTTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCCGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGaacg                                                                                         <  2:323135/231‑4 (MQ=255)
                                                                                                                                                  aaCGCTGTCGATGGTGTCTTTGATGCCAACTATAGTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTTATGGCGACGATACCATCAGCGTCGTAGGCTTGCGCATGCTGCGCCAGTTGGACCg                                                                                         <  1:164829/251‑1 (MQ=255)
                                                                                                                                                  aaCGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGTTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTGGACCg                                                                                         >  2:418855/1‑251 (MQ=255)
                                                                                                                                                    cgatGTCGATGGTGTCTCTGATGCCAACGAGATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCt                                                                                       <  2:172984/248‑1 (MQ=255)
                                                                                                                                                       tGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCg                                                                                    >  1:274576/1‑251 (MQ=255)
                                                                                                                                                           gATGGTGTCTTTGATGCCAACGATATTCTCGTGTTTCAGAGCCAGCCGCGTCACGGTTTCCGGGGTTAAGGCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGGCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAg                                                                                <  2:275736/251‑1 (MQ=255)
                                                                                                                                                               gtgtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGAGACCATCAGCGCCGTAGGCTTGCGCATGCTGCGTCAGTTTGACCGCTTCGTCAGTGGa                                                                            <  1:463962/251‑1 (MQ=255)
                                                                                                                                                                 gtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCg                                                                                         >  1:323135/1‑236 (MQ=255)
                                                                                                                                                                       gATGCCAACGATATTCTCGTGTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAAc                                                                    <  2:55543/251‑1 (MQ=255)
                                                                                                                                                                          gCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGa                                                                 <  1:514391/251‑1 (MQ=255)
                                                                                                                                                                           ccAACGATATACTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCAGCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGAc                                                                <  2:205381/251‑1 (MQ=255)
                                                                                                                                                                            cAACTATATTCTCGTCTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGAGTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACGTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACg                                                               <  2:380416/251‑1 (MQ=255)
                                                                                                                                                                                gatattcaagcattGCAGAGCCAGACGCGTCACGGCTTCCGGGGTAAAGTCCTGACCCGACAGATTCGGAACGTTGTACAGGATCACCGGTAGGGAGACGCTACGGGCGATCTGCTGGTAATAGTCGTCCCGATTTCGTGGAGAGACATTACAGTAGTAGGGGATGATGGCGACGATAACATCACCCACGTCGGCTTGCGCATGCTCCGCCAGGTTGACCGCTTCGTAAGGGGAAGGGGgacggacggcac                                                           >  2:75672/13‑247 (MQ=255)
                                                                                                                                                                                 atatTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTgg                                                                             <  2:169740/232‑1 (MQ=255)
                                                                                                                                                                                 atatTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTgg                                                                             >  1:169740/1‑232 (MQ=255)
                                                                                                                                                                                                   gccaacagcgTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCTGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGccc                                        <  1:3170/244‑1 (MQ=255)
                                                                                                                                                                                                              cACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATg                              <  1:482383/250‑1 (MQ=255)
                                                                                                                                                                                                              cACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGGGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATgg                             <  1:135072/251‑1 (MQ=255)
                                                                                                                                                                                                                  gattCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAAACGACGCCAATCAATACCGGCCCTTGCCCGCCGCCAATGGTTAc                         >  2:183491/3‑251 (MQ=255)
                                                                                                                                                                                                                   tttCCGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTAc                         >  2:77608/1‑250 (MQ=255)
                                                                                                                                                                                                                    ttcaggggTTAAGGCCGGACCCGTCAGATTCGGAAAGTTGTACTGGTACACTTGTAGGGTGACGCTACGGGCGTTCTGCTGTGTATTGTCGTCAAGATTACGTGGTGCGGCTTTCCAGTAGTAGGTGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCCATACCGGCACTCGCCCGTGGACAATGGTTACa                        <  2:310015/246‑1 (MQ=255)
                                                                                                                                                                                                                      ccGGGGTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGct                     >  1:260993/1‑251 (MQ=255)
                                                                                                                                                                                                                        ggggTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGTTAGGGTGACGCGACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTGCCAGTAGTAGGGGTTGATGGCGACGATACCATCATCGCCGTAGGCTTGCGCATGCTGCGCCCGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATCCCGGCACTCGCCCGTCGACAATGGTTACAGcttc                   <  1:75672/251‑1 (MQ=255)
                                                                                                                                                                                                                        ggggTTAAGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGcttc                   >  1:456157/1‑251 (MQ=255)
                                                                                                                                                                                                                              aaGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGc             >  1:183151/1‑251 (MQ=255)
                                                                                                                                                                                                                               aGTCCTGACCCGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCg            <  1:77608/251‑1 (MQ=255)
                                                                                                                                                                                                                               aGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGtt                                                                                               >  1:292114/1‑168 (MQ=255)
                                                                                                                                                                                                                               aGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCGACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGtt                                                                                               <  2:292114/168‑1 (MQ=255)
                                                                                                                                                                                                                                        ccGTCAGATTCGTAAAGTTGTACAGGAGCACCGGGAGGGTGACGCTACGTGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGGGCGACTTTCCAGTAGGCGGGGTTGATTTCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCTTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCGTCTTCGGCGAGTGCCa     <  2:61579/249‑1 (MQ=255)
                                                                                                                                                                                                                                         cGTCAGATTCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATgc  <  1:393045/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
AGTTCCGGAGCAAAGTTAGCGCTGGCGGTTATCGCACCGTCGCCGCCCAGCAGCATCGTATTCAGCAAATGATCATCGTAACCGCAGAATACCGAAAACGACGGGCGTACCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCTGGTAATAGTCGTCAAGATTTCGTGGTGCGACTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGC  >  NC_000913/4524355‑4524838

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