Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 2,062,549 G→A H282Y (CAT→TAT)  unknown ← Phosphoserine phosphatase (EC 3.1.3.3)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*12,062,5490GA100.0% 213.7 / NA 74H282Y (CAT→TAT) unknownPhosphoserine phosphatase (EC 3.1.3.3)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/40);  total (34/40)

AAATTCAATGGCGGCTTGTACGTTGCGGAAGTCCCAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATGATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCGTCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCAAGCCGAGCTGTTGCTGTAACACTTCGGCGA  >  1/2062315‑2062788
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
aaaTTCAATGGCGGCTTGTACGTTGCGGAAGTCCCAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCaaa                                                                                                                                                                                                                                 >  1:588769/1‑251 (MQ=255)
  aTTCAATGGCGGCTTGTACGTTGCGGAAGTCCCAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAg                                                                                                                                                                                                                                             <  1:177954/237‑1 (MQ=255)
  aTTCAATGGCGGCTTGTACGTTGCGGAAGTCCCAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAAt                                                                                                                                                                                                                                >  1:670915/1‑250 (MQ=255)
            ggcTTGTACGTTGCGGAAGTCCCAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCGCAAATCGGCGACCTg                                                                                                                                                                                                                      >  1:642651/1‑250 (MQ=255)
            ggcTTGTACGTTGCGGAAGTCCCAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTg                                                                                                                                                                                                                      >  1:121103/1‑250 (MQ=255)
             gcTTGTACGTTGCGGAAGTCCCAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTggg                                                                                                                                                                                                                    >  1:569973/1‑249 (MQ=255)
              cTTGTACGTTGCGGAAGTCCCAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGcat                                                                                                                                                                                                                   <  1:297347/251‑1 (MQ=255)
                tGTACGTTGCGGAAGTCCCAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGcatca                                                                                                                                                                                                                 <  1:690014/251‑1 (MQ=255)
                    cGTTGCGGAAGTCCCAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATgg                                                                                                                                                                                                              >  1:801861/1‑250 (MQ=255)
                                cccAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTgg                                                                                                                                                                                                  >  1:287509/1‑250 (MQ=255)
                                                      aTTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGgtg                                                                                                                                                                           >  1:228348/1‑251 (MQ=255)
                                                              tttCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTcggc                                                                                                                                                                   <  1:312086/251‑1 (MQ=255)
                                                               ttCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGa                                                                                                                                                                                                         <  1:4199/212‑1 (MQ=255)
                                                                      gCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAg                                                                                                                                                                                                                                             <  1:273615/169‑1 (MQ=255)
                                                                                   cGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTgcg                                                                                                                                               >  1:743688/1‑250 (MQ=255)
                                                                                        ccGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAgtt                                                                                                                                         >  1:334225/1‑251 (MQ=255)
                                                                                             gATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGc                                                                                                                                                                                                             <  1:702258/178‑1 (MQ=255)
                                                                                              aTTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTa                                                                                                                                   <  1:626820/251‑1 (MQ=255)
                                                                                                 tCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGa                                                                                                                                <  1:863024/251‑1 (MQ=255)
                                                                                                 tCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGa                                                                                                                                <  1:482459/251‑1 (MQ=255)
                                                                                                           tGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCttt                                                                                                                      >  1:704487/1‑251 (MQ=255)
                                                                                                            gTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGAt                                                                                                                                                                                                        <  1:213106/168‑1 (MQ=255)
                                                                                                                aaTGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTAt                                                                                                                  <  1:384283/250‑1 (MQ=255)
                                                                                                                aaTGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTAt                                                                                                                  <  1:863420/250‑1 (MQ=255)
                                                                                                                    atgGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGc                                                                                                             <  1:23725/251‑1 (MQ=255)
                                                                                                                      ggCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGa                                                                                                                                >  1:478745/1‑230 (MQ=255)
                                                                                                                      ggCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGt                                                                                                           >  1:705410/1‑251 (MQ=255)
                                                                                                                           ggTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGc                                                                                                                                                                                                             >  1:458183/1‑148 (MQ=255)
                                                                                                                            gTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCggcggc                                                                                                                                                                >  1:571822/1‑192 (MQ=255)
                                                                                                                            gTTTACTCCTTATTTGATTTTGTCG‑TAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACaa                                                                                                     <  1:10150/250‑1 (MQ=255)
                                                                                                                               tACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATgg                                                                                                  <  1:132801/251‑1 (MQ=255)
                                                                                                                                   ccTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCAc                                                                                               >  1:152726/1‑250 (MQ=255)
                                                                                                                                          tGATTTTGTCGTTAGCTCTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACg                                                                                       <  1:378381/251‑1 (MQ=255)
                                                                                                                                          tGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACg                                                                                       <  1:198925/251‑1 (MQ=255)
                                                                                                                                                     ttAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAg                                                                            <  1:531872/251‑1 (MQ=255)
                                                                                                                                                        gCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTcgg                                                                                                                                                                    >  1:340633/1‑160 (MQ=255)
                                                                                                                                                           cTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAACTTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGc                                                                       >  1:443605/1‑250 (MQ=255)
                                                                                                                                                             tAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGccat                                                                    <  1:114313/251‑3 (MQ=255)
                                                                                                                                                                 aTACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGAt                                                                                                                               <  1:693058/188‑1 (MQ=255)
                                                                                                                                                                 aTACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCAc                                                                                               <  1:65175/220‑1 (MQ=255)
                                                                                                                                                                 aTACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCgat                                                                >  1:58602/1‑251 (MQ=255)
                                                                                                                                                                   aCATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAAtt                                                                                                                                                     >  1:696801/1‑164 (MQ=255)
                                                                                                                                                                   aCATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTAt                                                                                                                  <  1:554993/199‑1 (MQ=255)
                                                                                                                                                                   aCATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACg                                                                                       <  1:249494/226‑1 (MQ=255)
                                                                                                                                                                   aCATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAg                                                                            <  1:868925/237‑1 (MQ=255)
                                                                                                                                                                   aCATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCgatga                                                              >  1:824736/1‑251 (MQ=255)
                                                                                                                                                                   aCATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCgatga                                                              >  1:782340/1‑251 (MQ=255)
                                                                                                                                                                   aCATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTCGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTAc                                                                                                       <  1:424311/210‑1 (MQ=255)
                                                                                                                                                                    cATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGa                                                                                                                                <  1:834593/184‑1 (MQ=255)
                                                                                                                                                                            tGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAAt                                                      <  1:248083/250‑1 (MQ=255)
                                                                                                                                                                            tGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAAt                                                      <  1:351297/250‑1 (MQ=255)
                                                                                                                                                                              cGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTAt                                                                                                                  >  1:155817/1‑188 (MQ=255)
                                                                                                                                                                                gTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGAt                                                                                                                                                                                                        <  1:431350/100‑1 (MQ=255)
                                                                                                                                                                                 ttAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAAttggtt                                                 <  1:317516/250‑1 (MQ=255)
                                                                                                                                                                                   aaaTCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGc                                              <  1:195244/251‑1 (MQ=255)
                                                                                                                                                                                         gTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCAc                                                                                               <  1:613637/196‑1 (MQ=255)
                                                                                                                                                                                              aTTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTAc                                                                                        >  1:426325/1‑198 (MQ=255)
                                                                                                                                                                                                 gatgatAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGa                                                                                                                                <  1:485611/155‑1 (MQ=255)
                                                                                                                                                                                                 gatgatAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTAt                                                                                                                  <  1:511552/169‑1 (MQ=255)
                                                                                                                                                                                                 gatgatAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCa                                >  1:434651/1‑251 (MQ=255)
                                                                                                                                                                                                 gatgatAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCa                                >  1:224017/1‑251 (MQ=255)
                                                                                                                                                                                                  atgatAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAgtt                                                                                                                                         >  1:337469/1‑145 (MQ=255)
                                                                                                                                                                                                         cTTGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGACAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCAAGCCGAg                         >  1:415998/1‑250 (MQ=255)
                                                                                                                                                                                                           tGCGCTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCAAGCCGAg                         >  1:708484/1‑248 (MQ=255)
                                                                                                                                                                                                            gcgcTAATTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCAAGCCGAg                         >  1:723413/1‑247 (MQ=255)
                                                                                                                                                                                                                  aTTGTTGCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCAAGCCGAGCTGTTGCTGt               <  1:127256/251‑1 (MQ=255)
                                                                                                                                                                                                                   ttgttgCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATcc                                 >  1:650956/1‑232 (MQ=255)
                                                                                                                                                                                                                   ttgttgCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCAAGCCGAGCTGTTGCTGTa              >  1:585935/1‑251 (MQ=255)
                                                                                                                                                                                                                   ttgttgCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCAAGCCGAGCTGTTGCTGTa              >  1:852178/1‑251 (MQ=255)
                                                                                                                                                                                                                       tgCACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCAAGCCGAGCTGTTGCTGTAacac          <  1:674314/251‑1 (MQ=255)
                                                                                                                                                                                                                       tgCACTTTCGATTTGGCATAATGGGCAGCGCCCAAATCGGCGGCCTGCAGCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCGACGGTGTTTTCGGCGGCAATGTGAAATTGTTACGCAAGTTGTTGTAAGTTTTGCGCTTTGTATTGTGCATCGACAATGACACCTTGTACGTTGCCGTTAAGTTTGCCATCGATGATTTCAAATTGGTTAGCAGCGATGAAATCCAAGCCGAGCTGTTGCTGTAacac          <  1:142934/251‑1 (MQ=255)
                                                                                                                                                                                                                         cACTTTCGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCAAGCCGAGCTGTTGCTGTAACACtt        <  1:456663/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGGTTTGGCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCATCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCAAGCCGAGCTGTTGCTGTAACACTTCGGCGa  <  1:434926/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gCATAATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCAc                                                                                               >  1:100634/1‑151 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                               
AAATTCAATGGCGGCTTGTACGTTGCGGAAGTCCCAGCGGTTGGTCAGCTCGCGATTGGCGTTTTCCACGGCGTTGCGGATTTCGTGTCCGCTGATTTCGGAAACAATGTCGAATGATGGCATGGTTTACTCCTTATTTGATTTTGTCGTTAGCACTTAAAATACATAAAAGTGCGGTTAAATCGGTGAAATTGATGATAACTTGCGCTAATTGTTGCACTTTCGGTTTGGCATGATAGGCAACGCCCAAATCGGCGACCTGCATCATGGCAAGATCGTTGGCGCCGTCGCCGATGGCAACGGTGTTTTCGGCGGCAATGTGAAATTGTTGCGCAAGTTGTTGTAAGATTTGCGCTTTGTATTGTGCGTCTACAATGGCACCTTGTACGTTGCCGGTAAGTTTGCCGTCGATGATTTCAAATTGGTTGGCAGCGATGAAATCCAAGCCGAGCTGTTGCTGTAACACTTCGGCGA  >  1/2062315‑2062788

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