Predicted mutation
evidence seq id position mutation annotation gene description
RA AM406671 1,588,245 2 bp→AT coding (837‑838/1170 nt) ugd ← UDP‑glucose 6‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM4066711,588,2450TA99.4% 3178.6 / NA 845T280S (ACA→TCA) ugdUDP‑glucose 6‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  major base A (418/416);  minor base C (2/3);  total (423/422)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*AM4066711,588,2460AT99.8% 3260.5 / NA 845N279K (AAT→AAAugdUDP‑glucose 6‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (420/421);  total (423/422)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 847 total aligned reads displayed.

TTGCAAGAGTTTCATCACATCAATAATCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTTGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGTATTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATTATGTGTCCCAATCCGTGAATCCAACCCAACAC  >  AM406671/1588102‑1588388
                                                                                                                                               ||                                                                                                                                              
   cAAGAGTTTCATCACATCAATAATCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACaa                                                                                                                                         >  1:6824647/1‑149 (MQ=255)
    aaGAGTTTCATCACATCAATAATCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAAt                                                                                                                                        >  1:6518146/1‑149 (MQ=255)
     agagTTTCATCACATCAATAATCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATa                                                                                                                                       >  2:10205295/1‑149 (MQ=255)
     agagTTTCATCACATCAATAATCGAACTCTATCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGACATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATa                                                                                                                                       <  1:8661048/149‑1 (MQ=255)
      gagTTTCATCACATCAATAATCGAACTCTGTCTGAAATTATCTGACTCAGCTTTCATGACCAACCGGTAAATTCCACCAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAg                                                                                                                                      <  2:3933529/149‑1 (MQ=255)
       agTTTCATCACATCAATAATCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGc                                                                                                                                     >  2:1326917/1‑149 (MQ=255)
        gTTTCATCACATCAATAATCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGcc                                                                                                                                    <  1:9267471/149‑1 (MQ=255)
         tttCATCACATCAATAATCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCg                                                                                                                                   >  2:2541373/1‑149 (MQ=255)
            catcacatcaATAATCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGtca                                                                                                                                <  1:7609582/149‑1 (MQ=255)
              tcacatcaATAATCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGtcatc                                                                                                                              >  1:3399199/1‑149 (MQ=255)
                    caATAATCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTcc                                                                                                                        <  1:8172154/149‑1 (MQ=255)
                       taatCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGa                                                                                                                     >  2:2682643/1‑149 (MQ=255)
                       taatCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGa                                                                                                                     <  1:8299494/149‑1 (MQ=255)
                        aatCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAg                                                                                                                    <  1:3560069/149‑1 (MQ=255)
                         atCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAgg                                                                                                                   <  1:4420366/149‑1 (MQ=255)
                            gAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAgg                                                                                                                   >  1:5269628/1‑146 (MQ=255)
                             aaCTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACg                                                                                                               >  2:4520303/1‑149 (MQ=255)
                              aCTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGc                                                                                                              >  1:9443029/1‑149 (MQ=255)
                                tctgtctgAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGccc                                                                                                            >  1:7672186/1‑149 (MQ=255)
                                tctgtctgAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGccc                                                                                                            <  2:482169/149‑1 (MQ=255)
                                 ctgtctgAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCt                                                                                                           <  1:1029111/149‑1 (MQ=255)
                                   gtctgAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGTATTTTTATCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCa                                                                                                         >  2:8046668/1‑149 (MQ=255)
                                    tctgAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGTTGTTTTATCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCAt                                                                                                        >  2:5937088/1‑149 (MQ=255)
                                     ctgAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATa                                                                                                       <  1:9832584/149‑1 (MQ=255)
                                        aaaTTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACt                                                                                                    <  2:6377232/149‑1 (MQ=255)
                                        aaaTTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACt                                                                                                    <  2:3376417/149‑1 (MQ=255)
                                          aTTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTgg                                                                                                  <  1:1533091/149‑1 (MQ=255)
                                          aTTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTgg                                                                                                  >  1:2794499/1‑149 (MQ=255)
                                            tATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATACTACGCTCAACAATGTGCTGTTCTCTGGTTTCGTTTGATTTGACCATAGCCGTCATCCACTCCTGAGGAACGCCCTCATAAGTGGcc                                                                                                >  2:3490663/1‑149 (MQ=255)
                                                  aCCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAgc                                                                                          >  1:2313471/1‑149 (MQ=255)
                                                  aCCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAgc                                                                                          >  1:3245399/1‑149 (MQ=255)
                                                       gCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCtt                                                                                     <  1:376248/149‑1 (MQ=255)
                                                        cTTTCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCttt                                                                                    <  2:5401097/149‑1 (MQ=255)
                                                         tttCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTg                                                                                   >  1:782152/1‑149 (MQ=255)
                                                         tttCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTg                                                                                   >  1:1293831/1‑149 (MQ=255)
                                                           tCATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAg                                                                                 <  2:7761454/149‑1 (MQ=255)
                                                            cATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGt                                                                                >  2:51814/1‑149 (MQ=255)
                                                            cATGACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGt                                                                                >  2:8182422/1‑149 (MQ=255)
                                                               gACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCtt                                                                             <  1:4877950/149‑1 (MQ=255)
                                                               gACCAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCtt                                                                             >  1:10019264/1‑149 (MQ=255)
                                                                 ccAACCGGTAAATTCCAACAGTTTTTTTTTTTTTTTCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTc                                                                           >  2:8864407/1‑149 (MQ=255)
                                                                 ccAACCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTc                                                                           >  1:5044702/1‑149 (MQ=255)
                                                                   aaCCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCgg                                                                         >  1:4658944/1‑149 (MQ=255)
                                                                    aCCGGTAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGc                                                                        <  1:1625784/149‑1 (MQ=255)
                                                                       ggTAAATTCCAACAGTTTTTGGAGATTTNGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCaaa                                                                     >  1:5808458/1‑149 (MQ=255)
                                                                         tAAATTCCAACAGTTTTTTTTTTTTTTTCTATAATCTGCTCAACAATGTGCTGTTTTCTTGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACa                                                                   >  2:6774927/1‑149 (MQ=255)
                                                                         tAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACa                                                                   <  2:6499092/149‑1 (MQ=255)
                                                                         tAAATTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACa                                                                   >  2:9321498/1‑149 (MQ=255)
                                                                          aaaTTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACaa                                                                  >  1:9994447/1‑149 (MQ=255)
                                                                          aaaTTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACaa                                                                  <  1:8993614/149‑1 (MQ=255)
                                                                           aaTTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAAt                                                                 <  1:410271/149‑1 (MQ=255)
                                                                           aaTTCCAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAAt                                                                 <  1:1061439/149‑1 (MQ=255)
                                                                               ccAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAacc                                                             <  1:8435310/149‑1 (MQ=255)
                                                                               ccAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAacc                                                             <  2:1632723/149‑1 (MQ=255)
                                                                                cAACAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAacca                                                            >  1:2581577/1‑149 (MQ=255)
                                                                                   cAGTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGCTTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAAccacca                                                         >  1:6173792/1‑149 (MQ=255)
                                                                                    aggttttgGAGATTTATCCATGATCTGCTCAACAATCTGCTGTTTTCTGATTTCGTTTGATTTGACAATAGCAGTCATCAACTCCTGAGGAACGTCCTCATAACTGGCCAAAAGCTGCTTTGTGTCTTTCGGCAAACAATAACCACCAt                                                        <  2:2133666/146‑1 (MQ=255)
                                                                                     gTTTTTGGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATa                                                       >  1:5350646/1‑149 (MQ=255)
                                                                                     gTTTTTGGAGAGTTAGCCATAAGCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACGCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATa                                                       <  2:9525153/149‑1 (MQ=255)
                                                                                      ttttttGAGATTTAGCCATAATCCGCTCAACAATGTGCTGTTTTCTTGTTTCGTTTGATTTGACAATAGCCGTCATCTACCCCTGAGGACCGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATaa                                                      <  2:417701/149‑1 (MQ=255)
                                                                                          tggagaTTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCaa                                                  >  2:3674998/3‑149 (MQ=255)
                                                                                           ggagaTTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCaaa                                                 <  2:1948523/148‑1 (MQ=255)
                                                                                           ggagaTTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCaaa                                                 <  1:9591852/148‑1 (MQ=255)
                                                                                           ggagaTTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCaaa                                                 >  1:198259/2‑149 (MQ=255)
                                                                                            gagaTTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCaaaa                                                <  1:8813696/149‑1 (MQ=255)
                                                                                              gaTTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACt                                              <  2:2389696/149‑1 (MQ=255)
                                                                                                tttAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTg                                            <  2:1783958/149‑1 (MQ=255)
                                                                                                tttAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTg                                            >  1:5616317/1‑149 (MQ=255)
                                                                                                tttAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTg                                            >  2:7273044/1‑149 (MQ=255)
                                                                                                 ttAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTgg                                           >  2:2941208/1‑149 (MQ=255)
                                                                                                   aGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGat                                         >  1:1190037/1‑149 (MQ=255)
                                                                                                    gCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGatt                                        <  1:4286222/149‑1 (MQ=255)
                                                                                                      cATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGattat                                      <  2:8389210/149‑1 (MQ=255)
                                                                                                       ttAATCTGCTCCACACTGCGCTGTTTTCTGGTTTCGTTCGATTTCACAATAGCCGTCATCAACTCCTTAGGAACGCCCTCATAACTGTCCACACGCTGCTTTGAGTCTTTCGGCAAACTATTTCCTCCATAACCACACCTTGGattatt                                     <  2:5055711/148‑1 (MQ=255)
                                                                                                         aaTCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGattattat                                   <  2:1042142/149‑1 (MQ=255)
                                                                                                            cTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGattattataat                                >  2:8740171/1‑146 (MQ=255)
                                                                                                              gCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGATAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAAtgt                              >  2:6304664/1‑149 (MQ=255)
                                                                                                              gCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAAtgt                              >  2:3156188/1‑149 (MQ=255)
                                                                                                               cTCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAAtgtg                             <  1:3993274/149‑1 (MQ=255)
                                                                                                                tCAACAATGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATgtgt                            >  2:9151972/1‑149 (MQ=255)
                                                                                                                   acaaTGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCGACTCCTGAGGAACGCCCTCATAACTGCCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTccc                         >  2:5787179/1‑149 (MQ=255)
                                                                                                                   acaaTGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTccc                         <  2:9855689/149‑1 (MQ=255)
                                                                                                                   acaaTGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTccc                         <  2:1242957/149‑1 (MQ=255)
                                                                                                                    caaTGTGCTGTTTTCTGGTTTCGTTTGGTTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGCCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCa                        <  1:6361636/149‑1 (MQ=255)
                                                                                                                    caaTGTGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCa                        <  1:2168822/149‑1 (MQ=255)
                                                                                                                     aaTGGGCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTGTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCaa                       <  2:6903493/149‑1 (MQ=255)
                                                                                                                       tgtgCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCAATc                     <  1:7245828/149‑1 (MQ=255)
                                                                                                                        gtgCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCAATcc                    <  1:3058348/149‑1 (MQ=255)
                                                                                                                         tgCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCAATCCg                   <  1:6638461/149‑1 (MQ=255)
                                                                                                                          gCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCAATCCGt                  >  1:2694474/1‑149 (MQ=255)
                                                                                                                          gCTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCAATCCGt                  >  1:1835009/1‑149 (MQ=255)
                                                                                                                           cTGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAGTGTGTCCCAATCCGTg                 >  1:7024112/1‑149 (MQ=255)
                                                                                                                            tGTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGACCGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTGGAAGTATTATAATGTGTCCCAATCCGTGa                <  2:7618869/149‑1 (MQ=255)
                                                                                                                             gTTTTCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCAATCCGTGaa               >  1:6028158/1‑149 (MQ=255)
                                                                                                                              ttttCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCAATCCGTGAAt              >  2:7868670/1‑149 (MQ=255)
                                                                                                                               tttCTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCAATCCGTGAATc             >  2:3800312/1‑149 (MQ=255)
                                                                                                                                  cTGGTTTCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCAATCCGTGAATccaa          <  1:7460112/149‑1 (MQ=255)
                                                                                                                                       tacGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTCTGAGGCTTTCGGCAAACAATAACCACCATAACCAAAACTCGGATTATTATAATGAGTCCCAATACGTCAATccaacccaa     >  2:7407015/3‑149 (MQ=255)
                                                                                                                                        tCGTTTGATTTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCAATCCGTGAATccaacccaac    <  2:436357/149‑1 (MQ=255)
                                                                                                                                         cGTTTGATTTGACAATCGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATAATGTGTCCCAATCCGTGAATCCAACCCAaca   <  1:801256/149‑1 (MQ=255)
                                                                                                                                               ||                                                                                                                                              
TTGCAAGAGTTTCATCACATCAATAATCGAACTCTGTCTGAAATTATCTGACCCAGCTTTCATGACCAACCGGTAAATTCCAACAGTTTTTTGAGATTTAGCCATAATCTGCTCAACAATGTGCTGTTTTCTGGTTTCGTTTGTATTGACAATAGCCGTCATCAACTCCTGAGGAACGCCCTCATAACTGGCCAAAAGCTGCTTTGAGTCTTTCGGCAAACAATAACCACCATAACCAAAACTTGGATTATTATTATGTGTCCCAATCCGTGAATCCAACCCAACAC  >  AM406671/1588102‑1588388

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: