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

TGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTT  >  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)
                                                tCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGCGTGTTGGGGGTACTGGGGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTTGCCTGGAGAAAGTTTATGGTCAGGCAAAACGTTGGCGGAATACTGTAGGCAGGCCAg                                                                                                                                                                                                  >  1:182732/1‑239 (MQ=255)
                                                      tggggtgggtgGTGGGGTTTTGGAAAATTTTTTGCAGGTTTCTTCTTCGGGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTAc                                                                                                                                                                                           <  2:771922/236‑1 (MQ=255)
                                                           tGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGGGGCATCGATGGGGTGATATCGCTTGGTAAACTTGGCCTGGGGAAAGGTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCa                                                                                                                                                                                      >  1:82110/1‑239 (MQ=255)
                                                             gCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTTTTTGCGGTACTGGAGATGGTGGACTCGCAGGTGGGGTTTCCCTTGTTTAAATTGTGCTTGGGGAAGATTAGTTTTCGGCAAAAAGTGGGACGAAAAAGTTAGGGGGGCCAGAAATAAAAGCAtt                                                                                                                                                                                    >  2:685232/1‑239 (MQ=255)
                                                                  tggggTTCTGGAAAACTTATTGGAGGTTTCCCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACg                                                                                                                                                                               <  1:319648/236‑1 (MQ=255)
                                                                  tgaGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTAc                                                                                                                                                                                <  1:583342/238‑1 (MQ=255)
                                                                                 cTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTTCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGTGTGATATCGCTTGATAAAATTGGCCTGGGGAAAGTTTATGTTCAGGCAAAAAGTTGGCCGAATACTGTAGGCAGGCCCGAATTACAGGCATTTTACGGGGCCCCGGCTGGGc                                                                                                                                                               >  1:152139/1‑239 (MQ=255)
                                                                                           gTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGGTATCGCCTGATAAACTTTGCCTGGGGAAAGGTTATGTTCCGGCACAACGTTTGCAGAATACTGTAGGCAGGCCCAAATTACAGGCAATTTAAGGCGCCCTGGCTTGGGAAAAAACGaa                                                                                                                                                   >  1:892692/1‑239 (MQ=255)
                                                                                            ttttcccTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGaa                                                                                                                                                   <  1:519468/233‑1 (MQ=255)
                                                                                              tgttctttcgggtttGTAGTCATTGTTCTGGTCGTTGTTCAGCGGCGGGTGTATGGCGGCCAGCGGTATGATTGTCAGCGGGTTGGGGGTACTGGAGATGTGGGCCTCGATGGTGTGATGTGGCTTGATAAACTTGGCCTGGAGAAATTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGaaaag                                                                                                                                                <  2:144648/227‑3 (MQ=255)
                                                                                                cccTTCGGGTTTTGTAGTCTTTTTTTTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCCCCGTGATGATTTGCAGCGTGGTGGCGGTTCTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTg                                                                                                                                              <  2:393618/238‑1 (MQ=255)
                                                                                                     cgcgCTTTGAAGTCCTTGTTCTGGATGTTTTGCATCGGCTGGGGTAAGGAGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCCTCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGgtg                                                                                                                                         <  1:717542/239‑1 (MQ=255)
                                                                                                             gAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTAc                                                                                                                                 <  1:59638/239‑1 (MQ=255)
                                                                                                                           ggATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTTAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACCGCGCACCTGCTGGGGAAAAAGGCGAAAGTGGGGTGGTTAATACCAATTTCGGGTTTa                                                                                                                   >  1:714162/1‑239 (MQ=255)
                                                                                                                            gATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGATGGTGGCATCGATGGGGGGATATCGCTTGGAAAACTTGGCCTGGGGAAAGTTTTTGTTTATGCAAAACGTTGGCAGAAAACTGTATGAAGGCCAGAATTTAAGGGAATTTTATGCGCACCTGTCGGGCAAAAAAGGGAAACAGGGGGGGTTATTTACCCCTCTTGGTTTTa                                                                                                                   >  2:777436/1‑239 (MQ=255)
                                                                                                                                 tttGCATCGCCTGGGGTATGGCGGCCACCGTGGTGATTTGCAGCGTGTTGGCGGTACTGGGGGTGGTGGCATCGCTGGGGGGGTATCTCTTGTGAAAATTTGCCTGGGGAAAGGTTTTGTTCCGGGAAAACTGTGGCCGCATACTGTAGGCCGGACAGCATTACCAGGCTTTTACGGCGGCCGACACGGGCAAAAAGCAGAAAAGGGGGGGGTTTATAACCACTTTTgttttattttta                                                                                                              >  2:470632/1‑228 (MQ=255)
                                                                                                                                             ggggTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGTCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACt                                                                                                <  1:794784/239‑1 (MQ=255)
                                                                                                                                                   tgtcggcCACCGTGATGATTTGCAGCGTGTTGGCGGCCCTGGGGACGGTGGACGCCATGGGGGGATATCGCTTGATTAAATTGGTCTGGATTAATTTTATTGTCAGGCACAAAGTGTGCTGAATAATGTATGCGAGGCCGAATTAAAGGGCATTTAAGGCCGACTGGCTGGGGCAAAAGGTGAACGGGGGGTGTTTTTTTACACATTTGGGTTTACTTTCTCAGGcgcggaattttccc                                                                                           >  2:648243/4‑229 (MQ=255)
                                                                                                                                                    ggcggcCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGGGTGGTATCGCTTGATAAAATTTGCCTGGAGAAAGTTTATGTTTAGGCAAAACGTTGGCAGAAAACTGTAGGCCGGCCAGACTTACAGGCATTTTTCGGCGCACTTGCCGGGGAAAAAGCGGAACATGGGGGGTTTTTTACCCCTTTTGGATTTACTTTTCAGGCGCGGGACCTTGCCCa                                                                                         >  2:783477/1‑239 (MQ=255)
                                                                                                                                                      cggcCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAGGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCaa                                                                                        <  1:648243/238‑1 (MQ=255)
                                                                                                                                                          cACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGt                                                                                   <  1:470632/239‑1 (MQ=255)
                                                                                                                                                           aCCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATAGCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGt                                                                                   <  1:184541/238‑1 (MQ=255)
                                                                                                                                                             cGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGa                                                                                <  1:186069/239‑1 (MQ=255)
                                                                                                                                                                   gttttGCCGCGTGGTGGGGGTTCTGGAGAGGGTGGCATCGGTTGGGTGATAACGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGTACTGGCTGGTCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTa                                                                           <  1:507668/236‑1 (MQ=255)
                                                                                                                                                                      ttGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTTTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCTGGCAAAACGTTGGCAGAATACTGTTGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACATGGGGGGTTTATTACCACTATTGGATTTATTTCTCCGGCGCGTGAATTGTCCCCATcccctcgggtttagga                                                                       >  2:507668/1‑225 (MQ=255)
                                                                                                                                                                       tGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGGTTGATAAACTTGGCCTGGAGAAAGGTTATGTTCAGGCCAAACGTTGGTAGAATACTGTAGGCAGGCCTGAATTACCGGCCTTTTACGGCCAACTGGCTTGGGAACAAGGGGACAGTGGGGTGTTTATTAACCCTTCGGGATTTACTTTTACGGCGCGTGAATTTggcagacaagtcctgggttggggg                                                                      >  2:663656/1‑216 (MQ=255)
                                                                                                                                                                       tGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGgtg                                                                      <  1:783477/239‑1 (MQ=255)
                                                                                                                                                                       tGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGGGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGCCGAGGGTATgggg                                                                      >  1:354448/1‑237 (MQ=255)
                                                                                                                                                                         cAGCGTGTTGGCGGTCCTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGCATACTGTAGGCAGGCCAGATTTACAGGCATTTTACGGCGCCCTGGCTGGGCAAAAAACGAAACGTGGGGGGGTTATTACCACTCTTGGATTTAATTCTCAGGCGCGGGACTTTTCCCACACCGCCGCGGGTATGgttgt                                                                    >  2:836302/1‑236 (MQ=255)
                                                                                                                                                                             gttttGGGGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTgttggt                                                                 <  2:619747/235‑1 (MQ=255)
                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTGAGGCAAGACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGTAAAAAGGGAAAAGTGGGGGGTTTATTCGCCATTTTGGCTTTCCTTCTAAGCGGCGGGACTTTGCCCAATCCGTGGAggggaggggatgtgttg                                                               >  2:304931/1‑225 (MQ=255)
                                                                                                                                                                              tgtTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGGGTTTATTACCACTTCTGGGTTTACTTTTCAGGCGCGTGCCTTTGCCCAAACCGCCGAGGGTGTGGGgtgggttg                                                               >  2:552278/1‑239 (MQ=255)
                                                                                                                                                                               gtTGGCGGTACTGGAGAGTGTGGCATCGATGGTGTGATAGCGCTTGAGTAACTTGGCCGGGAGAAAGGTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGa                                                              <  1:739068/239‑1 (MQ=255)
                                                                                                                                                                                gttgCGGTACGGGTGTTGGGGGCATCGATGGTGTGATTTCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAGAACGTTGGCAGAATTCTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTTGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGAt                                                             <  2:182732/236‑1 (MQ=255)
                                                                                                                                                                                    cgtttcttgggttggtggCATTGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                         <  2:225937/227‑1 (MQ=255)
                                                                                                                                                                                     ggTATTGGGGATGGGGGGATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                         <  2:152139/238‑1 (MQ=255)
                                                                                                                                                                                     ggTACGGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTTGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATgggg                                                         <  1:762077/238‑1 (MQ=255)
                                                                                                                                                                                      gttcTGGAGATGGGGGAATCGATTGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGa                                                        <  2:128861/235‑1 (MQ=255)
                                                                                                                                                                                        aagggAGATGGTCGCATCGATGGTTCGATATCGTTTGATAAACTTGGCCTGGGGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCGCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAc                                                      <  1:735418/235‑1 (MQ=255)
                                                                                                                                                                                         cTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTggg                                                                                                                                            <  1:301937/152‑1 (MQ=255)
                                                                                                                                                                                         cTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTggg                                                                                                                                            >  2:301937/1‑152 (MQ=255)
                                                                                                                                                                                             agaTGGTGGCATTGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCg                                                                                 <  2:56524/206‑1 (MQ=255)
                                                                                                                                                                                             agaTGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCg                                                                                 >  1:56524/1‑206 (MQ=255)
                                                                                                                                                                                                tggtggGATCGATGGTGTGATATCGGTTGAAAAACTTGGCCTGGAGAAAGTTTAAGTTCAGGCAAAACGGTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  2:212083/238‑1 (MQ=255)
                                                                                                                                                                                                tggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  1:552278/238‑1 (MQ=255)
                                                                                                                                                                                                tggtggCATCGATGGTGTGATATCCCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTg                                              <  1:304931/238‑1 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTc                                                                                  <  1:744421/201‑1 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTc                                                                                  >  2:744421/1‑201 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGGGACTTTGCCCAATCCGTCGCGGGAATGGTGTTGGTTGATGGGGAACGCCTGGGGCa                                            >  1:198335/1‑239 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGGGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGc                                             >  1:417280/1‑238 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGAGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGGGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAAGGCCGGGTg                                              >  1:450720/1‑237 (MQ=255)
                                                                                                                                                                                                 ggtggCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGCGAGGGGTTTATTACCCCTTCTGGATTTCCTTCTCAGGCGCGTGACTTTGCCCAATCCGCCGAGGGTATGGGGTTGGTTGATGGGGGACGCCGGGGGc                                             >  2:703988/1‑238 (MQ=255)
                                                                                                                                                                                                      cATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCAAATCCGCCGGGGCGTTGGTGTTGGTTGATGGGGacagccgggggcccctg                                        >  2:735418/1‑222 (MQ=255)
                                                                                                                                                                                                               ttggATATCGCTTGATAAACTTGGCCGGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTaa                                       <  1:843142/227‑1 (MQ=255)
                                                                                                                                                                                                               tgtgATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGCATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGaaa                              >  1:862724/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGGAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCAGTTTACGGCCCCCTGGCTGGGCAAAAAGCGAACCGTGGGGTGTTTACTACCACTTCTGGATTTACTTCTCAGGAGCGTCACTTTGCCCCATCCGTCGGGGGTATGGCGTTTGTGACAGGGGAAACCCTGGTGCCCCTACAGAACCAAAACGAaggagg                    >  2:37168/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCCTGCCAGAATTACAGGCATTTTACGGCTCACTGGCTGGGCAAAAATCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTTCTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGGCGGTCTGTGTTTGGTTGGTGGGGCACGCCTGGTGCACTTAATGATAGAAAACGAAGTAgg                    >  2:776788/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGGGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAgg                    >  1:301741/1‑239 (MQ=255)
                                                                                                                                                                                                                         cTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTAAGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCAATTCTGGATTTACTTCTCAGGCCCGTGACTTTTCCCAATCCGTTGGGGGAATGGGTATGGTGGATCGGGAACGCCTGGCGCCCGTaaacatgcaaaaccttgtcct                    >  2:843142/1‑220 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTgggg                                                                                                                                           <  2:606338/117‑1 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTgggg                                                                                                                                           >  1:606338/1‑117 (MQ=255)
                                                                                                                                                                                                                             aTAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGtt                >  1:687926/1‑239 (MQ=255)
                                                                                                                                                                                                                                    tGGCCTGGGGAAATTTTTTTTTCAGGTAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCAc         <  2:354448/239‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCGGGAGAAAGGTTAAGTTCGGGCAGAACGTTGGCAGAATCCTGTAGGCAGGCCAGGATTACAGGCATTTTGCGGCGCACTGGCTGGGCAAAAAGCGAAACGGGGGGTGTTTATTACCTGTTCTGGATTTACTTCTCAGGTGCGTGACTTTGCGCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGt       <  1:37168/239‑1 (MQ=255)
                                                                                                                                                                                                                                       ccTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGtt      >  1:140634/1‑239 (MQ=255)
                                                                                                                                                                                                                                       ccTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGGTTTGGTTGATGGGGAACGCCTGGTGCCCTTAATGATCGAAAACGAAAAAGGGGTGTTTTCAAGtt      >  2:710044/1‑239 (MQ=255)
                                                                                                                                                                                                                                           gagaAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATTGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGCAAACGAACTAGGGGTTTCTTCCCGTttgtt  >  1:226147/1‑239 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   
TGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGCAAAACGTTGGCAGAATACTGTAGGCAGGCCAGAATTACAGGCATTTTACGGCGCACTGGCTGGGCAAAAAGCGAAACGTGGGGTGTTTATTACCACTTCTGGATTTACTTCTCAGGCGCGTGACTTTGCCCAATCCGTCGAGGGTATGGTGTTGGTTGATGGGGAACGCCTGGTGCACTTAATGATCGAAAACGAAGTAGGGGTTTCTTCACGTTTGTT  >  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: