Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,587,543 A→G 13.0% G198G (GGA→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,5430AG13.0% 211.6 / 2.6 84G198G (GGA→GGG) mrrmethylated adenine and cytosine restriction protein
Reads supporting (aligned to +/- strand):  ref base A (13/58);  new base G (8/5);  total (21/63)
Fisher's exact test for biased strand distribution p-value = 2.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTGTTGGTTGATGGGG  >  NC_000913/4587308‑4587772
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
ccctCGTTTTAGATCCGAAACCCGTCTTCTCCGGTTATGAAGAACTTGCAAAGAGCAGCCCGGACGGTCGGTTAGAACAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGAt                                                                                                                                                                                                                                         <  1:547755/236‑1 (MQ=255)
    cGATTTAGAGCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCACGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGAtggt                                                                                                                                                                                                                                      <  2:966230/238‑1 (MQ=255)
           gATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  2:560915/239‑1 (MQ=255)
             tCCCAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGCCGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  2:6516/237‑1 (MQ=255)
             gccGAAAGCCGACTCGCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  2:887284/236‑1 (MQ=255)
                 aaaGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGACGTCATTGTTCTGGCTG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGTTGATTTGCATC‑GTGTTGGCGGTACTGGGGAGGGGGGCACCGgtggtgg                                                                                                                                                                                                                        >  2:1297053/1‑238 (MQ=255)
                 aaaGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCCCCTTGTTGATTTGCAGC‑GGTTTTGCGGTACTGTAGATGGTGGGACCGCTggggg                                                                                                                                                                                                                        >  2:550188/1‑236 (MQ=255)
                       gaatctcCCGAACATGAAGAACTTGCAAAGGGCCGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATc                                                                                                                                                                                                                  <  2:781587/236‑1 (MQ=255)
                         ttcgtccgactctGAAGAACGTGAAAAGGGCCGGCCGGGAGATCGGTTAGATCAGGCCGTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCTGGATTCTGCTTCGCGCTTTGAAGTCATTGTTGTGGATG‑TTTTTCTTCGCCTGGGGTATGGCGGCCAACGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGc                                                                                                                                                                                                                <  1:1100579/227‑1 (MQ=255)
                         ctctCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCCCCGGGATGATTTTCAGC‑GGGTTGGCGGTACTGGGGATGGTGGCATCGCTGGGGTGaattc                                                                                                                                                                                                                  >  2:1116275/1‑233 (MQ=255)
                         accccccggcgaTGGAGGACTTGCAAGAGGTAGCCCGGAGCATCTGTGAGATCAGGGTGTAAATGGTCTCCGTGCGGGTGGGGCTGATGATGTTCTGGAAAACTTATGTCAGGTTTCTCCTTTGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCGCCGTGATGATTTGGAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGc                                                                                                                                                                                                                <  1:729183/228‑1 (MQ=255)
                             cccGATCATGAAGAACTTGCAAAGGGCCGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTTATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGa                                                                                                                                                                                                            <  1:283800/239‑1 (MQ=255)
                                  tCATTAAGAACTTGCAAAGAGCAGCCCGGACGATCGGGTTGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACGTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAAc                                                                                                                                                                                                       <  1:845614/239‑1 (MQ=255)
                                      gaagaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGGGATGGTGGCCTCGATGGTGTGATATCGCTTGATAAACttt                                                                                                                                                                                                    >  2:454329/1‑237 (MQ=255)
                                       aagtaCTTGCAAAGGGCCGTCCGGACGATGGGTTAGATCGGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGc                                                                                                                                                                                                  <  1:732877/239‑1 (MQ=255)
                                       aaggaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGGTAAAAGAGCGTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGc                                                                                                                                                                                                  <  2:214822/239‑1 (MQ=255)
                                        agaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTTGCGGTACTGGAGATGGGGGCATCGATGGTGGGGTATGGCTTGATAAAATTGGcc                                                                                                                                                                                                 >  2:156397/1‑239 (MQ=255)
                                        agaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGcc                                                                                                                                                                                                 <  2:659652/239‑1 (MQ=255)
                                        agaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGc                                                                                                                                                                                                  <  2:605578/238‑1 (MQ=255)
                                            cTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCCCCGTGATGATTTTCAGC‑GGGTTGGCGGTACTGGAGATGGTGGCAGCGATGGGGTGATATCGCTTTGTAAACTTGTGCTggg                                                                                                                                                                                             >  2:114992/1‑238 (MQ=255)
                                            cTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCAACGGAGGTGTGAGATCGCTTGTTAAACTTGGCCCggg                                                                                                                                                                                             >  2:611709/1‑238 (MQ=255)
                                            attGCAAAGGGGAGCCGGGTGGGTTGGGTAGATCAGGAGCTAAAAGGGCTTCGTGATGCGGGGGCTGATGAGGTTCTGGAAAACTTATTTCAGGTTTCCTCGTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTGGCATCGCCTGGGGTAGGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGga                                                                                                                                                                                             <  1:626841/238‑1 (MQ=255)
                                              tGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGGGATGGTGGACTCGATGGGGGGATATCGCTTGGTAAACTTGGCCTGggga                                                                                                                                                                                           >  2:555542/1‑239 (MQ=255)
                                               gCAAAGGGCAGCCCGGACGATCGGTTGGATCAGGCGCTAAAAGTGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGaa                                                                                                                                                                                          <  2:950007/239‑1 (MQ=255)
                                                         gcCCGGGCTTTCGGTTAGGTCAGGTGGTAAAAGTGCTGCGTGTGGTGTGGGCTGATGTGTTTCTGGAAACCTTATTGCAGGTTTCTCCTTCGTGCATGGAAGGCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGtt                                                                                                                                                                               <  1:446124/239‑1 (MQ=255)
                                                         ccccgggcgtttcgTTAGATGAGGCGTTAAAAGGGGTTCGTGATGCGGTGGCTGTTGAGGTTCTGGAAAACTGAGTGCAGGTTTCTGCTTCCCGGTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGGCTGGGGTATGGCGGCCCCCGTGATGATTTGCAGC‑GTGTGGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGtt                                                                                                                                                                               <  1:1100177/226‑1 (MQ=255)
                                                          cccGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GGGTTGGCGGTAACGGGGATGGGGGCATCGGTGGGGTGTATTCGCTTGGTAAACTTGGGCTGGGGAAA‑GTTTAttttt                                                                                                                                                                              >  2:748034/1‑238 (MQ=255)
                                                           ccGGACGAGCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATTCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGtt                                                                                                                                                                               <  2:669781/237‑1 (MQ=255)
                                                                     ggTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGt                                                                                                                                                                   <  2:322365/239‑1 (MQ=255)
                                                                     ggTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGt                                                                                                                                                                   <  2:1113843/239‑1 (MQ=255)
                                                                      gTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGGGATGGTGGCATCGATGGTGTGATATCg                                                                                                                                                                                                                 <  2:441920/193‑1 (MQ=255)
                                                                        ttgATTAGGGGCTAAAAGGGCTTCGTTAGGCGGTGGCTGATGGTGTTCTGGAAAACTTATTGCAGGTTTCTCCTTGGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTgg                                                                                                                                                                <  1:1270284/237‑1 (MQ=255)
                                                                         aGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCAACGTGATGATTTGCAGC‑GTGTTGGAGGGACTGGAGATGGTGGCATCCATGGTGTGATATCGATTGTTTAACTTGGCCTGGGGAAA‑GTTTATTGTCAGGGaaagagatggc                                                                                                                                                               >  2:262356/1‑231 (MQ=255)
                                                                          ggtcgggCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCa                                                                                                                                                              <  2:61357/234‑1 (MQ=255)
                                                                            tCAGGCGCTAAAAGGGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  2:1047488/239‑1 (MQ=255)
                                                                            tCAGGAGCTAAAGGAGCTTCGTGATGCCGTGGCTGCGGAGGTTCGGGAAAACTTATTGCAGGTTTCTCCTTCGCCCTGTGAAGTCATTGTTCTGTATG‑GTTTACATCGCCTGGGGTATGGCGGGCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  2:251036/239‑1 (MQ=255)
                                                                             cAGGGGCTTACAGGGTTTCGTGGTGGGGGGGGGGATGGGGTTCGTGAAACTTGTTTGCTGGTTTTCCCTTCGTGCTTTTAAGCCATTTTTCTGGAGG‑TTTGGCATCGCCTGGGGTATTGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  1:921915/238‑1 (MQ=255)
                                                                                      aaagggtttttttgttggtgggCTGTTGAGGTTTTTGGAAATTTATTGCTGGTTTCTACTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTagg                                                                                                                                                 <  1:591629/219‑1 (MQ=255)
                                                                                      aaaGATATTCGTGATGCGGTGGCTGAGGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCGGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTagg                                                                                                                                                 <  2:763572/239‑1 (MQ=255)
                                                                                       aaGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTaggc                                                                                                                                                <  2:879715/239‑1 (MQ=255)
                                                                                           gCTTCGGTATGGTGTGGTTGATGGGGTTGTGGAATACTTATTGCAGGTTTCTCCGTCGGTCTTTGAACTCATTGTCTTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTaggcaggc                                                                                                                                            <  1:609057/239‑1 (MQ=255)
                                                                                            cTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGTGGGCATCGATGGTGGGGTATCGCTTGATAAACTTGGCCTGGGGAAA‑GTTTATGTTCAGGCAAACAGTTGGCAGAGA‑ACTGTAGGCAGGcc                                                                                                                                           >  2:839014/1‑239 (MQ=255)
                                                                                             ttCGTGATGGGGTGGCTTATGGTGTTCTGGAAAACTTTTTGCAGGTTTCACCTTCGAGTTTGTGAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGcc                                                                                                                                           <  1:395170/238‑1 (MQ=255)
                                                                                                               ttgtggTTTTGGAAAACTTTTTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTAc                                                                                                                       <  1:839014/235‑1 (MQ=255)
                                                                                                                     ttCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              >  2:34950/1‑133 (MQ=255)
                                                                                                                     ttCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  1:34950/133‑1 (MQ=255)
                                                                                                                      tttGGAAAACGTATTGCCGGGTTCTCCTTCGCGCTTTGAAGTCATTGTTCGGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GGGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCa                                                                                                                 <  2:1249943/236‑1 (MQ=255)
                                                                                                                      tgtgagaagtttATTGGCGGTTTCGCTTTCGGGCTTTGACGGTATTGTTCTGGATG‑GTTTGCATGGGCTGGGTTATGGCGGCAGCCGTGATGAGTTGCAGC‑GTGTTGACGGTACTGGAGATGGTGGCATCGATGGTGTGATCTCGCTTGATAACCTTGTCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTGGGCAGAAT‑ACTGGAGGCAGGCAAGAA‑TTACAGGCATTTTACGGCGCa                                                                                                                 <  1:790845/228‑1 (MQ=255)
                                                                                                                            aaaCTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCActggct                                                                                                           <  2:1049650/238‑1 (MQ=255)
                                                                                                                               gttttgtcccttttCTCCTTGGCCCTTTGGATTCATTGTT‑TGGTGGTTTTTGGATCGGCTGGGGAAGTGCGGCGCCCGTGATGATGTGCAGC‑GGGTTGGCGGTAGTGGGGATGGTGGCATCGATGGTGGGATGTCGCTTGATAAACTTGGCCTGGAGAAA‑GGTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTTCGGAGCTCTGGCTGGGc                                                                                                       <  1:1203954/228‑1 (MQ=255)
                                                                                                                                         gTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGGGATGGTGGCATCGATGGGGTGATATCGCTTGATAAACTTGGGCTGGGGAAA‑GGTTATGTTCAGGCAAAACGTTTGCCGAAT‑ACTTTAGGCAGGCCCAAA‑TTTAAGGCCTTTTACGGCGCCCTGGCCGGGCaaaaaaggaa                                                                                             >  2:981761/1‑234 (MQ=255)
                                                                                                                                          tttCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGCTTTGCAGG‑GGGTTGGCGGTACTGGAGAGGGTGGGATCGATGGGGTGATATCGCTTGGTAAACTTGGCCTGGAGAAA‑GGTTATGTTTGGGGAAAACGTTTGGAGAAA‑AATGTAGGGAGGGCCGAATTTAGAGGCTTTTTAGGGCGC‑CTGGGTGGGGAAAAAACGaaa                                                                                            >  1:975280/1‑239 (MQ=255)
                                                                                                                                              gccTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGTCGGCCACCGTGATGATGTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTg                                                                                       <  2:842937/238‑1 (MQ=255)
                                                                                                                                                 ttCGGGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTgggg                                                                                    <  2:975280/239‑1 (MQ=255)
                                                                                                                                                 ttCGCTTTTGGAAGTCTTTGTTTTGGATG‑TTTTGGATCGCTTGGGGTAGGGGGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGGTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTgggg                                                                                    <  1:454329/239‑1 (MQ=255)
                                                                                                                                                     cgcTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGTAGC‑GTGTTGGCGGTACTGGAGAGGGTGGCACCGATGGGGTGATATCGCTTGATAAAATTGGCCTGGAGGAA‑GTTTTTGTTCAGGGAAAACGTTGGCAGAAT‑ACCTTGGGCAGGGCCAAA‑TTACCGGGATTTTATGGCGCACAGGTTTGGGAAAAAAGGAAAC‑AGgggggggt                                                                                >  1:1080649/1‑235 (MQ=255)
                                                                                                                                                                cATTGTTCTGGATG‑GTTTGCATCGACGGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACtt                                                                     <  2:1185873/239‑1 (MQ=255)
                                                                                                                                                                   tttgttggggtggttggCATCGCCTGGGGTAGGGCGGCCTACGGGATGATTTGCAGG‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCt                                                                   <  1:1116275/223‑1 (MQ=255)
                                                                                                                                                                    gTTCTGGATT‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTgg                                                                 <  2:1229833/239‑1 (MQ=255)
                                                                                                                                                                      tCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAAATTGGCCTGGAGAAA‑GGTTATGTTCAGGCAAAAAGTTGGCAGGAT‑ACTGTATGCCGGCCAGAA‑TTACAGGCATTTTACGGGGCCCTGGCTGGGCAAAAAACGCAAC‑GGGGGGGGTTTTTTCTCCCTTCTgggt                                                               >  2:144665/1‑237 (MQ=255)
                                                                                                                                                                        tGGATA‑TTTTGCGTCGCCT‑GGGTATGTCAGCCACCGTGATGATTTGCAGC‑GTGTGGGCGGTGCTGGAGATGGTGGCATCGATGGAGTGATATCGATTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGAAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTa                                                            <  2:728280/239‑1 (MQ=255)
                                                                                                                                                                          ggttgtttttgcttggcTGGGGGGTTGCGGCCACCGGGTGGGTTTGCGGCAGTG‑TGGGGGTTATTGAGATGTTGGGTTCGATGGTGTGATATCGCTTGATAAACTTGGGCTGGAGAAA‑GTGTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTTGAtttt                                                            <  1:755741/223‑2 (MQ=255)
                                                                                                                                                                                tttgcaggggcgggggTATGGGGGCCACCGGGATTATTTGCAGG‑GGGTTGGTGGTACTGGAGGTGGTGGCATCGATGGTGTGAAATCGCTTGATAAACTTGGGCTGGAGAAA‑GTTTATGGTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTtctc                                                      <  1:1297053/227‑1 (MQ=255)
                                                                                                                                                                                tttGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTTCTGCGGATGGTGGCCGCCATTGGGTGATATCGCTTGATAAACTTGGCCTTGGGAAAGGTTTAGGTTCGGGCAAAAC‑TTGGGAGAATAACTG‑AGGCCGGGCCGGA‑TATCAGGGCTTATTCGGCGGACTGGGGGGGAGCAAAGCAAAACACGGGGGTGTTT‑TTTCCCCTTCTGGTTTTatgttcc                                                      >  1:914026/1‑233 (MQ=255)
                                                                                                                                                                                                tATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGTGCACTGGCTGGGCAAAAAGCGAAAC‑GTg                                                                                       <  2:23430/190‑1 (MQ=255)
                                                                                                                                                                                                      gccccccGTGGTGATTGGCGGC‑GGGGTGGGGGTACTGGGGATGGTGGGAGCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTAGGTTCAGGCAAAACGGTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACGTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                 <  1:144665/233‑1 (MQ=255)
                                                                                                                                                                                                          aCCGTGATGATTTGCAGG‑GTGTTGGCGGTACTGGAGATGGTGGCCTCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTc                            <  2:1157002/239‑1 (MQ=255)
                                                                                                                                                                                                            cgggctgatTTGCAGC‑GTGTTGGCGTTACTGGAGATGGTGGGTTCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTAGGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑Ga                         <  2:370865/234‑1 (MQ=255)
                                                                                                                                                                                                            cGTTATGTTTTGCAGC‑GTTTTGGCGGTACTGGAGATGGTTGCATCGGTGGTGTGATATCGCTTGATAAACTTGGGCTGGAGAAA‑GTTTATGTTCAGGCAAAACGGTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑Ga                         <  1:981761/239‑1 (MQ=255)
                                                                                                                                                                                                            cGTGATTCTTTGCAGC‑GTGTTGGCGGTATTGGGGATGGTGTCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCTGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑Ga                         <  1:262356/239‑1 (MQ=255)
                                                                                                                                                                                                                aagatTTGCAGC‑CAGCTGCCGGAACGTGGGATGGTGGCATCGATGGCGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACCGGCATTTTACCGCCCACTGGCTGGGCAAAAAGCGAAAC‑GAGGGGGAGTTCTTACCACTTCCGGAGTTACTTCTCAGGCGCGGGAATGTTCCCAAACCCTC‑Ggggga                     >  1:70861/3‑238 (MQ=255)
                                                                                                                                                                                                                 tgatTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGAATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTa                    <  2:871297/239‑1 (MQ=255)
                                                                                                                                                                                                                  gatTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTTGGATAGCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTa                    <  2:1284461/238‑1 (MQ=255)
                                                                                                                                                                                                                       ttcagtgggtTGGCGTTACTGGAGATGGTGGCATTGATGGTGTGAAATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGgtg               <  2:1143635/231‑1 (MQ=255)
                                                                                                                                                                                                                       tggaggatgtTTGCGGTGCTGGAAATGGTGGCATCGATGGTTTGATATCGCTTGATAAACGTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGgtg               <  2:673356/232‑1 (MQ=255)
                                                                                                                                                                                                                          gC‑GTGTTGGCGGTACTGGAGATTGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            >  1:616496/1‑98 (MQ=255)
                                                                                                                                                                                                                          gC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  2:616496/98‑1 (MQ=255)
                                                                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGGGCACTGGCTGGGCAAAAAGGGAAAC‑GTGGGGTGTTTTTTACCACTCCTGGATTTACTTCTCCGGGGCGGGGCTTTGGCCAATCCGCCTGGGGGTAGGGTgtgggtt         >  1:963941/1‑239 (MQ=255)
                                                                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTggg                                                                                     >  1:226663/1‑163 (MQ=255)
                                                                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTggg                                                                                     <  2:226663/163‑1 (MQ=255)
                                                                                                                                                                                                                               tttGGCGGTTCTGGGGATGGGGGCATGGATGGGTGGTTATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGAAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTgttggttg        <  1:114992/237‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTATTGGAGGTGGGGGCATCGATGGTGTGATACCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTGTTGGTTGATgggg  <  1:156397/238‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTGTTGGTTGATgggg  <  2:70861/238‑1 (MQ=255)
                                                                                                                                                                                                                                       tACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTgttgg           <  2:168593/227‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
CCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATG‑TTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGC‑GTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAA‑GTTTATGTTCAGGCAAAACGTTGGCAGAAT‑ACTGTAGGCAGGCCAGAA‑TTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAAC‑GTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTC‑GAGGGTATGGTGTTGGTTGATGGGG  >  NC_000913/4587308‑4587772

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