Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,548,618 A→G K12E (AAA→GAA)  unknown → Exopolyphosphatase (EC 3.6.1.11)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,548,6180AG100.0% 144.6 / NA 53K12E (AAA→GAA) unknownExopolyphosphatase (EC 3.6.1.11)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/23);  total (30/23)

AACCTGAGTACCGCCGGGGTGGGCGACAGATCTAAAGTAGTGGTCAGCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACAAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAACGCGGGGTAAATTGTCTGGCGTTATTTGCCGAACGCTTGCAGGGATTTTC  >  1/1548376‑1548853
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
aaCCTGAGTACCGCCGGGGTGGGCGACAGATCTAAAGTAGTGGTCAGCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTAc                                                                                                                                                                                                                                     <  1:693572/251‑1 (MQ=255)
   cTGAGTACCGCCGGGGTGGGCGACAGATCTAAAGTAGTGGTCAGCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGcc                                                                                                                                                                                                                                  <  1:791506/251‑1 (MQ=255)
                 ggTGGGCGACAGATCTAAAGTAGTGGTCAGCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAgg                                                                                                                                                                                                                     <  1:191791/250‑1 (MQ=255)
                            gATCTAAAGTAGTGGTCAGCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGa                                                                                                                                                                                                                       <  1:706841/237‑1 (MQ=255)
                            gATCTAAAGTAGTGGTCAGCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAAt                                                                                                                                                                                                           >  1:476198/1‑249 (MQ=255)
                            gATCTAAAGTAGTGGTCAGCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATcg                                                                                                                                                                                                         >  1:27136/1‑251 (MQ=255)
                             aTCTAAAGTAGTGGTCAGCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAgg                                                                                                                                                                                                                     >  1:86982/1‑238 (MQ=255)
                              tCTAAAGTAGTGGTCAGCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTAGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATcgc                                                                                                                                                                                                        >  1:482997/1‑250 (MQ=255)
                                          gTCAGCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACtt                                                                                                                                                                                            <  1:279680/250‑1 (MQ=255)
                                              gCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGtt                                                                                                                                                                                       >  1:715968/1‑251 (MQ=255)
                                                    ccATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGaa                                                                                                                                                                                                             <  1:431358/223‑1 (MQ=255)
                                                    ccATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAAt                                                                                                                                                                                                           <  1:161022/225‑1 (MQ=255)
                                                    ccATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAg                                                                                                                                                                                          >  1:603162/1‑242 (MQ=255)
                                                    ccATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGt                                                                                                                                                                                        >  1:75268/1‑244 (MQ=255)
                                                    ccATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACa                                                                                                                                                                                 >  1:222203/1‑251 (MQ=255)
                                                    ccATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACa                                                                                                                                                                                 >  1:376152/1‑251 (MQ=255)
                                                                     cAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGatcat                                                                                                                                                                 >  1:523259/1‑250 (MQ=255)
                                                                            aaTACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCa                                                                                                                                                         >  1:135925/1‑251 (MQ=255)
                                                                                  gCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAg                                                                                                                                                                                >  1:322434/1‑222 (MQ=255)
                                                                                            aTTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGTTTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGtt                                                                                                                                         >  1:41533/1‑251 (MQ=255)
                                                                                                 aCGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAgg                                                                                                                                    <  1:286985/251‑1 (MQ=255)
                                                                                                 aCGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAg                                                                                                                                     <  1:650356/250‑1 (MQ=255)
                                                                                                    gTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAgg                                                                                                                                                                                                                     >  1:371580/1‑167 (MQ=255)
                                                                                                    gTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGtat                                                                                                                                 >  1:892259/1‑251 (MQ=255)
                                                                                                         gCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATcg                                                                                                                            >  1:673787/1‑251 (MQ=255)
                                                                                                               gTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAg                                                                                                                                                                                >  1:848254/1‑193 (MQ=255)
                                                                                                                                     tGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCgg                                                                                                >  1:633641/1‑251 (MQ=255)
                                                                                                                                       ccACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGcc                                                                                                                                                             <  1:367008/188‑1 (MQ=255)
                                                                                                                                          caACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAgg                                                                                            >  1:490791/1‑250 (MQ=255)
                                                                                                                                            aCTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCtt                                                                                         >  1:394324/1‑251 (MQ=255)
                                                                                                                                              tGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTgg                                                                                       >  1:258407/1‑251 (MQ=255)
                                                                                                                                               gTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGa                                                                                      >  1:261859/1‑251 (MQ=255)
                                                                                                                                                                  ggTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAg                                                                           >  1:241439/1‑243 (MQ=255)
                                                                                                                                                                   gtgtTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGc                                                                   >  1:696075/1‑250 (MQ=255)
                                                                                                                                                                   gtgtTAATCAAGCCTCTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAgg                                                                                                                                    <  1:455968/185‑1 (MQ=255)
                                                                                                                                                                     gtTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGta                                                                                                                                  >  1:300721/1‑185 (MQ=255)
                                                                                                                                                                         aTCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTATGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGgcag                                                            <  1:438926/251‑1 (MQ=255)
                                                                                                                                                                             aGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATcg                                                                                                                            <  1:612809/183‑1 (MQ=255)
                                                                                                                                                                               ccTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATcg                                                                                                                            <  1:137281/181‑1 (MQ=255)
                                                                                                                                                                                 tttGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAActc                                                                                                                                                                                                                                        >  1:496648/1‑71 (MQ=255)
                                                                                                                                                                                 tttGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGaa                                                    <  1:675261/251‑1 (MQ=255)
                                                                                                                                                                                  ttGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGaa                                                    >  1:621056/1‑250 (MQ=255)
                                                                                                                                                                                     aTCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCATATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAAcgc                                                 <  1:874031/250‑1 (MQ=255)
                                                                                                                                                                                     aTCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAAcgc                                                 <  1:380537/250‑1 (MQ=255)
                                                                                                                                                                                         aCTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGatc                                                                                                                                                                   <  1:768754/132‑1 (MQ=255)
                                                                                                                                                                                               tCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAACGCGGGGTAAATTg                                      <  1:298909/251‑1 (MQ=255)
                                                                                                                                                                                                      gAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAACGCGGGGTAAATTGTCTGGCg                               >  1:348537/1‑251 (MQ=255)
                                                                                                                                                                                                       aaTTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAACGCGGGGTAAATTGTCTGGCg                               <  1:101872/250‑1 (MQ=255)
                                                                                                                                                                                                                   gAATAACGAAAATTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAACGCGGGGTAAATTGTCTGGCGTTATTTGCCGaa                   >  1:593186/1‑250 (MQ=255)
                                                                                                                                                                                                                           aaaaTTTACTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAACGCGGGGTAAATTGTCTGGCGTTATTTGCCGAACGCTTGCAg          <  1:518992/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cTAAAGAAAGCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAACGCGGGGTAAATTGTCTGGCGTTATTTGCCGAACGCTTGCAGGGATTTTc  <  1:22966/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aaagaaagCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAACGCGGGGTAAATTg                                      <  1:234388/213‑1 (MQ=255)
                                                                                                                                                                                                                                     aaagaaagCTAACGAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAACGCGGGGTAAATTGt                                     >  1:850695/1‑214 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                           
AACCTGAGTACCGCCGGGGTGGGCGACAGATCTAAAGTAGTGGTCAGCGGCACCATTACACCGCGCTTACAGGTAAAATACGGCGTTGGGTTATTTGACGGTTTAGCCGAAGTGACTCTGCGTTACCGATTATTGCCACAACTGTATGTTCAATCCGTTTCGGGTGTTAATCAAGCCTTTGATCTACTCTATCAATTTGAATTTTAATTATGAATAACGAAAATTTACTAAAGAAAGCTAACAAACTCTACGCCCGACGTGGCGAGGTGCGTGAAATCGCCGCGATCGACTTAGGTTCCAACAGCTTTCACATGATCATCGCCCGCATCGTCAACGGCACGTTACAGGTATTATCGCGGTTGAAGCAGAAAGTCCAGCTTGCGGAAGGCTTGGATGAAAACCAAGTGTTAAGCCAGGCAGCAATTGAACGCGGGGTAAATTGTCTGGCGTTATTTGCCGAACGCTTGCAGGGATTTTC  >  1/1548376‑1548853

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