Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,870,135 G→A P446L (CCG→CTG)  unknown ← ATP‑dependent DNA helicase UvrD/PcrA (EC 3.6.4.12)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,870,1350GA100.0% 191.3 / NA 66P446L (CCG→CTG) unknownATP‑dependent DNA helicase UvrD/PcrA (EC 3.6.4.12)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/31);  total (35/31)

TCTTCAGGTTTGATAAATTCACGGGCGGCAGACACCAATTCTTCCAGGTTTTCAATACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCGGCGCGCCCGGCGAGTTTGTTTTCCTGAAGCGCCGGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATCTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCATAAATGCGATACGGAATTTGCGAACGAATCAGCGCTTCTTCAATCACT  >  1/1869892‑1870380
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
tCTTCAGGTTTGATAAATTCACGGGCGGCAGACACCAATTCTTCCAGGTTTTCAATACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATtggt                                                                                                                                                                                                                                                 >  1:847941/1‑251 (MQ=255)
              aaaTTCACGGGCGGCAGACACCAATTCTTCCAGGTTTTCAATACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAgg                                                                                                                                                                                                                                   <  1:594602/251‑1 (MQ=255)
                       ggcggcAGACACCAATTCTTCCAGGTTTTCAATACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGAt                                                                                                                                                                                                                          >  1:753859/1‑251 (MQ=255)
                           gcAGACACCAATTCTTCCAGGTTTTCAATACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGtt                                                                                                                                                                                                                       >  1:40844/1‑250 (MQ=255)
                                cacCAATTCTTCCAGGTTTTCAATACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATtggtgg                                                                                                                                                                                                                                               <  1:129479/221‑1 (MQ=255)
                                           ccAGGTTTTCAATACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTgcgc                                                                                                                                                                                                      <  1:427804/251‑1 (MQ=255)
                                                    cAATACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATg                                                                                                                                                                                                                         <  1:835388/223‑1 (MQ=255)
                                                    cAATACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATcc                                                                                                                                                                                             <  1:460408/251‑1 (MQ=255)
                                                      aTACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCaa                                                                                                                                                                                           <  1:183427/251‑1 (MQ=255)
                                                      aTACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCaa                                                                                                                                                                                           >  1:353191/1‑251 (MQ=255)
                                                       tACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑c                                                                                                                                                                                         <  1:630018/251‑1 (MQ=255)
                                                             ccTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGt                                                                                                                                                                                                 >  1:722176/1‑238 (MQ=255)
                                                                   cTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑Cg                                                                                                                                                                                        >  1:807027/1‑240 (MQ=255)
                                                                           cccTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGa                                                                                                                                                                                                                           <  1:878331/198‑1 (MQ=255)
                                                                                   cTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCg                                                                                                                                                             <  1:79071/251‑1 (MQ=255)
                                                                                     gTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTa                                                                                                                                                           >  1:427438/1‑251 (MQ=255)
                                                                                           aCATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCGTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGGTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATg                                                                                                                                                     <  1:840868/251‑1 (MQ=255)
                                                                                              tttGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATCAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGAcc                                                                                                                                                  <  1:435083/251‑1 (MQ=255)
                                                                                                 gATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGAt                                                                                                                                                                                                           >  1:758745/1‑192 (MQ=255)
                                                                                                  aTACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCg                                                                                                                                                >  1:418074/1‑249 (MQ=255)
                                                                                                       aCCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCt                                                                                                                                              >  1:364086/1‑246 (MQ=255)
                                                                                                       aCCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAg                                                                                                                                         >  1:567281/1‑251 (MQ=255)
                                                                                                         ccGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGgc                                                                                                                                       <  1:373549/251‑1 (MQ=255)
                                                                                                               tGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATg                                                                                                                                                                                                                         <  1:192281/164‑1 (MQ=255)
                                                                                                                 cTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCAt                                                                                                                                  >  1:231379/1‑248 (MQ=255)
                                                                                                                       cacaAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTc                                                                                                                              >  1:349504/1‑246 (MQ=255)
                                                                                                                        acaAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCg                                                                                                                        <  1:98879/251‑1 (MQ=255)
                                                                                                                                          aacaacGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAAc                                                                                                           >  1:500166/1‑246 (MQ=255)
                                                                                                                                            caacGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAAt                                                                                                    <  1:37712/251‑1 (MQ=255)
                                                                                                                                                 gCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAgg                                                                                                                                                                                                                                   <  1:867139/120‑1 (MQ=255)
                                                                                                                                                 gCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGc                                                                                                >  1:323250/1‑250 (MQ=255)
                                                                                                                                                         tcCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGa                                                                                                                       >  1:437205/1‑219 (MQ=255)
                                                                                                                                                              gtCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTaa                                                                                  >  1:183582/1‑251 (MQ=255)
                                                                                                                                                                      gCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCg                                                                                                                        <  1:710104/205‑1 (MQ=255)
                                                                                                                                                                           gAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGc                                                                                          >  1:594951/1‑230 (MQ=255)
                                                                                                                                                                           gAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGtt                                                                      >  1:440010/1‑250 (MQ=255)
                                                                                                                                                                           gAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGtt                                                                      >  1:36794/1‑250 (MQ=255)
                                                                                                                                                                           gAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGtt                                                                      >  1:563893/1‑250 (MQ=255)
                                                                                                                                                                                aattaaCTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAACCGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAAt                                                                                                    <  1:546034/216‑1 (MQ=255)
                                                                                                                                                                                aattaaCTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAaaga                                                                >  1:1026/1‑251 (MQ=255)
                                                                                                                                                                                aattaaCTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAaaga                                                                <  1:740229/251‑1 (MQ=255)
                                                                                                                                                                                aattaaCTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAaaga                                                                <  1:181327/251‑1 (MQ=255)
                                                                                                                                                                                aattaaCTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTAGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAaaga                                                                >  1:475482/1‑251 (MQ=255)
                                                                                                                                                                                   taaCTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGAt                                                                                                                                                                                                           >  1:394952/1‑110 (MQ=255)
                                                                                                                                                                                     aCTCCATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAAt                                                                                                                                                                           >  1:10060/1‑139 (MQ=255)
                                                                                                                                                                                        ccATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTTAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAAc                                                                                                           <  1:10284/200‑1 (MQ=255)
                                                                                                                                                                                        ccATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGAcc                                                                                                                                                  >  1:820927/1‑161 (MQ=255)
                                                                                                                                                                                        ccATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGcc                                                                                               >  1:817124/1‑212 (MQ=255)
                                                                                                                                                                                        ccATAAAGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAAAGAAGCGCAt                                                         >  1:617846/1‑250 (MQ=255)
                                                                                                                                                                                             aaGCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGc                                                                                          <  1:27885/212‑1 (MQ=255)
                                                                                                                                                                                               gCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTAt                                                                                                                                                          >  1:256989/1‑146 (MQ=255)
                                                                                                                                                                                               gCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAg                                                                                                                                         >  1:13307/1‑163 (MQ=255)
                                                                                                                                                                                               gCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGgc                                                                                                                                       <  1:223491/165‑1 (MQ=255)
                                                                                                                                                                                               gCGTAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCAt                                                 >  1:342275/1‑251 (MQ=255)
                                                                                                                                                                                                  tAACAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCATaa                                               <  1:383298/250‑1 (MQ=255)
                                                                                                                                                                                                     cAATGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCATAAATGc                                           <  1:63751/251‑1 (MQ=255)
                                                                                                                                                                                                        tGCCGTAGCAGCGCGCCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCATAAATGCGAt                                        >  1:710622/1‑251 (MQ=255)
                                                                                                                                                                                                              agcagcgcgcCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCg                                                                                                                        <  1:840591/161‑1 (MQ=255)
                                                                                                                                                                                                                     cgcCAGGCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCATAAATGCGATACGGAATTTGcg                            <  1:573583/250‑1 (MQ=255)
                                                                                                                                                                                                                       ccaggCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGAcc                                                                                                                                                  <  1:448066/127‑1 (MQ=255)
                                                                                                                                                                                                                       ccaggCGAGTTTGTTTTCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCATAAATGCGATACGGAATTTGcgaa                          <  1:338972/247‑1 (MQ=255)
                                                                                                                                                                                                                                      ttCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTgg                                                                                                                                                               <  1:650721/102‑1 (MQ=255)
                                                                                                                                                                                                                                      ttCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCATAAATGCGATACGGAATTTGCGAACGAATCAGCGcttctt          >  1:623435/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ttCCTGAAGCGCCAGATTGGTGGCTTGCCATAAGGGGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAG‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCATAAATGCGATACGGAATTTGCGAACGAATCAGCGcttctt          >  1:902020/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaGCGCCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCATAAATGCGATACGGAATTTGCGAACGAATCAGCGCTTCTTCAATCa    <  1:29099/251‑1 (MQ=255)
                                                                                                                                                                                                                                               cgcCAGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATTTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCATAAATGCGATACGGAATTTGCGAACGAATCAGCGCTTCTTCAATCACt  >  1:469069/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
TCTTCAGGTTTGATAAATTCACGGGCGGCAGACACCAATTCTTCCAGGTTTTCAATACGCACCTCCCCTTTTTCGCCCTTTTCCTGTTTGTACATTTGATACAACCCGGAATGCTTAATCACAAAATCCGTTTGCGCAAACAACGGCATTTCTTCCGTGTCCCGCTGCAAAGAATTAATTAACTCCATAAAGCGTAACAATGCCGTAGCGGCGCGCCCGGCGAGTTTGTTTTCCTGAAGCGCCGGATTGGTGGCTTGCCATAAGGTGATTTGATGTTCACGGGTAAGATTGCGCAAGGTATCCAA‑CGTGCGATCGCCAATGCCGCGCGGTGGCGTATTAATGACCCGCTCAAAGGCAGCATCGTCCTGTCGATTGGCGATTAACCGTAAATAAGCCAAGGCATCTTTAATTTCCTGCCGTTCAAAGAAGCGCATCCCACCATAAATGCGATACGGAATTTGCGAACGAATCAGCGCTTCTTCAATCACT  >  1/1869892‑1870380

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