Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,617,928 G→A P202P (CCG→CCA)  deoC → deoxyribose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,617,9280GA92.4% 182.5 / 1.7 67P202P (CCG→CCA) deoCdeoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base A (30/31);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCA‑CGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCG  >  NC_000913/4617685‑4618168
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
ggTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCTACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAATCCagcgggc                                                                                                                                                                                                                                             <  1:174079/250‑1 (MQ=255)
ggTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCagcgggc                                                                                                                                                                                                                                             <  1:578196/250‑1 (MQ=255)
  tGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGggcggc                                                                                                                                                                                                                                          >  1:326521/1‑251 (MQ=255)
          tGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTAc                                                                                                                                                                                                                                  >  1:290454/1‑251 (MQ=255)
                gagGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGa                                                                                                                                                                                                                            <  1:509213/251‑1 (MQ=255)
                                  aatttACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCGGAAATCTCCATCAAAGCTGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATTGGCGTAGAAAAAACCGTTGGTTGCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATAtct                                                                                                                                                                                                          <  1:75174/247‑1 (MQ=255)
                                               aaGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGAt                                                                                                                                                                                             <  1:294121/251‑1 (MQ=255)
                                                  tGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGAGCGGACTTCATCAAAACCTCTACAGGGAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGGAGATGCGCAGAAATATCTCGACATCGCAGATGa                                                                                                                                                                                           >  1:275555/1‑250 (MQ=255)
                                                        acGAAACCGGCGAACTGAAAGACGAGGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTGCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTc                                                                                                                                                                                    <  1:427331/250‑1 (MQ=255)
                                                              ccGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCCGCGGGCGGCGTGCGTACTGCGGGAGATGCGCAGAAATATCTCGCCATCGCAGATGAACGGTTCGGAGCt                                                                                                                                                                              >  1:553632/1‑251 (MQ=255)
                                                               cGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTg                                                                                                                                                                             >  1:597324/1‑251 (MQ=255)
                                                                        gAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGc                                                                                                                                                                      >  1:639995/1‑249 (MQ=255)
                                                                         aaaGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAACAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGTAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGc                                                                                                                                                                      >  2:44935/1‑248 (MQ=255)
                                                                         aaaGACGAAGAGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGCAAAGCGCGCGCATCATGATGGAAGTGATCCGTGCTAAGGGCGTAGAAAAAACCGATGGTTTCAAACCAGCGGGCGGAGGGCGTACTTAGGAAGACGAGCAGAAATATCTCGCCAACCCAGGTGAACTGTGCGGTGCCGACTgagaaa                                                                                                                                                                    >  2:364265/1‑245 (MQ=255)
                                                                          aagaagaagCGCTGATCAGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGCAGACGGACAGCGAGCGCATCATGATGGAAGTGATCCGTGGTATGGGCGTAGAAAAAACCGTTGGTGTCAAACCACCGGGAGGCGTGCGTACTGCGGAAGATGCGCAGAAAATTCCCGCCATCGCAGACGAACTATGTGGTGTGGACTGGGCAgcg                                                                                                                                                                  >  2:150351/1‑249 (MQ=255)
                                                                              cGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAACGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGCTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTAGGGGCTGACTGGGCAGATgcgc                                                                                                                                                              >  1:571436/1‑251 (MQ=255)
                                                                                          ccGTAAAGCGTCTGAAATCTCCATCAAAGCGGGCGCGGACTTCATCAAAACCTCTACCGGTAAAGCGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACAACCGCt                                                                                                                                                  >  1:526718/1‑251 (MQ=255)
                                                                                           cGTAAAGCGTCTTAAATCTACAGTAAAGCGGGTGCGGACTTCAACAAAACCTCTACCGGTAAATTGGCTGTGAACGAGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAAGATCTCGCCAGCGCAGATGAACTTTTCAGTGCTTACTGGGCAGATGCGCGTCACTACCGCtt                                                                                                                                                 <  2:593478/251‑1 (MQ=255)
                                                                                            gTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCttt                                                                                                                                                <  2:248797/251‑1 (MQ=255)
                                                                                             tAAAGCGTCTGAAATCTCCATCAAAGCGAGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGAAGGAAGTGATCCGTGATATGGTCGTAGAAAAAACCGATGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGAGCCGAAATATCTCGCCAACGCAGATGACCTTTTCGGTGCAGACTAGGCAGATACGCGGCACTCACGCtttt                                                                                                                                               >  2:24585/1‑250 (MQ=255)
                                                                                                aGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTggtg                                                                                                                                            >  1:281475/1‑249 (MQ=255)
                                                                                                 gCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTggtgc                                                                                                                                           <  1:483050/251‑4 (MQ=255)
                                                                                                       gAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCa                                                                                                                                      <  1:20271/250‑1 (MQ=255)
                                                                                                            ctcCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCTCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACt                                                                                                                                                                                                                                 <  2:14240/154‑1 (MQ=255)
                                                                                                            ctcCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCTCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACt                                                                                                                                                                                                                                 >  1:14240/1‑154 (MQ=255)
                                                                                                              ccATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCctg                                                                                                                                 <  1:644608/248‑1 (MQ=255)
                                                                                                                agtaaaGCGGGTGCGGACTTCATCAAAACCTCTGCCGGTAAAGTGGCTGTGAACGCGACGCCTGAAAGCGCGAGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAGCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAAGATCTCGCCATCGCAGATGAACTGTTCGGGGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGc                                                                                                                            <  2:103645/248‑1 (MQ=255)
                                                                                                                    aaGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGc                                                                                                                        <  1:394674/251‑1 (MQ=255)
                                                                                                                            gCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTTCTGGCAAGCCTGTTGa                                                                                                                 >  2:132326/1‑250 (MQ=255)
                                                                                                                              ggACTGCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGGGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCg                                                                                                                                                    <  2:79741/213‑1 (MQ=255)
                                                                                                                                     atcaAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTggg                                                                                                       >  1:400001/1‑251 (MQ=255)
                                                                                                                                       aaaacccTCTACCGGTATAGTGGCTGGTAACGCGACGCCGGAAAGCGCGCGCCTCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGTCGGCGCGCGTACTGCGGAAGATGCGCAGAAATAACTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGTAGACGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTgggcc                                                                                                     <  1:354609/246‑3 (MQ=255)
                                                                                                                                          gctcgtccgagctAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGAGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑Cg                                                                                                 <  1:375903/239‑1 (MQ=255)
                                                                                                                                                  ccGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGACGTGAACCGTGATATGGGCGTAGAAAAAACCGTTCGCTTCACACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGGAGACATATCTCGCCATCCCAGCTGAACTGATCGGTGCTGACAGGGCAGATGCCCGTAACTACCGCTTTCGTGCTTACAGCCTGCTGGCAAGCACGTTGACAGCGCTgggccccggtgccggg                                                                                          >  1:650303/1‑238 (MQ=255)
                                                                                                                                                     gTAAAGTGTCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAag                                                                                      <  1:58816/251‑1 (MQ=255)
                                                                                                                                                       aaaGTGGGTGTGAACGCGACGCCGGAAAGCGCGAGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGGTTAAAACCGGCGGGCGGCGTGCGTACTGCGGTAGATGAGCAGAAATATATCGCTATCGCAGATGAACTGTTCGGTGCGGACTGGGCAGATGCGTGTCACTTCCGCTTTGGTGCTTCCAGCCTGCGGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAagag                                                                                    <  2:193063/251‑1 (MQ=255)
                                                                                                                                                       aaaGAGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGCACTGTTCGGTGCTGACTGGGCAGATGCGCGGCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTgggccccgagagcggaagcgc                                                                                     >  2:375903/1‑233 (MQ=255)
                                                                                                                                                         aGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGTCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAgcg                                                                                  <  1:219610/251‑1 (MQ=255)
                                                                                                                                                            ggCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGcc                                                                                >  1:45776/1‑250 (MQ=255)
                                                                                                                                                            ggCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGGAGAAAAAACCGTTTGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTGCGGGGCTGAGTGGGCAGATGCGCGTCACTACCGCTTTGGGGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCca                                                                               <  1:276554/251‑1 (MQ=255)
                                                                                                                                                                           gCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACtaagtaa                                                                <  1:87709/251‑1 (MQ=255)
                                                                                                                                                                               gcaaGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGAGCTTCCAGCCTGCTGGCAAGCCTGTTGACAGCGCTGGGCCA‑CGGTGACGGTACGAACGCCAGCAGCTACTAAGTAAGATg                                                            >  2:430494/3‑251 (MQ=255)
                                                                                                                                                                               gAAAGCGCGCGCATCATGATGGCAGGGATCCGTGATCTGGGCGTAGAAACAACCGTTGGGTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCCCTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCT‑CGGTGCCGGTACGAGCACCCGCAGCTCCtaagaacaatg                                                            >  2:388835/1‑244 (MQ=255)
                                                                                                                                                                               gAAAGCGCGCGCATCATGATGGAAGTGATCCTTGATATGGGCGTAGAAAAAACCGTTGTTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAGGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTGCGGTTCTGACTGGTCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAATTAAGa                                                              <  2:557682/249‑1 (MQ=255)
                                                                                                                                                                                 aaGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCt                                                          <  1:523746/251‑1 (MQ=255)
                                                                                                                                                                                   gcgcgcgcATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTggtg                                                                                                                                            <  2:251935/168‑3 (MQ=255)
                                                                                                                                                                                   gcgcgcgcATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTggtg                                                                                                                                            >  1:251935/1‑166 (MQ=255)
                                                                                                                                                                                           atcatGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGa                                                <  1:430494/251‑1 (MQ=255)
                                                                                                                                                                                                gatgGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgc                                           <  1:71953/251‑1 (MQ=255)
                                                                                                                                                                                                 atgGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgct                                          >  1:466643/1‑251 (MQ=255)
                                                                                                                                                                                                   ggAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑Cggtga                                                                                             >  1:483439/1‑195 (MQ=255)
                                                                                                                                                                                                   ggAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTGGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCCCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTCCCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑Cggtga                                                                                             <  2:483439/198‑4 (MQ=255)
                                                                                                                                                                                                       gTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTc                                                                                                                                                                                    >  1:472100/1‑108 (MQ=255)
                                                                                                                                                                                                       gTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGAGGCTCAGAAATATCTCGCCATCGCAGATGAACTGTTc                                                                                                                                                                                    <  2:472100/108‑1 (MQ=255)
                                                                                                                                                                                                          aTCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATc                                 >  1:158222/1‑251 (MQ=255)
                                                                                                                                                                                                             cGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAgg                              >  1:445083/1‑251 (MQ=255)
                                                                                                                                                                                                              gTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGCAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAACGCTGGGCCA‑CGGTGACGGTAAGCGCGCCAGCCGCTACTCAGGAAGATGCTTTACGCCTGATCCGCTGCGCTTATCAGGc                             >  2:5349/1‑251 (MQ=255)
                                                                                                                                                                                                               ggATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGcc                            <  1:16165/250‑1 (MQ=255)
                                                                                                                                                                                                                gATATGGGCGTAGAAAAAACAGTTGGATTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGAGGCGAGTCACTACCGCTTTGGTGCTTACAACCTGATGGCAAGCCTGTTGAATGCGCTGG‑TCAACGGTGACGGTAAGAGCGCAAGCAGCTACTAAGGAAGTTGCTTTACGGCTGATGAGCTGCGCTTATCacgctt                           >  2:116612/1‑246 (MQ=255)
                                                                                                                                                                                                                       gCGTAGAAAAAACCGTGGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGAc                    <  1:388835/251‑1 (MQ=255)
                                                                                                                                                                                                                           aGAAAAAACCGTTGTTTTCAAACCGCCGGGAGGCGTGCGTACGGCGGAAGATGCGCAGAAATATCTCGCCAGTGCAGATGATCTGGTCGGTGCTGACTGGGCAGATACGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCA‑CGGCGACGGTAAGAGCGCCAGCAGCTGCTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTAt                <  2:16722/251‑1 (MQ=255)
                                                                                                                                                                                                                           aGAAAAAACCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTa                 <  1:132326/250‑1 (MQ=255)
                                                                                                                                                                                                                              aaaaaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCt              >  1:361568/1‑250 (MQ=255)
                                                                                                                                                                                                                                 aaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTAccc          <  1:495002/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aaCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTAccc          >  1:649923/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aCCGTTGGTTTCAAACCAGCGGGCGGCGTGCGTCCTGCGGAAGATGCGCAGAAATATCGCGCCATCGCAGATTAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCg         <  2:571436/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gTTGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTAGCTACCCGTAg      <  2:649923/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ttGGTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGGTGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTCCCCGTAgg     >  1:416489/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gTTTCAAACCAGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATCGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGTGCTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCA‑CGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCg  <  1:213778/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
GGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAAACCGGCGAACTGAAAGACGAAGCGCTGATCCGTAAAGCGTCTGAAATCTCCATCAAAGCGGGTGCGGACTTCATCAAAACCTCTACCGGTAAAGTGGCTGTGAACGCGACGCCGGAAAGCGCGCGCATCATGATGGAAGTGATCCGTGATATGGGCGTAGAAAAAACCGTTGGTTTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCA‑CGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCG  >  NC_000913/4617685‑4618168

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: ‑