Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 519,975 C→T R518Q (CGG→CAG)  unknown ← Efflux ABC transporter for glutathione/L‑cysteine, essential for assembly of bd‑type respiratory oxidases => CydC subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1519,9750CT100.0% 206.9 / NA 77R518Q (CGG→CAG) unknownEfflux ABC transporter for glutathione/L‑cysteine, essential for assembly of bd‑type respiratory oxidases => CydC subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/40);  total (37/40)

TCAGAATTTTCGTGTGCCGAGATCAAGCCGGGAAGGTTTCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCCGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACTTGTTGTAGCACGGCAAGCATTTTTTCATCGCTCAACGGAGCATGATTGGCA  >  1/519734‑520217
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
tCAGAATTTTCGTGTGCCGAGATCAAGCCGGGAAGGTTTCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAAt                                                                                                                                                                                                                                             <  1:287517/249‑1 (MQ=255)
    aaTTTTCGTGTGCCGAGATCAAGCCGGGAAGGTTTCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGc                                                                                                                                                                                                                                        <  1:63894/250‑1 (MQ=255)
     aTTTTCGTGTGCCGAGATCAAGCCGGGAAGGTTTCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGt                                                                                                                                                                                                                                      >  1:538683/1‑251 (MQ=255)
             gtgCCGAGATCAAGCCGGGAAGGTTTCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGt                                                                                                                                                                                                                                      <  1:7234/243‑1 (MQ=255)
                           ccGGGAAGGTTTCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAAcc                                                                                                                                                                                                                  <  1:433003/249‑1 (MQ=255)
                           ccGGGAAGGTTTCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAAcc                                                                                                                                                                                                                  <  1:245158/249‑1 (MQ=255)
                                 aGGTTTCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCCTAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTg                                                                                                                                                                                                          >  1:432628/1‑251 (MQ=255)
                                  ggTTTCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTgg                                                                                                                                                                                                         >  1:740365/1‑251 (MQ=255)
                                   gTTTCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTgg                                                                                                                                                                                                         >  1:106383/1‑250 (MQ=255)
                                    tttCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCt                                                                                                                                                                                                       <  1:560141/251‑1 (MQ=255)
                                    tttCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCt                                                                                                                                                                                                       <  1:326763/251‑1 (MQ=255)
                                        cACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATc                                                                                                                                                                                                   <  1:747234/251‑1 (MQ=255)
                                          ccTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTa                                                                                                                                                                                                 >  1:634963/1‑251 (MQ=255)
                                                      tCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCg                                                                                                                                                                                     >  1:658129/1‑251 (MQ=255)
                                                         cAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAg                                                                                                                                                                                             <  1:755310/240‑1 (MQ=255)
                                                                  gCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTa                                                                                                                                                                         >  1:139362/1‑251 (MQ=255)
                                                                      tgtgTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGTGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAg                                                                                                                                                                     <  1:484926/251‑1 (MQ=255)
                                                                      tgtgTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAg                                                                                                                                                                     <  1:111216/251‑1 (MQ=255)
                                                                               tGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGc                                                                                                                                                            <  1:495365/251‑1 (MQ=255)
                                                                                 cTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGcc                                                                                                                                                           <  1:144053/250‑1 (MQ=255)
                                                                                          gCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAg                                                                                                                                                                     >  1:283634/1‑231 (MQ=255)
                                                                                          gCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTACCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCaaa                                                                                                                                                  >  1:409952/1‑250 (MQ=255)
                                                                                           ccATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAAc                                                                                                                                                 <  1:53396/250‑1 (MQ=255)
                                                                                                         gATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACgg                                                                                                                                               <  1:315829/238‑1 (MQ=255)
                                                                                                          aTTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCAccgc                                                                                                                                  <  1:630934/250‑1 (MQ=255)
                                                                                                              aaTCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGAt                                                                                                                             <  1:46798/251‑1 (MQ=255)
                                                                                                               aTCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACC‑TCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAg                                                                                                                                                                     >  1:839660/1‑209 (MQ=255)
                                                                                                                    gATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCAc                                                                                                                                                                                                                            <  1:188346/150‑1 (MQ=255)
                                                                                                                    gATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGga                                                                                                                       >  1:355903/1‑251 (MQ=255)
                                                                                                                          cACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGgacgac                                                                                                                   <  1:411699/249‑1 (MQ=255)
                                                                                                                           aCCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGAccg                                                                                                                 >  1:750876/1‑250 (MQ=255)
                                                                                                                              ttGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTc                                                                                                                                       <  1:806310/225‑1 (MQ=255)
                                                                                                                                gTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCAc                                                                                                                                                                                                                            >  1:299648/1‑138 (MQ=255)
                                                                                                                                   cAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCg                                                                                                                                                                                                                                       >  1:570676/1‑124 (MQ=255)
                                                                                                                                      ttaattGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAAc                                                                                                     <  1:199933/251‑1 (MQ=255)
                                                                                                                                        aattGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGc                                                                                                                                                                    <  1:351813/186‑1 (MQ=255)
                                                                                                                                            gattgTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATaaa                                                                                               >  1:765268/1‑251 (MQ=255)
                                                                                                                                              ttgTTATCAATCATGACCAACTGATAAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATc                                                                                             <  1:496397/251‑1 (MQ=255)
                                                                                                                                                gTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATc                                                                                             >  1:548425/1‑249 (MQ=255)
                                                                                                                                                 ttATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCaa                                                                                          >  1:878886/1‑251 (MQ=255)
                                                                                                                                                    tcaatcaTGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGCTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGcc                                                                                       <  1:252654/251‑1 (MQ=255)
                                                                                                                                                      aatcaTGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATc                                                                                             >  1:104834/1‑243 (MQ=255)
                                                                                                                                                       atcaTGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCg                                                                                                        >  1:109573/1‑231 (MQ=255)
                                                                                                                                                       atcaTGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTg                                                                                    >  1:81708/1‑251 (MQ=255)
                                                                                                                                                       atcaTGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTg                                                                                    >  1:51995/1‑251 (MQ=255)
                                                                                                                                                            gACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATc                                                                                             >  1:174280/1‑237 (MQ=255)
                                                                                                                                                               cAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGAt                                                                            <  1:257950/251‑1 (MQ=255)
                                                                                                                                                                    gATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAAcc                                                                                                                                                                                                                  >  1:564727/1‑112 (MQ=255)
                                                                                                                                                                    gATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTa                                                                                                                                                                                                 >  1:527556/1‑129 (MQ=255)
                                                                                                                                                                    gATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCAcc                                                                                                                                                                 <  1:889095/161‑1 (MQ=255)
                                                                                                                                                                    gATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAAc                                                                                                                                                 >  1:188213/1‑177 (MQ=255)
                                                                                                                                                                    gATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACg                                                                       >  1:234365/1‑251 (MQ=255)
                                                                                                                                                                    gATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACg                                                                       >  1:818665/1‑251 (MQ=255)
                                                                                                                                                                    gATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCACAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAg                                                                                                                                                                     <  1:391533/157‑1 (MQ=255)
                                                                                                                                                                           ttGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCAccgc                                                                                                                                  >  1:447903/1‑185 (MQ=255)
                                                                                                                                                                             gTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTc                                                               >  1:43663/1‑250 (MQ=255)
                                                                                                                                                                                      cGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCa                                                       >  1:507936/1‑249 (MQ=255)
                                                                                                                                                                                      cGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCa                                                       >  1:152015/1‑249 (MQ=255)
                                                                                                                                                                                        gTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACtt                                                    <  1:12685/250‑1 (MQ=255)
                                                                                                                                                                                         ttAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACttg                                                   <  1:597605/250‑1 (MQ=255)
                                                                                                                                                                                           aaGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACttgttg                                                <  1:605724/251‑1 (MQ=255)
                                                                                                                                                                                            aGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTc                                                                                 <  1:840663/217‑1 (MQ=255)
                                                                                                                                                                                                 tGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACTTGTTGTAGCa                                           <  1:72880/250‑1 (MQ=255)
                                                                                                                                                                                                     tCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAACCCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCCCCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACTTGTTGTAGCa                                           <  1:115009/246‑1 (MQ=255)
                                                                                                                                                                                                        ccATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCgg                                                                                                                                                                                    <  1:93484/106‑1 (MQ=255)
                                                                                                                                                                                                            aattaatGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGTATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCt                                                                                                                                                                                                       >  1:689948/1‑83 (MQ=255)
                                                                                                                                                                                                            aattaatGTTTTCCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAAt                                                                                                                                                                       <  1:249500/115‑1 (MQ=255)
                                                                                                                                                                                                                aatGTTTTTCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACTTGTTGTAGCACGGCAAGCATTTTTTc                           <  1:884153/251‑1 (MQ=255)
                                                                                                                                                                                                                       ttCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAg                                                                                                                                                                                             <  1:893039/82‑1 (MQ=255)
                                                                                                                                                                                                                       ttCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATcc                                                                                            <  1:302664/179‑1 (MQ=255)
                                                                                                                                                                                                                       ttCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGc                                                                                        >  1:456847/1‑183 (MQ=255)
                                                                                                                                                                                                                       ttCCACGGCTGTGTGCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACttg                                                   >  1:808668/1‑220 (MQ=255)
                                                                                                                                                                                                                                   gtgCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACTTGTTGTAGCACGGCAAGCATTTTTTCATCGCTCAACGGAGCATGa        >  1:21455/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gtgCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACTTGTTGTAGCACGGCAAGCATTTTTTCATCGCTCAACGGAGCATGa        >  1:85137/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gtgCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACTTGTTGTAGCACGGCAAGCATTTTTTCATCGCTCAACGGAGCATGa        >  1:46/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tgCCAAAATTGCCTGCAGAATGTGGCGTTCCGTTTCGCGATCCAACCCTTCCGTGGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCCCCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACTTGTTGTAGCACGGCAAGCATTTTTTCATCGCTCAACGGAGCATGAt       <  1:669603/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaTTGCCTGCAGAATGTGGCGTTCCGTTTCACGATCCAACCCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACTTGTTGTAGCACGGCAAGCATTTTTTCATCGCTCAACGGAGCATGATTGGCa  <  1:343260/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  
TCAGAATTTTCGTGTGCCGAGATCAAGCCGGGAAGGTTTCCACCTTTCCAGCCTTCCCAATTAAATGCGTTGTGTCAATTGCTTGAAAAAGCCATTTTCTTTACTGATTAAATCCAGATAAGCACCTTGTTCAATTAATTGATTGTTATCAATCATGACCAACTGATCAAATTGTTCAATAGCGGTTAAGCGGTGCGTCACCATAATTAATGTTTTTCCACGGCTGTGTGCCAAAATTGCCCGCAGAATGTGGCGTTCCGTTTCACGATCCAAACCTTCCGTTGGCTCATCTAACAGTAAAATCGGTGCGTCATTTAATAGCACCCGAGCCAAGCCCAAACGGCGTTGTTCACCGCCGGATAACGGACGACCGCCGTCACCGAACCATAAATCCAAGCCTTGTTCCTGATTCACGAGTTTTTCCAAGCCCACTTGTTGTAGCACGGCAAGCATTTTTTCATCGCTCAACGGAGCATGATTGGCA  >  1/519734‑520217

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