Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,986,499 A→G D271D (GAT→GAC)  unknown ← 16S rRNA (guanine(1207)‑N(2))‑methyltransferase (EC 2.1.1.172)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,986,4990AG96.3% 198.6 / ‑5.3 81D271D (GAT→GAC) unknown16S rRNA (guanine(1207)‑N(2))‑methyltransferase (EC 2.1.1.172)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (33/45);  minor base C (3/0);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 8.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

AACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCCGTGCACCGAATAAACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCATCATGAAACGGCGGGTTGGAAATGATCAGATCAAATTTGCCGTCAATA‑TGGGAAAACACATCACTGGCAACCACTGTTCCCTGTAGCCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCATGAATATCTGTCATGATTAATCGGGTGTTCGGGTTGTGCCGTTGAATATAACTGCCGATCACTCCCGCGCCGCAGCCCACGTCCAACACCCCGCCTGCGATGTTTTT  >  1/1986261‑1986738
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 
aactaaAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGGCAAAGCCGTCATGaaaaagg                                                                                                                                                                                                                                       >  1:796712/1‑247 (MQ=255)
aactaaAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGTAGCGGTGTCAATGCCGTCATGaaaaagg                                                                                                                                                                                                                                       >  1:310541/1‑247 (MQ=255)
  ctaaAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGGCGTCATGaaaaggtgg                                                                                                                                                                                                                                     >  1:649459/1‑245 (MQ=255)
  ctaaAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGCGTACGCGGGGTCAATACCGTCCTGaaaaagcgg                                                                                                                                                                                                                                     >  1:315432/1‑245 (MQ=255)
          cTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGaaa                                                                                                                                                                                                                             <  1:18630/251‑1 (MQ=255)
             cGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGGCAATGCCGGCATGaaaaagcgggatgggaaaga                                                                                                                                                                                                                          >  1:731724/1‑234 (MQ=255)
             cGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGGCAAAGCCCGCATGAAAAAGCGGGTTgggaaaga                                                                                                                                                                                                                          >  1:578588/1‑245 (MQ=255)
               cTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATgatc                                                                                                                                                                                                                        <  1:268844/251‑1 (MQ=255)
               cTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATgatc                                                                                                                                                                                                                        <  1:309317/251‑1 (MQ=255)
                       gCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGAAGCGGGGGCCATGCCGTCATGAAAAAGCGGAGTGGGAATgatcagatagaa                                                                                                                                                                                                                >  1:597404/1‑247 (MQ=255)
                                  tGCACACTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGCGTAAGCGGGGTCCATGCCGTCATGAAAAGGCGGGTTGGGAAAgaacagatc                                                                                                                                                                                                                   >  1:460007/1‑237 (MQ=255)
                                      cacTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCGGCATGAAAAAGCGGGGTGGAAATgatc                                                                                                                                                                                                                        >  1:719455/1‑228 (MQ=255)
                                      cacTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGGCAATGCCGGCCTGGAAAGGGGGGTTGGGAATGATTAGATCGAATTTggcggcaaaatg                                                                                                                                                                                                 >  1:675427/1‑240 (MQ=255)
                                      cacTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGCACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGAGGGGGCCATGCCCTCCTGAAAAAGCGGATTTGAAATGATTAGATCGAATTTTCCGTCCATA‑Tg                                                                                                                                                                                                >  1:623715/1‑251 (MQ=255)
                                        cTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑Tggg                                                                                                                                                                                              <  1:66558/251‑1 (MQ=255)
                                        cTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑Tggg                                                                                                                                                                                              <  1:829198/251‑1 (MQ=255)
                                        cTTAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCCTCACTAAAAGGCGGATTGGAAATGATCAGATCGAATTTTCCGTCAATA‑Tggg                                                                                                                                                                                              >  1:702141/1‑251 (MQ=255)
                                          tAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATgatc                                                                                                                                                                                                                        <  1:31835/224‑1 (MQ=255)
                                          tAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGaa                                                                                                                                                                                            <  1:729665/251‑1 (MQ=255)
                                          tAACTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGaa                                                                                                                                                                                            <  1:800969/251‑1 (MQ=255)
                                             cTCATCTGCGCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAac                                                                                                                                                                                         <  1:739006/251‑1 (MQ=255)
                                                 tctgcgCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGAAGCGGGGTCAATGGCGTCATGAAAAAGCGGGTTGGGAATGATCAGATCGAATTTGGCGGCAATA‑TGGGaaaaaacgg                                                                                                                                                                                     >  1:318/6‑246 (MQ=255)
                                                      gCACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGAAGGGGTGGCAATGCCGTCATggaaaagcgg                                                                                                                                                                                                                                     >  1:352990/1‑190 (MQ=255)
                                                       cACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTgg                                                                                                                                                                               <  1:740646/251‑1 (MQ=255)
                                                       cACCGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTgg                                                                                                                                                                               <  1:475818/251‑1 (MQ=255)
                                                         ccGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑Tggg                                                                                                                                                                                              <  1:782848/234‑1 (MQ=255)
                                                         ccGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCGTCATGAAAAGGCGGGTTGGGAATGATCAGATCGAAATTGCCGTCAATA‑TGGGGAAAAACGTCACTGGCa                                                                                                                                                                             >  1:327429/1‑251 (MQ=255)
                                                         ccGAATATACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGAAGCGGTGTCCATGCCGTCCTTAAAAAGCGGGTTGGAAATTATCAGATCGAATTTTCCGTCCATA‑TTGGAAAACACACCACTGGCa                                                                                                                                                                             >  1:887063/1‑251 (MQ=255)
                                                              tataCCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCAc                                                                                                                                                                        <  1:235068/248‑1 (MQ=255)
                                                                 aCCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCGTCCTGGAAAGGGGGGTTGGGAATGATCAGATCGAATTTGCCGGCAATATTGGGAAAA‑AAGTCACTGGGAACCAccgg                                                                                                                                                                     >  1:224597/1‑248 (MQ=255)
                                                                   cTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCAccgtgc                                                                                                                                                                   <  1:836938/251‑6 (MQ=255)
                                                                               aTTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGTAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCAccgtgccttgcaaccgat                                                                                                                                                       <  1:448429/251‑18 (MQ=255)
                                                                                     cGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATgatcagat                                                                                                                                                                                                                    <  1:608175/185‑1 (MQ=255)
                                                                                     cGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCAc                                                                                                                                                                        <  1:630124/228‑1 (MQ=255)
                                                                                            aCGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAg                                                                                                                                          <  1:386539/250‑1 (MQ=255)
                                                                                            aCGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAg                                                                                                                                          <  1:506731/250‑1 (MQ=255)
                                                                                             cGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑Tggg                                                                                                                                                                                              <  1:121830/198‑1 (MQ=255)
                                                                                                     cTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGaa                                                                                                                                                                                            <  1:621521/192‑1 (MQ=255)
                                                                                                              tGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑Tggg                                                                                                                                                                                              <  1:334201/181‑1 (MQ=255)
                                                                                                                  gATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTgg                                                                                                                                                                                                                                <  1:314511/144‑1 (MQ=255)
                                                                                                                                    aaGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTc                                                                                                 <  1:626801/251‑1 (MQ=255)
                                                                                                                                          gAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAat                                                                                                        <  1:601032/238‑1 (MQ=255)
                                                                                                                                          gAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCGTCATGAAAAAGCGGATTGGGAATGATCAGATCGAATTTTCCGTCAATA‑TGGGAAAAAACGTCACTGGCAACCACCGGGCCTTGCAACCGATTTTTCTG‑TAAAg                                                                                                                                          >  1:77823/1‑204 (MQ=255)
                                                                                                                                           aaaCGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGa                                                                                             <  1:806713/248‑1 (MQ=255)
                                                                                                                                           aaaCGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGAtt                                                                                           <  1:762075/250‑1 (MQ=255)
                                                                                                                                                cATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑Tggg                                                                                                                                                                                              <  1:803955/147‑1 (MQ=255)
                                                                                                                                                     gCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCGTCATGAAAAGGCGGGTTGGGAATGATCAGATCGAATTTGCCGTCAATA‑TGGGGAAACAAGTCACTGGCAACCACCCTGCCTTGCAACCGATTTTTCTG‑TAAAATGCGCTTTGCGGAAGCCAACGCCAAGGCGGGAATATCCGTCATGATTAATCGGGGGtt                                                                                >  1:597229/1‑251 (MQ=255)
                                                                                                                                                      ccACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCAc                                                                                                                                                                        <  1:654530/163‑1 (MQ=255)
                                                                                                                                                         cGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCggg                                                                            <  1:25843/251‑1 (MQ=255)
                                                                                                                                                           aTGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGt                                                                           <  1:2235/250‑1 (MQ=255)
                                                                                                                                                           aTGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGGAAGCGGGGTCAAAGGCGTCATGAAAAGGGGGGATGGGAAAGATTAAAACGAATTTGGCGTTAATA‑TGGGGAAACACGGTACCGGGAACCACCCGGCCCTGCAGCCGATTTTTCTGTTAAAG‑GCGCCTTGCGGGAGCCAACGCCAAGGCGGGGAAATTCGTCATGATTAATCggggggttggggt                                                                          >  1:456001/1‑241 (MQ=255)
                                                                                                                                                             gcgcAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAAGCGGGTTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAAAACGTCACTGGGAACCACCGTGCCTTGCAACCGCTTTTCCTG‑TAAAGTGCGCTTTGCGGGAGCCAACGCCAAGGCGTGGATATCCGTCCTGATTAATctgggggtcgggg                                                                           >  1:652145/1‑236 (MQ=255)
                                                                                                                                                             gcgcAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCGTCATGAAAAGGCGGATTGGGAATGATCAGATTGAATTTTCCGGCAATA‑TGGGGAAAAACGGCAATGGCAACCACCCTGCCTTGCAACCGATTTTTCTG‑TAAAGTGCGCTTTGCGGGAGCCAACGCCCAGGGGGGAATATTCCTCATGATTAATCGGGGGTTCggggtgg                                                                        >  1:56099/1‑247 (MQ=255)
                                                                                                                                                                cAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGc                                                                      <  1:116878/250‑1 (MQ=255)
                                                                                                                                                                cAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGc                                                                      <  1:725967/250‑1 (MQ=255)
                                                                                                                                                                  gCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAACGGGGGCAATGCCGTCCTTAAAAAGCGGGTTGGGAAAGATCAGATCGAATTTGCCGGCCATA‑TGGGGAAAAACGTCACTGGCCACCCCCGGGCCTTTCAACCGATTTTTCTG‑TAAAGAGCGCTTTGCGGGAGCCCACGCCAAGGCGTGAATATCCGTCCTGATTAATCGGGGGTTCGGGGTGTGctgg                                                                   >  1:889829/1‑248 (MQ=255)
                                                                                                                                                                          cccGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAatat                                                            <  1:294368/250‑1 (MQ=255)
                                                                                                                                                                            cGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAAc                                                         <  1:763838/251‑1 (MQ=255)
                                                                                                                                                                                ttCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCGTCCTGAAAAGGCGGGTTGGGAATGATCAGATCGAATTTTCCGTCAATA‑TGGGGAAAAACGTCACTGGCAACCACCCTGCCTTGCAAACGATTTTCCTG‑TAAAATGCGCTTTTCGGAAGACAACGCCAAGGCGGGAATATCCGTCATGATTAAACGGGGGTTCGGGGTGTGCTGTTGAATATAAATGc                                                      >  1:830181/1‑250 (MQ=255)
                                                                                                                                                                                  cAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAAAGCCGTCATGAAAAGGGGGATTGGGAATGAAAAGATCGAATTTGCCGTCAATA‑TGGGAAAACAAGACACTGGCAACCACCGTGCCCTGCAAACGATTTTCCTG‑TAAAGTGCGCTTTTTGGAAGCCAACGCCAAGGCGGGGATATCCTTCATGGTTAAATGGGTGGTCGGGTTGTGCTGTTGTATATAAATGccaa                                                   >  1:220091/1‑249 (MQ=255)
                                                                                                                                                                                     aTGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCggg                                                                                    <  1:673398/215‑1 (MQ=255)
                                                                                                                                                                                        ccATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCAc                                                                                                                                                                        <  1:631068/129‑1 (MQ=255)
                                                                                                                                                                                         cATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCAcccc                                            <  1:549221/251‑4 (MQ=255)
                                                                                                                                                                                         cATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCAcccc                                            <  1:850484/251‑4 (MQ=255)
                                                                                                                                                                                         cATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGTCGTCCTGAAAAGGCGGGTTGGGAAAGATCAGATCGAATTTTTCGTCAATA‑TGGGAAAAAACGTCACTGGCAACCACCGCGCCTTGCAAACGATTTTTCTG‑TTAAGTGCGCTTTTCGGGAGCCAAAGCCAAGGCGGGAATATCCGTCATGATTAATCGGGGGTTCGGGGTGTGCTGTTGAATATAAATGCCAATCAccc                                             >  1:479260/1‑248 (MQ=255)
                                                                                                                                                                                                   gTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCGTCATGAAAAGGCGGGTTGGGAAAGATCAGATCGAATTTGCCGTCAATA‑TGGGGAAACACGTCCCTGGCAAACACCCTGCCTTGCAACCGATTTTTCTG‑TAAAGGGCGCTTTGCGGAAGCCAACGCCAAGGCGGGAATATCCGTCCTGATTAATCGGGTGTTCGGGGTGTGCTGTTGGatat                                                            >  1:316686/1‑225 (MQ=255)
                                                                                                                                                                                                   gTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAAAGCCGTCATGAAAAAGCGGGTTGGAAAAGAACAAATCGAATTTGCCGTCAAAA‑AGGGGAAAAACGTCACTGGCAAACACCGTGCCCTGCAACCGATTTTTCTG‑TAAAGAGCGCTTTGCGGAAGCCAACGCCCAGGCGTGAATATCCGGCATGATTAATCGGGGGTTCGGGGTGGGCTGTTGAATATAAATGCCAATCAACCCTGCGCCGCAg                                  >  1:482006/1‑251 (MQ=255)
                                                                                                                                                                                                      ttATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGccc                               <  1:358242/251‑1 (MQ=255)
                                                                                                                                                                                                       tATCAGTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCGTCCTGAAAAAGCGGGTTGGGAATGATCAGATCGAATTTTCCGTCAATA‑TTGGGAAACACGTCACTGGCAACCACCCTGCCCTGCAACCGATTTTTCTT‑TAAAATGCGCTTTTCGGGAGCCAACGCCCAGGCGGGGATATCCCTCCTGATTAATCTGGGGTTCGGGGTTTGCTGTTGAATATAAATGCCCATCACCCCCGcgcccccgccc                               >  1:140971/1‑243 (MQ=255)
                                                                                                                                                                                                            gTTCGTTAACGGCGGTGTAAGCGGGGTCAATGCCGTCATGAAAAAGCGGATTGGGAATGATCAGATCGAATTTGCCGTCCATA‑TGGGGAAACACGTCACTGGCAACCACCCTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGGAGCCAACGCCCAGGCGTGAATATCCGTCATGATTAATCGGGGGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCCGCGCCGCAGCCCCCGTc                          >  1:849997/1‑250 (MQ=255)
                                                                                                                                                                                                                 ttAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAacac                    <  1:362688/251‑1 (MQ=255)
                                                                                                                                                                                                                 ttAACGGCGGTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAacac                    <  1:86979/251‑1 (MQ=255)
                                                                                                                                                                                                                     cggcggTGTAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTg                                                       >  1:323483/1‑212 (MQ=255)
                                                                                                                                                                                                                          gtgtAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTgcgc                                                                                                                                     >  1:461547/1‑129 (MQ=255)
                                                                                                                                                                                                                          gtgtAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGc                                                                                                                    <  1:556513/146‑1 (MQ=255)
                                                                                                                                                                                                                          gtgtAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCggg                                                                                    >  1:379572/1‑178 (MQ=255)
                                                                                                                                                                                                                           tgtAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCAcccc                                            <  1:440843/217‑4 (MQ=255)
                                                                                                                                                                                                                            gtAAGCGGTGTCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAacac                    >  1:589367/1‑240 (MQ=255)
                                                                                                                                                                                                                                   gtgtCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGATTAATCGGGTGTTCGGGTTGTGCTGTTGAATATAACTGCCAATCACCCCTGCGCCGCAGCCCACGTCCAACACCTCGCCCGCGATGttttt  >  1:274544/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gtCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATg                                                                                              <  1:525035/157‑1 (MQ=255)
                                                                                                                                                                                                                                     gtCAATGCCGTCATGAAAAGGCGGATTGGAAATGATCAGATCGAATTTGCCGTCAATA‑TGGGAAAACACGTCACTGGCAACCACCGTGCCTTGCAACCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCGTGAATATCCGTCATGAtt                                                                                           >  1:877830/1‑160 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                 
AACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCCGTGCACCGAATAAACCTTAAACTTGTTATTTCTCGCCAGGACGTCATGACTGCCGAAATGTTGATTCAGCCAGTCGGCGTAAGGCAGAAACGCATTTGCCACGATGCGCAGCTCACCGCCCGGTTTCAAATGCCATTTTGCCTGTTTTATCAGTTCGTTAACGGCGGTGTAAGCGGTGTCAATGCCATCATGAAACGGCGGGTTGGAAATGATCAGATCAAATTTGCCGTCAATA‑TGGGAAAACACATCACTGGCAACCACTGTTCCCTGTAGCCGATTTTCCTG‑TAAAGTGCGCTTTGCGGAAGCCAACGCCAAGGCATGAATATCTGTCATGATTAATCGGGTGTTCGGGTTGTGCCGTTGAATATAACTGCCGATCACTCCCGCGCCGCAGCCCACGTCCAACACCCCGCCTGCGATGTTTTT  >  1/1986261‑1986738

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