Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,587,590 A→G 14.8% 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,5900AG14.8% 200.3 / 5.4 81E214G (GAG→GGG) mrrmethylated adenine and cytosine restriction protein
Reads supporting (aligned to +/- strand):  ref base A (29/39);  new base G (11/2);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTT  >  NC_000913/4587354‑4587827
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
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)
                                  gcgcTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTTCAGCGTGTTGGCGGGACTGGGGGGGGGGGCATCGATGGTGTGGTATCGTTGTGTAAACTTGGGCCGGGAAAAATTTATGtttggagaaaaaggtgggagaaaaa                                                                                                                                                                                                                 >  2:186069/1‑216 (MQ=255)
                                       aaaaGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTTGAGATGGTGGCATCGATGGGGTGATATCGCTTGTTAAACTTGGCCTGGGGAAAGTTTTTGTTCAGGGAAAACTGTGGCAGAATACTGTag                                                                                                                                                                                                            >  1:225937/1‑239 (MQ=255)
                                        aaaGAGCTTCGTGATGCGGTGGCTGATGAGGTTCAGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGGATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTagg                                                                                                                                                                                                           <  1:103904/239‑1 (MQ=255)
                                                tCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCCCCGTGATGATTTGCAGCGTGTTGGCGGGACTGGAGATGGGTGCATCCATGGTGTGATATCGCTTTGTAAAACTGGCCGGGGGGAAGTTTATTTTAAGGCAAAACGTTggtggaaaaattgaggcacgcaca                                                                                                                                                                                                   >  2:59638/1‑217 (MQ=255)
                                                tCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGCGTGTTGGGGGTACTGGGGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTTGCCTGGAGAAAGTTTATGGTCAGGCAAAACGTTGGCGGAATACTGTAGGC‑AGGCCAg                                                                                                                                                                                                  >  1:182732/1‑239 (MQ=255)
                                                      tggggtgggtgGTGGGGTTTTGGAAAATTTTTTGCAGGTTTCTTCTTCGGGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑c                                                                                                                                                                                           <  2:771922/236‑1 (MQ=255)
                                                           tGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGGGGCATCGATGGGGTGATATCGCTTGGTAAACTTGGCCTGGGGAAAGGTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCa                                                                                                                                                                                      >  1:82110/1‑239 (MQ=255)
                                                             gCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTTTTTGCGGTACTGGAGATGGTGGACTCGCAGGTGGGGTTTCCCTTGTTTAAATTGTGCTTGGGGAAGATTAGTTTTCGGCAAAAAGTGGGACGAAAAAGTTAGGG‑GGGCCAGAAATA‑AAAGCAtt                                                                                                                                                                                    >  2:685232/1‑239 (MQ=255)
                                                                  tggggTTCTGGAAAACTTATTGGAGGTTTCCCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACg                                                                                                                                                                               <  1:319648/236‑1 (MQ=255)
                                                                  tgaGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTAc                                                                                                                                                                                <  1:583342/238‑1 (MQ=255)
                                                                                 cTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTTCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGTGTGATATCGCTTGATAAAATTGGCCTGGGGAAAGTTTATGTTCAGGCAAAAAGTTGGCCGAATACTGTAGGC‑AGGCCCGAATTA‑CAGGCATTTTACGGGGCC‑CCGGCTGGGc                                                                                                                                                               >  1:152139/1‑239 (MQ=255)
                                                                                           gTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGGTATCGCCTGATAAACTTTGCCTGGGGAAAGGTTATGTTCCGGCACAACGTTTGCAGAATACTGTAGGC‑AGGCCCAAATTA‑CAGGCAATTTAAGGCGCC‑CTGGCTTGGGAAAAA‑AC‑Gaa                                                                                                                                                   >  1:892692/1‑239 (MQ=255)
                                                                                            ttttcccTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑Gaa                                                                                                                                                   <  1:519468/233‑1 (MQ=255)
                                                                                              tgttctttcgggtttGTAGTCATTGTTCTGGTCGTTGTTCAGCGGCGGGTGTATGGCGGCCAGCGGTATGATTGTCAGCGGGTTGGGGGTACTGGAGATGTGGGCCTCGATGGTGTGATGTGGCTTGATAAACTTGGCCTGGAGAAATTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑Gaaaag                                                                                                                                                <  2:144648/227‑3 (MQ=255)
                                                                                                cccTTCGGGTTTTGTAGTCTTTTTTTTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCCCCGTGATGATTTGCAGCGTGGTGGCGGTTCTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTg                                                                                                                                              <  2:393618/238‑1 (MQ=255)
                                                                                                     cgcgCTTTGAAGTCCTTGTTCTGGATGTTTTGCATCGGCTGGGGTAAGGAGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCCTCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGgtg                                                                                                                                         <  1:717542/239‑1 (MQ=255)
                                                                                                             gAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTAc                                                                                                                                 <  1:59638/239‑1 (MQ=255)
                                                                                                                           ggATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTTAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACCGCGCA‑CCTGCTGGGGAAAAAGGC‑GAAA‑GTGGGGTGGTTAATACCAATTTCGGGTTTa                                                                                                                   >  1:714162/1‑239 (MQ=255)
                                                                                                                            gATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGGGGGATATCGCTTGGAAAACTTGGCCTGGGGAAAGTTTTTGTTTATGCAAAACGTTGGCAGAAAACTGTATGA‑AGGCCAGAATTT‑AAGGGAATTTTATGCGCACCTGTCGGGCAAAAAA‑GG‑GAAACAGGGGGGGTTATTTACCCCTCTTGGTTTTa                                                                                                                   >  2:777436/1‑239 (MQ=255)
                                                                                                                                 tttGCATCGCCTGGGGTATGGCGGCCACCGTGGTGATTTGCAGCGTGTTGGCGGTACTGGGGGTGGTGGCATCGCTGGGGGGGTATCTCTTGTGAAAATTTGCCTGGGGAAAGGTTTTGTTCCGGGAAAACTGTGGCCGCATACTGTAGGC‑CGGACAGCATTACCAGGC‑TTTTACGGCGGC‑CGACACGGGCAAAAA‑GCAGAAAAGGGGGGGGTTTATAACCACTTTTgttttattttta                                                                                                              >  2:470632/1‑228 (MQ=255)
                                                                                                                                             ggggTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGTCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACt                                                                                                <  1:794784/239‑1 (MQ=255)
                                                                                                                                                   tgtcggcCACCGTGATGATTTGCAGCGTGTTGGCGGCCCTGGGGACGGTGGACGCCATGGGGGGATATCGCTTGATTAAATTGGTCTGGATTAATTTTATTGTCAGGCACAAAGTGTGCTGAATAATGTATGCGAGGCC‑GAATTA‑AAGGGCATTTAAGGCCGA‑CTGGCTGGGGCAAAA‑GG‑TGAACGGGGGGTGTTTTTTTACACATTTGGGTTTACTTTCTCAGGcgcggaattttccc                                                                                           >  2:648243/4‑229 (MQ=255)
                                                                                                                                                    ggcggcCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGGGTGGTATCGCTTGATAAAATTTGCCTGGAGAAAGTTTATGTTTAGGCAAAACGTTGGCAGAAAACTGTAGGC‑CGGCCAGACTTA‑CAGGCATTTTTCGGCGCA‑CTTGCCGGGGAAAAA‑GC‑GGAACATGGGGGGTTTTTTACCCCTTTTGGATTTAC‑TTTTCAGGCGCGGGACCTTGCCCa                                                                                         >  2:783477/1‑239 (MQ=255)
                                                                                                                                                      cggcCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAGGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCaa                                                                                        <  1:648243/238‑1 (MQ=255)
                                                                                                                                                          cACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGt                                                                                   <  1:470632/239‑1 (MQ=255)
                                                                                                                                                           aCCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATAGCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGt                                                                                   <  1:184541/238‑1 (MQ=255)
                                                                                                                                                             cGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGa                                                                                <  1:186069/239‑1 (MQ=255)
                                                                                                                                                                   gttttGCCGCGTGGTGGGGGTTCTGGAGAGGGTGGCATCGGTTGGGTGATAACGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGTA‑CTGGCTGGTCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTa                                                                           <  1:507668/236‑1 (MQ=255)
                                                                                                                                                                      ttGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTTTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCTGGCAAAACGTTGGCAGAATACTGTTGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACATGGGGGGTTTATTACCACTATTGGATTTAT‑TTCTCCGGCGCGTGAATTGTCCCCATcccctcgggtttagga                                                                       >  2:507668/1‑225 (MQ=255)
                                                                                                                                                                       tGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGGTTGATAAACTTGGCCTGGAGAAAGGTTATGTTCAGGCCAAACGTTGGTAGAATACTGTAGGC‑AGGCCTGAATTA‑CCGGCCTTTTACGGCCAA‑CTGGCTTGGGAACAA‑GG‑GGACAGTGGGGTGTTTATTAACCCTTCGGGATTTAC‑TTTTACGGCGCGTGAATTTggcagacaagtcctgggttggggg                                                                      >  2:663656/1‑216 (MQ=255)
                                                                                                                                                                       tGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGgtg                                                                      <  1:783477/239‑1 (MQ=255)
                                                                                                                                                                       tGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGGGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGCCGAGGGTATgggg                                                                      >  1:354448/1‑237 (MQ=255)
                                                                                                                                                                         cAGCGTGTTGGCGGTCCTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGCATACTGTAGGC‑AGGCCAGATTTA‑CAGGCATTTTACGGCGCC‑CTGGCTGGGCAAAAA‑AC‑GAAACGTGGGGGGGTTATTACCACTCTTGGATTTAA‑TTCTCAGGCGCGGGACTTTTCCCACACCGCCGCGGGTATGgttgt                                                                    >  2:836302/1‑236 (MQ=255)
                                                                                                                                                                             gttttGGGGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTgttggt                                                                 <  2:619747/235‑1 (MQ=255)
                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTGAGGCAAGACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGTAAAAA‑GG‑GAAAAGTGGGGGGTTTATTCGCCATTTTGGCTTTCC‑TTCTAAGCGGCGGGACTTTGCCCAATCCGTGGAggggaggggatgtgttg                                                               >  2:304931/1‑225 (MQ=255)
                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGGGTTTATTACCACTTCTGGGTTTAC‑TTTTCAGGCGCGTGCCTTTGCCCAAACCGCCGAGGGTGTGGGgtgggttg                                                               >  2:552278/1‑239 (MQ=255)
                                                                                                                                                                               gtTGGCGGTACTGGAGAGTGTGGCATCGATGGTGTGATAGCGCTTGAGTAACTTGGCCGGGAGAAAGGTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGa                                                              <  1:739068/239‑1 (MQ=255)
                                                                                                                                                                                gttgCGGTACGGGTGTTGGGGGCATCGATGGTGTGATTTCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAGAACGTTGGCAGAATTCTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTTGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGAt                                                             <  2:182732/236‑1 (MQ=255)
                                                                                                                                                                                    cgtttcttgggttggtggCATTGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                         <  2:225937/227‑1 (MQ=255)
                                                                                                                                                                                     ggTATTGGGGATGGGGGGATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                         <  2:152139/238‑1 (MQ=255)
                                                                                                                                                                                     ggTACGGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTTGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                         <  1:762077/238‑1 (MQ=255)
                                                                                                                                                                                      gttcTGGAGATGGGGGAATCGATTGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGa                                                        <  2:128861/235‑1 (MQ=255)
                                                                                                                                                                                        aagggAGATGGTCGCATCGATGGTTCGATATCGTTTGATAAACTTGGCCTGGGGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCGCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAc                                                      <  1:735418/235‑1 (MQ=255)
                                                                                                                                                                                         cTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTggg                                                                                                                                            <  1:301937/152‑1 (MQ=255)
                                                                                                                                                                                         cTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTggg                                                                                                                                            >  2:301937/1‑152 (MQ=255)
                                                                                                                                                                                             agaTGGTGGCATTGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCg                                                                                 <  2:56524/206‑1 (MQ=255)
                                                                                                                                                                                             agaTGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCg                                                                                 >  1:56524/1‑206 (MQ=255)
                                                                                                                                                                                                tggtggGATCGATGGTGTGATATCGGTTGAAAAACTTGGCCTGGAGAAAGTTTAAGTTCAGGCAAAACGGTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  2:212083/238‑1 (MQ=255)
                                                                                                                                                                                                tggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  1:552278/238‑1 (MQ=255)
                                                                                                                                                                                                tggtggCATCGATGGTGTGATATCCCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  1:304931/238‑1 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTc                                                                                  <  1:744421/201‑1 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTc                                                                                  >  2:744421/1‑201 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGGGACTTTGCCCAATCCGTCGCGGGAATGGTGTTGGTTGATGGGGAACGCCTGGGGCa                                            >  1:198335/1‑239 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGGGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGc                                             >  1:417280/1‑238 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGAGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGGGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAGGCCGGGTg                                              >  1:450720/1‑237 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGCGAGGGGTTTATTACCCCTTCTGGATTTCC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGCCGAGGGTATGGGGTTGGTTGATGGGGGACGCCGGGGGc                                             >  2:703988/1‑238 (MQ=255)
                                                                                                                                                                                                      cATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCAAATCCGCCGGGGCGTTGGTGTTGGTTGATGGGGacagccgggggcccctg                                        >  2:735418/1‑222 (MQ=255)
                                                                                                                                                                                                               ttggATATCGCTTGATAAACTTGGCCGGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTaa                                       <  1:843142/227‑1 (MQ=255)
                                                                                                                                                                                                               tgtgATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGCATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGaaa                              >  1:862724/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGGAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCAGTTTACGGCCCC‑CTGGCTGGGCAAAAA‑GC‑GAACCGTGGGGTGTTTACTACCACTTCTGGATTTAC‑TTCTCAGGAGCGTCACTTTGCCCCATCCGTCGGGGGTATGGCGTTTGTGACAGGGGAAACCCTGGTGCCCCTACAGAACCAAAACGAaggagg                    >  2:37168/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑CTGCCAGAATTA‑CAGGCATTTTACGGCTCA‑CTGGCTGGGCAAAAA‑TC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTTC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGGCGGTCTGTGTTTGGTTGGTGGGGCACGCCTGGTGCACTTAATGATAGAAAACGAAGTAgg                    >  2:776788/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGGGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAgg                    >  1:301741/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTAAGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCAATTCTGGATTTAC‑TTCTCAGGCCCGTGACTTTTCCCAATCCGTTGGGGGAATGGGTATGGTGGATCGGGAACGCCTGGCGCCCGTaaacatgcaaaaccttgtcct                    >  2:843142/1‑220 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTgggg                                                                                                                                           <  2:606338/117‑1 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTgggg                                                                                                                                           >  1:606338/1‑117 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGtt                >  1:687926/1‑239 (MQ=255)
                                                                                                                                                                                                                                    tGGCCTGGGGAAATTTTTTTTTCAGGTAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCAc         <  2:354448/239‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCGGGAGAAAGGTTAAGTTCGGGCAGAACGTTGGCAGAATCCTGTAGGC‑AGGCCAGGATTA‑CAGGCATTTTGCGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGGGGGGTGTTTATTACCTGTTCTGGATTTAC‑TTCTCAGGTGCGTGACTTTGCGCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGt       <  1:37168/239‑1 (MQ=255)
                                                                                                                                                                                                                                       ccTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGtt      >  1:140634/1‑239 (MQ=255)
                                                                                                                                                                                                                                       ccTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGGTTTGGTTGATGGGGAACGCCTGGTGCCCTTAATGATCGAAAACGAAAAAGGGGTGTTTTCAAGtt      >  2:710044/1‑239 (MQ=255)
                                                                                                                                                                                                                                           gagaAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATTGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGCAAACGAACTAGGGGTTTCTTCCCGTttgtt  >  1:226147/1‑239 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
TGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGC‑AGGCCAGAATTA‑CAGGCATTTTACGGCGCA‑CTGGCTGGGCAAAAA‑GC‑GAAACGTGGGGTGTTTATTACCACTTCTGGATTTAC‑TTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTT  >  NC_000913/4587354‑4587827

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