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→GGGmrrmethylated 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

CCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGG  >  NC_000913/4587308‑4587772
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
ccctCGTTTTAGATCCGAAACCCGTCTTCTCCGGTTATGAAGAACTTGCAAAGAGCAGCCCGGACGGTCGGTTAGAACAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGAt                                                                                                                                                                                                                                         <  1:547755/236‑1 (MQ=255)
    cGATTTAGAGCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCACGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGAtggt                                                                                                                                                                                                                                      <  2:966230/238‑1 (MQ=255)
           gATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  2:560915/239‑1 (MQ=255)
             tCCCAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGCCGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  2:6516/237‑1 (MQ=255)
             gccGAAAGCCGACTCGCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  2:887284/236‑1 (MQ=255)
                 aaaGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGACGTCATTGTTCTGGCTGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGTTGATTTGCATCGTGTTGGCGGTACTGGGGAGGGGGGCACCGgtggtgg                                                                                                                                                                                                                        >  2:1297053/1‑238 (MQ=255)
                 aaaGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCCCCTTGTTGATTTGCAGCGGTTTTGCGGTACTGTAGATGGTGGGACCGCTggggg                                                                                                                                                                                                                        >  2:550188/1‑236 (MQ=255)
                       gaatctcCCGAACATGAAGAACTTGCAAAGGGCCGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATc                                                                                                                                                                                                                  <  2:781587/236‑1 (MQ=255)
                         ttcgtccgactctGAAGAACGTGAAAAGGGCCGGCCGGGAGATCGGTTAGATCAGGCCGTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCTGGATTCTGCTTCGCGCTTTGAAGTCATTGTTGTGGATGTTTTTCTTCGCCTGGGGTATGGCGGCCAACGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGc                                                                                                                                                                                                                <  1:1100579/227‑1 (MQ=255)
                         ctctCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCCCCGGGATGATTTTCAGCGGGTTGGCGGTACTGGGGATGGTGGCATCGCTGGGGTGaattc                                                                                                                                                                                                                  >  2:1116275/1‑233 (MQ=255)
                         accccccggcgaTGGAGGACTTGCAAGAGGTAGCCCGGAGCATCTGTGAGATCAGGGTGTAAATGGTCTCCGTGCGGGTGGGGCTGATGATGTTCTGGAAAACTTATGTCAGGTTTCTCCTTTGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCGCCGTGATGATTTGGAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGc                                                                                                                                                                                                                <  1:729183/228‑1 (MQ=255)
                             cccGATCATGAAGAACTTGCAAAGGGCCGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTTATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGa                                                                                                                                                                                                            <  1:283800/239‑1 (MQ=255)
                                  tCATTAAGAACTTGCAAAGAGCAGCCCGGACGATCGGGTTGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACGTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAAc                                                                                                                                                                                                       <  1:845614/239‑1 (MQ=255)
                                      gaagaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCCTCGATGGTGTGATATCGCTTGATAAACttt                                                                                                                                                                                                    >  2:454329/1‑237 (MQ=255)
                                       aagtaCTTGCAAAGGGCCGTCCGGACGATGGGTTAGATCGGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGc                                                                                                                                                                                                  <  1:732877/239‑1 (MQ=255)
                                       aaggaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGGTAAAAGAGCGTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGc                                                                                                                                                                                                  <  2:214822/239‑1 (MQ=255)
                                        agaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTTGCGGTACTGGAGATGGGGGCATCGATGGTGGGGTATGGCTTGATAAAATTGGcc                                                                                                                                                                                                 >  2:156397/1‑239 (MQ=255)
                                        agaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGcc                                                                                                                                                                                                 <  2:659652/239‑1 (MQ=255)
                                        agaaCTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGc                                                                                                                                                                                                  <  2:605578/238‑1 (MQ=255)
                                            cTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCCCCGTGATGATTTTCAGCGGGTTGGCGGTACTGGAGATGGTGGCAGCGATGGGGTGATATCGCTTTGTAAACTTGTGCTggg                                                                                                                                                                                             >  2:114992/1‑238 (MQ=255)
                                            cTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCAACGGAGGTGTGAGATCGCTTGTTAAACTTGGCCCggg                                                                                                                                                                                             >  2:611709/1‑238 (MQ=255)
                                            attGCAAAGGGGAGCCGGGTGGGTTGGGTAGATCAGGAGCTAAAAGGGCTTCGTGATGCGGGGGCTGATGAGGTTCTGGAAAACTTATTTCAGGTTTCCTCGTCGCGCTTTGAAGTCATTGTTCTGGATGTTTGGCATCGCCTGGGGTAGGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGga                                                                                                                                                                                             <  1:626841/238‑1 (MQ=255)
                                              tGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGACTCGATGGGGGGATATCGCTTGGTAAACTTGGCCTGggga                                                                                                                                                                                           >  2:555542/1‑239 (MQ=255)
                                               gCAAAGGGCAGCCCGGACGATCGGTTGGATCAGGCGCTAAAAGTGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGaa                                                                                                                                                                                          <  2:950007/239‑1 (MQ=255)
                                                         gcCCGGGCTTTCGGTTAGGTCAGGTGGTAAAAGTGCTGCGTGTGGTGTGGGCTGATGTGTTTCTGGAAACCTTATTGCAGGTTTCTCCTTCGTGCATGGAAGGCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGtt                                                                                                                                                                               <  1:446124/239‑1 (MQ=255)
                                                         ccccgggcgtttcgTTAGATGAGGCGTTAAAAGGGGTTCGTGATGCGGTGGCTGTTGAGGTTCTGGAAAACTGAGTGCAGGTTTCTGCTTCCCGGTTTGAAGTCATTGTTCTGGATGTTTTGCATCGGCTGGGGTATGGCGGCCCCCGTGATGATTTGCAGCGTGTGGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGtt                                                                                                                                                                               <  1:1100177/226‑1 (MQ=255)
                                                          cccGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGGGTTGGCGGTAACGGGGATGGGGGCATCGGTGGGGTGTATTCGCTTGGTAAACTTGGGCTGGGGAAAGTTTAttttt                                                                                                                                                                              >  2:748034/1‑238 (MQ=255)
                                                           ccGGACGAGCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATTCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGtt                                                                                                                                                                               <  2:669781/237‑1 (MQ=255)
                                                                     ggTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGt                                                                                                                                                                   <  2:322365/239‑1 (MQ=255)
                                                                     ggTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGt                                                                                                                                                                   <  2:1113843/239‑1 (MQ=255)
                                                                      gTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGTGTGATATCg                                                                                                                                                                                                                 <  2:441920/193‑1 (MQ=255)
                                                                        ttgATTAGGGGCTAAAAGGGCTTCGTTAGGCGGTGGCTGATGGTGTTCTGGAAAACTTATTGCAGGTTTCTCCTTGGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTgg                                                                                                                                                                <  1:1270284/237‑1 (MQ=255)
                                                                         aGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCAACGTGATGATTTGCAGCGTGTTGGAGGGACTGGAGATGGTGGCATCCATGGTGTGATATCGATTGTTTAACTTGGCCTGGGGAAAGTTTATTGTCAGGGaaagagatggc                                                                                                                                                               >  2:262356/1‑231 (MQ=255)
                                                                          ggtcgggCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCa                                                                                                                                                              <  2:61357/234‑1 (MQ=255)
                                                                            tCAGGCGCTAAAAGGGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  2:1047488/239‑1 (MQ=255)
                                                                            tCAGGAGCTAAAGGAGCTTCGTGATGCCGTGGCTGCGGAGGTTCGGGAAAACTTATTGCAGGTTTCTCCTTCGCCCTGTGAAGTCATTGTTCTGTATGGTTTACATCGCCTGGGGTATGGCGGGCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  2:251036/239‑1 (MQ=255)
                                                                             cAGGGGCTTACAGGGTTTCGTGGTGGGGGGGGGGATGGGGTTCGTGAAACTTGTTTGCTGGTTTTCCCTTCGTGCTTTTAAGCCATTTTTCTGGAGGTTTGGCATCGCCTGGGGTATTGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  1:921915/238‑1 (MQ=255)
                                                                                      aaagggtttttttgttggtgggCTGTTGAGGTTTTTGGAAATTTATTGCTGGTTTCTACTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTagg                                                                                                                                                 <  1:591629/219‑1 (MQ=255)
                                                                                      aaaGATATTCGTGATGCGGTGGCTGAGGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCGGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTagg                                                                                                                                                 <  2:763572/239‑1 (MQ=255)
                                                                                       aaGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTaggc                                                                                                                                                <  2:879715/239‑1 (MQ=255)
                                                                                           gCTTCGGTATGGTGTGGTTGATGGGGTTGTGGAATACTTATTGCAGGTTTCTCCGTCGGTCTTTGAACTCATTGTCTTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTaggcaggc                                                                                                                                            <  1:609057/239‑1 (MQ=255)
                                                                                            cTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGTGGGCATCGATGGTGGGGTATCGCTTGATAAACTTGGCCTGGGGAAAGTTTATGTTCAGGCAAACAGTTGGCAGAGAACTGTAGGCAGGcc                                                                                                                                           >  2:839014/1‑239 (MQ=255)
                                                                                             ttCGTGATGGGGTGGCTTATGGTGTTCTGGAAAACTTTTTGCAGGTTTCACCTTCGAGTTTGTGAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGcc                                                                                                                                           <  1:395170/238‑1 (MQ=255)
                                                                                                               ttgtggTTTTGGAAAACTTTTTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTAc                                                                                                                       <  1:839014/235‑1 (MQ=255)
                                                                                                                     ttCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              >  2:34950/1‑133 (MQ=255)
                                                                                                                     ttCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGa                                                                                                                                                                                                                              <  1:34950/133‑1 (MQ=255)
                                                                                                                      tttGGAAAACGTATTGCCGGGTTCTCCTTCGCGCTTTGAAGTCATTGTTCGGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGGGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCa                                                                                                                 <  2:1249943/236‑1 (MQ=255)
                                                                                                                      tgtgagaagtttATTGGCGGTTTCGCTTTCGGGCTTTGACGGTATTGTTCTGGATGGTTTGCATGGGCTGGGTTATGGCGGCAGCCGTGATGAGTTGCAGCGTGTTGACGGTACTGGAGATGGTGGCATCGATGGTGTGATCTCGCTTGATAACCTTGTCCTGGAGAAAGTTTATGTTCAGGCAAAACGTGGGCAGAATACTGGAGGCAGGCAAGAATTACAGGCATTTTACGGCGCa                                                                                                                 <  1:790845/228‑1 (MQ=255)
                                                                                                                            aaaCTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCActggct                                                                                                           <  2:1049650/238‑1 (MQ=255)
                                                                                                                               gttttgtcccttttCTCCTTGGCCCTTTGGATTCATTGTTTGGTGGTTTTTGGATCGGCTGGGGAAGTGCGGCGCCCGTGATGATGTGCAGCGGGTTGGCGGTAGTGGGGATGGTGGCATCGATGGTGGGATGTCGCTTGATAAACTTGGCCTGGAGAAAGGTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTTCGGAGCTCTGGCTGGGc                                                                                                       <  1:1203954/228‑1 (MQ=255)
                                                                                                                                         gTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGGGTGATATCGCTTGATAAACTTGGGCTGGGGAAAGGTTATGTTCAGGCAAAACGTTTGCCGAATACTTTAGGCAGGCCCAAATTTAAGGCCTTTTACGGCGCCCTGGCCGGGCaaaaaaggaa                                                                                             >  2:981761/1‑234 (MQ=255)
                                                                                                                                          tttCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGCTTTGCAGGGGGTTGGCGGTACTGGAGAGGGTGGGATCGATGGGGTGATATCGCTTGGTAAACTTGGCCTGGAGAAAGGTTATGTTTGGGGAAAACGTTTGGAGAAAAATGTAGGGAGGGCCGAATTTAGAGGCTTTTTAGGGCGCCTGGGTGGGGAAAAAACGaaa                                                                                            >  1:975280/1‑239 (MQ=255)
                                                                                                                                              gccTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGTCGGCCACCGTGATGATGTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTg                                                                                       <  2:842937/238‑1 (MQ=255)
                                                                                                                                                 ttCGGGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTgggg                                                                                    <  2:975280/239‑1 (MQ=255)
                                                                                                                                                 ttCGCTTTTGGAAGTCTTTGTTTTGGATGTTTTGGATCGCTTGGGGTAGGGGGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGGTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTgggg                                                                                    <  1:454329/239‑1 (MQ=255)
                                                                                                                                                     cgcTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGTAGCGTGTTGGCGGTACTGGAGAGGGTGGCACCGATGGGGTGATATCGCTTGATAAAATTGGCCTGGAGGAAGTTTTTGTTCAGGGAAAACGTTGGCAGAATACCTTGGGCAGGGCCAAATTACCGGGATTTTATGGCGCACAGGTTTGGGAAAAAAGGAAACAGgggggggt                                                                                >  1:1080649/1‑235 (MQ=255)
                                                                                                                                                                cATTGTTCTGGATGGTTTGCATCGACGGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACtt                                                                     <  2:1185873/239‑1 (MQ=255)
                                                                                                                                                                   tttgttggggtggttggCATCGCCTGGGGTAGGGCGGCCTACGGGATGATTTGCAGGGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCt                                                                   <  1:1116275/223‑1 (MQ=255)
                                                                                                                                                                    gTTCTGGATTTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTgg                                                                 <  2:1229833/239‑1 (MQ=255)
                                                                                                                                                                      tCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAAATTGGCCTGGAGAAAGGTTATGTTCAGGCAAAAAGTTGGCAGGATACTGTATGCCGGCCAGAATTACAGGCATTTTACGGGGCCCTGGCTGGGCAAAAAACGCAACGGGGGGGGTTTTTTCTCCCTTCTgggt                                                               >  2:144665/1‑237 (MQ=255)
                                                                                                                                                                        tGGATATTTTGCGTCGCCTGGGTATGTCAGCCACCGTGATGATTTGCAGCGTGTGGGCGGTGCTGGAGATGGTGGCATCGATGGAGTGATATCGATTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGAAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTa                                                            <  2:728280/239‑1 (MQ=255)
                                                                                                                                                                          ggttgtttttgcttggcTGGGGGGTTGCGGCCACCGGGTGGGTTTGCGGCAGTGTGGGGGTTATTGAGATGTTGGGTTCGATGGTGTGATATCGCTTGATAAACTTGGGCTGGAGAAAGTGTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTTGAtttt                                                            <  1:755741/223‑2 (MQ=255)
                                                                                                                                                                                tttgcaggggcgggggTATGGGGGCCACCGGGATTATTTGCAGGGGGTTGGTGGTACTGGAGGTGGTGGCATCGATGGTGTGAAATCGCTTGATAAACTTGGGCTGGAGAAAGTTTATGGTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTtctc                                                      <  1:1297053/227‑1 (MQ=255)
                                                                                                                                                                                tttGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTTCTGCGGATGGTGGCCGCCATTGGGTGATATCGCTTGATAAACTTGGCCTTGGGAAAGGTTTAGGTTCGGGCAAAACTTGGGAGAATAACTGAGGCCGGGCCGGATATCAGGGCTTATTCGGCGGACTGGGGGGGAGCAAAGCAAAACACGGGGGTGTTTTTTCCCCTTCTGGTTTTatgttcc                                                      >  1:914026/1‑233 (MQ=255)
                                                                                                                                                                                                tATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGTGCACTGGCTGGGCAAAAAGCGAAACGTg                                                                                       <  2:23430/190‑1 (MQ=255)
                                                                                                                                                                                                      gccccccGTGGTGATTGGCGGCGGGGTGGGGGTACTGGGGATGGTGGGAGCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTAGGTTCAGGCAAAACGGTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACGTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAAt                                 <  1:144665/233‑1 (MQ=255)
                                                                                                                                                                                                          aCCGTGATGATTTGCAGGGTGTTGGCGGTACTGGAGATGGTGGCCTCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTc                            <  2:1157002/239‑1 (MQ=255)
                                                                                                                                                                                                            cgggctgatTTGCAGCGTGTTGGCGTTACTGGAGATGGTGGGTTCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTAGGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGa                         <  2:370865/234‑1 (MQ=255)
                                                                                                                                                                                                            cGTTATGTTTTGCAGCGTTTTGGCGGTACTGGAGATGGTTGCATCGGTGGTGTGATATCGCTTGATAAACTTGGGCTGGAGAAAGTTTATGTTCAGGCAAAACGGTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGa                         <  1:981761/239‑1 (MQ=255)
                                                                                                                                                                                                            cGTGATTCTTTGCAGCGTGTTGGCGGTATTGGGGATGGTGTCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCTGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGa                         <  1:262356/239‑1 (MQ=255)
                                                                                                                                                                                                                aagatTTGCAGCCAGCTGCCGGAACGTGGGATGGTGGCATCGATGGCGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACCGGCATTTTACCGCCCACTGGCTGGGCAAAAAGCGAAACGAGGGGGAGTTCTTACCACTTCCGGAGTTACTTCTCAGGCGCGGGAATGTTCCCAAACCCTCGgggga                     >  1:70861/3‑238 (MQ=255)
                                                                                                                                                                                                                 tgatTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGAATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTa                    <  2:871297/239‑1 (MQ=255)
                                                                                                                                                                                                                  gatTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTTGGATAGCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTa                    <  2:1284461/238‑1 (MQ=255)
                                                                                                                                                                                                                       ttcagtgggtTGGCGTTACTGGAGATGGTGGCATTGATGGTGTGAAATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGgtg               <  2:1143635/231‑1 (MQ=255)
                                                                                                                                                                                                                       tggaggatgtTTGCGGTGCTGGAAATGGTGGCATCGATGGTTTGATATCGCTTGATAAACGTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGgtg               <  2:673356/232‑1 (MQ=255)
                                                                                                                                                                                                                          gCGTGTTGGCGGTACTGGAGATTGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            >  1:616496/1‑98 (MQ=255)
                                                                                                                                                                                                                          gCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGa                                                                                                                                                            <  2:616496/98‑1 (MQ=255)
                                                                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGGGCACTGGCTGGGCAAAAAGGGAAACGTGGGGTGTTTTTTACCACTCCTGGATTTACTTCTCCGGGGCGGGGCTTTGGCCAATCCGCCTGGGGGTAGGGTgtgggtt         >  1:963941/1‑239 (MQ=255)
                                                                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTggg                                                                                     >  1:226663/1‑163 (MQ=255)
                                                                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTggg                                                                                     <  2:226663/163‑1 (MQ=255)
                                                                                                                                                                                                                               tttGGCGGTTCTGGGGATGGGGGCATGGATGGGTGGTTATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGAAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTgttggttg        <  1:114992/237‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTATTGGAGGTGGGGGCATCGATGGTGTGATACCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg  <  1:156397/238‑1 (MQ=255)
                                                                                                                                                                                                                                     ggTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg  <  2:70861/238‑1 (MQ=255)
                                                                                                                                                                                                                                       tACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTgttgg           <  2:168593/227‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          
CCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGG  >  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: