Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,587,589 A→G 14.0% E214G (GAG→GGG)  mrr → methylated adenine and cytosine restriction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,587,5890AG14.0% 185.9 / 3.1 71E214G (GAG→GGG) mrrmethylated adenine and cytosine restriction protein
Reads supporting (aligned to +/- strand):  ref base A (23/37);  new base G (9/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTT  >  NC_000913/4587353‑4587826
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
tGCAAAGAGCAGCCCGGGCGATCGGTTAAATCAGGCGGTAAAAGAGCTTCGTGATGCAGTGGCTGATGAGGTTCTGTAAAACTTATTGCATGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGgaga                                                                                                                                                                                                                                                 <  1:462349/239‑1 (MQ=255)
tGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTTGCGGGACTTGAGATGGGGGGATCGCTGGTGTGGTATCCCTTTGTAAACTTGGCCTGggga                                                                                                                                                                                                                                                 >  1:158948/1‑239 (MQ=255)
            cccGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCCCCGGGATGGTTTTCAGCGTGTTGGGGGGACTGGGGATGGGGGGATCCCTGGGGTGATAACGCTGGTAAAAAAAGGCCTGGGGAAAAATTATGttt                                                                                                                                                                                                                                     >  2:319648/1‑238 (MQ=255)
             ctggACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGtt                                                                                                                                                                                                                                      <  1:560156/235‑1 (MQ=255)
                     tCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTTCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGGTAAAATTGGCCTGGAGAAAGGTTATGTTTAGGCAAAAc                                                                                                                                                                                                                            >  1:128861/1‑239 (MQ=255)
                      cGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGGTGGGGTGGTATCGCTTGGTAAACTTGGCCTGGGGAAAGTTTTTGTTGTGGGAAAACg                                                                                                                                                                                                                           >  1:393618/1‑239 (MQ=255)
                       ggTTTGAGCTGGCGCTAAAAGAGCTTCGTGACGCGGTGGCTGAGGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGGTTTGAAGGCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGt                                                                                                                                                                                                                          <  2:421022/239‑1 (MQ=255)
                       ggTTAGATCAGGCGCTAAAAGGGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAATCTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGt                                                                                                                                                                                                                          <  1:913755/239‑1 (MQ=255)
                        gttaagtCAGGCGCTAAATGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGGCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCGACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTTAGGCAAAACGt                                                                                                                                                                                                                          <  1:227767/232‑1 (MQ=255)
                            gATCAGGCGCTAAAAGTGCTTCGTGATGCGGTGGTTTAGTGGGTTCTGGAAATCTTATTGCAGGGTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCa                                                                                                                                                                                                                     <  1:435386/239‑1 (MQ=255)
                              gtaGGGGCTAAAAGAGGTTCGGGTTTCGGTGGCTGGTGAGGTTCTGGAAAACTTATGGCAGGTTTCTGCTTGGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGGGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                                                                                   <  2:667103/237‑1 (MQ=255)
                               cAGGCGCTAAAAGGTGTTCGTTATGCGGTGGCTGATGAGGTTCTGGAAAACTTAGTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                                                                                   <  1:797719/238‑1 (MQ=255)
                               cAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCACCGCCTGGGGAATGGCGGCCACCGTGTGGATTTGCAGGCTCTTGGCGGTACTGGGGACGGGGGCATGAAGGGGGGGATCTGGCGTGATAAACTTGCACTGGAGGAAGTGTTGGTTGGGGCAAAAAATTGGCGGaa                                                                                                                                                                                                                  >  2:762077/1‑239 (MQ=255)
                                       aaaaGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTTGAGATGGTGGCATCGATGGGGTGATATCGCTTGTTAAACTTGGCCTGGGGAAAGTTTTTGTTCAGGGAAAACTGTGGCAGAATACTGTag                                                                                                                                                                                                          >  1:225937/1‑239 (MQ=255)
                                        aaaGAGCTTCGTGATGCGGTGGCTGATGAGGTTCAGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGGATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTagg                                                                                                                                                                                                         <  1:103904/239‑1 (MQ=255)
                                                tCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGCGTGTTGGGGGTACTGGGGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTTGCCTGGAGAAAGTTTATGGTCAGGCAAAACGTTGGCGGAATACTGTAGGCAGGCCAg                                                                                                                                                                                                 >  1:182732/1‑239 (MQ=255)
                                                      tggggtgggtgGTGGGGTTTTGGAAAATTTTTTGCAGGTTTCTTCTTCGGGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑c                                                                                                                                                                                          <  2:771922/236‑1 (MQ=255)
                                                           tGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGGGGCATCGATGGGGTGATATCGCTTGGTAAACTTGGCCTGGGGAAAGGTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCa                                                                                                                                                                                     >  1:82110/1‑239 (MQ=255)
                                                             gCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTTTTTGCGGTACTGGAGATGGTGGACTCGCAGGTGGGGTTTCCCTTGTTTAAATTGTGCTTGGGGAAGATTAGTTTTCGGCAAAAAGTGGGACGAAAAAGTTAGGGGGGCCAGAAATA‑AAAGCAtt                                                                                                                                                                                   >  2:685232/1‑239 (MQ=255)
                                                                  tggggTTCTGGAAAACTTATTGGAGGTTTCCCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACg                                                                                                                                                                              <  1:319648/236‑1 (MQ=255)
                                                                  tgaGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTAc                                                                                                                                                                               <  1:583342/238‑1 (MQ=255)
                                                                                 cTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTTCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGTGTGATATCGCTTGATAAAATTGGCCTGGGGAAAGTTTATGTTCAGGCAAAAAGTTGGCCGAATACTGTAGGCAGGCCCGAATTA‑CAGGCATTTTACGGGGCC‑CCGGCTGGGc                                                                                                                                                              >  1:152139/1‑239 (MQ=255)
                                                                                           gTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGGTATCGCCTGATAAACTTTGCCTGGGGAAAGGTTATGTTCCGGCACAACGTTTGCAGAATACTGTAGGCAGGCCCAAATTA‑CAGGCAATTTAAGGCGCC‑CTGGCTTGGGAAAAA‑AC‑Gaa                                                                                                                                                  >  1:892692/1‑239 (MQ=255)
                                                                                            ttttcccTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑Gaa                                                                                                                                                  <  1:519468/233‑1 (MQ=255)
                                                                                                cccTTCGGGTTTTGTAGTCTTTTTTTTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCCCCGTGATGATTTGCAGCGTGGTGGCGGTTCTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTg                                                                                                                                             <  2:393618/238‑1 (MQ=255)
                                                                                                     cgcgCTTTGAAGTCCTTGTTCTGGATGTTTTGCATCGGCTGGGGTAAGGAGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCCTCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGgtg                                                                                                                                        <  1:717542/239‑1 (MQ=255)
                                                                                                             gAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTAc                                                                                                                                <  1:59638/239‑1 (MQ=255)
                                                                                                                           ggATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTTAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACCGCGCA‑CCTGCTGGGGAAAAAGGC‑GAAA‑GTGGGGTGGTTAATACCAATTTCGGGTTTa                                                                                                                  >  1:714162/1‑239 (MQ=255)
                                                                                                                            gATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGGGGGATATCGCTTGGAAAACTTGGCCTGGGGAAAGTTTTTGTTTATGCAAAACGTTGGCAGAAAACTGTATGAAGGCCAGAATTT‑AAGGGAATTTTATGCGCACCTGTCGGGCAAAAAA‑GG‑GAAACAGGGGGGGTTATTTACCCCTCTTGGTTTTa                                                                                                                  >  2:777436/1‑239 (MQ=255)
                                                                                                                                 tttGCATCGCCTGGGGTATGGCGGCCACCGTGGTGATTTGCAGCGTGTTGGCGGTACTGGGGGTGGTGGCATCGCTGGGGGGGTATCTCTTGTGAAAATTTGCCTGGGGAAAGGTTTTGTTCCGGGAAAACTGTGGCCGCATACTGTAGGCCGGACAGCATTACCAGGC‑TTTTACGGCGGC‑CGACACGGGCAAAAA‑GCAGAAAAGGGGGGGGTTTATAACCACTTTTgttttattttta                                                                                                             >  2:470632/1‑228 (MQ=255)
                                                                                                                                             ggggTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGTCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACt                                                                                                <  1:794784/239‑1 (MQ=255)
                                                                                                                                                    ggcggcCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGGGTGGTATCGCTTGATAAAATTTGCCTGGAGAAAGTTTATGTTTAGGCAAAACGTTGGCAGAAAACTGTAGGCCGGCCAGACTTA‑CAGGCATTTTTCGGCGCA‑CTTGCCGGGGAAAAA‑GC‑GGAACATGGGGGGTTTTTTACCCCTTTTGGATTTACTTTTCAGGCGCGGGACCTTGCCCa                                                                                         >  2:783477/1‑239 (MQ=255)
                                                                                                                                                      cggcCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAGGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCaa                                                                                        <  1:648243/238‑1 (MQ=255)
                                                                                                                                                          cACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGt                                                                                   <  1:470632/239‑1 (MQ=255)
                                                                                                                                                           aCCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATAGCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGt                                                                                   <  1:184541/238‑1 (MQ=255)
                                                                                                                                                             cGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGa                                                                                <  1:186069/239‑1 (MQ=255)
                                                                                                                                                                   gttttGCCGCGTGGTGGGGGTTCTGGAGAGGGTGGCATCGGTTGGGTGATAACGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGTA‑CTGGCTGGTCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTa                                                                           <  1:507668/236‑1 (MQ=255)
                                                                                                                                                                       tGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGgtg                                                                      <  1:783477/239‑1 (MQ=255)
                                                                                                                                                                       tGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGGGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGCCGAGGGTATgggg                                                                      >  1:354448/1‑237 (MQ=255)
                                                                                                                                                                         cAGCGTGTTGGCGGTCCTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGCATACTGTAGGCAGGCCAGATTTA‑CAGGCATTTTACGGCGCC‑CTGGCTGGGCAAAAA‑AC‑GAAACGTGGGGGGGTTATTACCACTCTTGGATTTAATTCTCAGGCGCGGGACTTTTCCCACACCGCCGCGGGTATGgttgt                                                                    >  2:836302/1‑236 (MQ=255)
                                                                                                                                                                             gttttGGGGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTgttggt                                                                 <  2:619747/235‑1 (MQ=255)
                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGGGTTTATTACCACTTCTGGGTTTACTTTTCAGGCGCGTGCCTTTGCCCAAACCGCCGAGGGTGTGGGgtgggttg                                                               >  2:552278/1‑239 (MQ=255)
                                                                                                                                                                               gtTGGCGGTACTGGAGAGTGTGGCATCGATGGTGTGATAGCGCTTGAGTAACTTGGCCGGGAGAAAGGTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGa                                                              <  1:739068/239‑1 (MQ=255)
                                                                                                                                                                                gttgCGGTACGGGTGTTGGGGGCATCGATGGTGTGATTTCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAGAACGTTGGCAGAATTCTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTTGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGAt                                                             <  2:182732/236‑1 (MQ=255)
                                                                                                                                                                                     ggTATTGGGGATGGGGGGATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                         <  2:152139/238‑1 (MQ=255)
                                                                                                                                                                                     ggTACGGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTTGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                         <  1:762077/238‑1 (MQ=255)
                                                                                                                                                                                      gttcTGGAGATGGGGGAATCGATTGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGa                                                        <  2:128861/235‑1 (MQ=255)
                                                                                                                                                                                        aagggAGATGGTCGCATCGATGGTTCGATATCGTTTGATAAACTTGGCCTGGGGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCGCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAc                                                      <  1:735418/235‑1 (MQ=255)
                                                                                                                                                                                         cTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTggg                                                                                                                                           <  1:301937/152‑1 (MQ=255)
                                                                                                                                                                                         cTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTggg                                                                                                                                           >  2:301937/1‑152 (MQ=255)
                                                                                                                                                                                             agaTGGTGGCATTGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCg                                                                                 <  2:56524/206‑1 (MQ=255)
                                                                                                                                                                                             agaTGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCg                                                                                 >  1:56524/1‑206 (MQ=255)
                                                                                                                                                                                                tggtggGATCGATGGTGTGATATCGGTTGAAAAACTTGGCCTGGAGAAAGTTTAAGTTCAGGCAAAACGGTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  2:212083/238‑1 (MQ=255)
                                                                                                                                                                                                tggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  1:552278/238‑1 (MQ=255)
                                                                                                                                                                                                tggtggCATCGATGGTGTGATATCCCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  1:304931/238‑1 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTc                                                                                  >  2:744421/1‑201 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTc                                                                                  <  1:744421/201‑1 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGGGACTTTGCCCAATCCGTCGCGGGAATGGTGTTGGTTGATGGGGAACGCCTGGGGCa                                            >  1:198335/1‑239 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGGGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGc                                             >  1:417280/1‑238 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGAGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGGGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAGGCCGGGTg                                              >  1:450720/1‑237 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGCGAGGGGTTTATTACCCCTTCTGGATTTCCTTCTCAGGCGCGTGACTTTGCCCAATCCGCCGAGGGTATGGGGTTGGTTGATGGGGGACGCCGGGGGc                                             >  2:703988/1‑238 (MQ=255)
                                                                                                                                                                                                               ttggATATCGCTTGATAAACTTGGCCGGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTaa                                       <  1:843142/227‑1 (MQ=255)
                                                                                                                                                                                                               tgtgATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGCATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGaaa                              >  1:862724/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGGAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCAGTTTACGGCCCC‑CTGGCTGGGCAAAAA‑GC‑GAACCGTGGGGTGTTTACTACCACTTCTGGATTTACTTCTCAGGAGCGTCACTTTGCCCCATCCGTCGGGGGTATGGCGTTTGTGACAGGGGAAACCCTGGTGCCCCTACAGAACCAAAACGAaggagg                    >  2:37168/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCCTGCCAGAATTA‑CAGGCATTTTACGGCTCA‑CTGGCTGGGCAAAAA‑TC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTTCTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGGCGGTCTGTGTTTGGTTGGTGGGGCACGCCTGGTGCACTTAATGATAGAAAACGAAGTAgg                    >  2:776788/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGGGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAgg                    >  1:301741/1‑239 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTgggg                                                                                                                                          <  2:606338/117‑1 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTgggg                                                                                                                                          >  1:606338/1‑117 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGtt                >  1:687926/1‑239 (MQ=255)
                                                                                                                                                                                                                                    tGGCCTGGGGAAATTTTTTTTTCAGGTAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCAc         <  2:354448/239‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCGGGAGAAAGGTTAAGTTCGGGCAGAACGTTGGCAGAATCCTGTAGGCAGGCCAGGATTA‑CAGGCATTTTGCGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGGGGGGTGTTTATTACCTGTTCTGGATTTACTTCTCAGGTGCGTGACTTTGCGCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGt       <  1:37168/239‑1 (MQ=255)
                                                                                                                                                                                                                                       ccTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGtt      >  1:140634/1‑239 (MQ=255)
                                                                                                                                                                                                                                       ccTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGGTTTGGTTGATGGGGAACGCCTGGTGCCCTTAATGATCGAAAACGAAAAAGGGGTGTTTTCAAGtt      >  2:710044/1‑239 (MQ=255)
                                                                                                                                                                                                                                           gagaAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATTGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGCAAACGAACTAGGGGTTTCTTCCCGTttgtt  >  1:226147/1‑239 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                 
TGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTT  >  NC_000913/4587353‑4587826

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