Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,481,792 2 bp→AG coding (398‑399/504 nt) unknown → Flp pilus assembly protein RcpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,481,7920CA100.0% 154.3 / NA 55P133Q (CCC→CAG) unknownFlp pilus assembly protein RcpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (24/31);  total (24/31)
*11,481,7930CG100.0% 149.3 / NA 55P133Q (CCC→CAG) unknownFlp pilus assembly protein RcpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (24/31);  total (24/31)

ATGTCAGTTATAGTAATTTATTTAATTCCGTTATTCGTGATATAGGTAGTAACAACCAAAAGCAAATCAATATAATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGAAAATTCGTACAAAATTGATTAAAAATGGTATTACCTCTAAATCTATTTACTTAGAGAAAAAAGCACATAAAGCAAAAACATATCCAATCTATATTGAAGTAGCTCAAATTGCGCCCAAACCGGCACCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCTCGCACTATTACTATAGTGTCAAGCCGTGATGATATCCAAGGTGAAGTTGCACAAATTCTACGTACT  >  1/1481567‑1482021
                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                    
aTGTCAGTTATAGTAATTTATTCAATTCTGTCATTCGTGATATAGGTAGTAACAACCAAAAACAAATCAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGt                                                                                                                                                                                                              <  1:1760/251‑1 (MQ=255)
         ataGTAATTTATTCAATTCTGTCATTCGTGATATAGGTAGTAACAACCAAAAACAAATCAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCCTGTGCATATCGTACCGGAGaa                                                                                                                                                                                                     >  1:475489/1‑251 (MQ=255)
              aaTTTATTCAATTCTGTCATTCGTGATATAGGTAGTAACAACCAAAAACAAATCAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAagtagt                                                                                                                                                                                                >  1:696406/1‑251 (MQ=255)
                        attctgtcaTTCGTGATATAGGTAGTAACAACCAAAAACAAATCAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGt                                                                                                                                                                                                              <  1:259751/219‑1 (MQ=255)
                                     tGATATAGGTAGTAACAACCAAAAACAAATCAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATc                                                                                                                                                                                                                >  1:680285/1‑212 (MQ=255)
                                     tGATATAGGTAGTAACAACCAAAAACAAATCAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGActt                                                                                                                                                                         >  1:751871/1‑251 (MQ=255)
                                         ataGGTAGTAACAACCAAAAACAAATCAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTtg                                                                                                                                                                     <  1:385612/251‑1 (MQ=255)
                                            ggTAGTAACAACCAAAAACAAATCAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGc                                                                                                                                                                                                                     <  1:27961/200‑1 (MQ=255)
                                            ggTAGTAACAACCAAAAACAAATCAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTtgtg                                                                                                                                                                   >  1:46376/1‑250 (MQ=255)
                                            ggTAGTAACAACCAAAAACAAATCAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTtgtg                                                                                                                                                                   >  1:403827/1‑250 (MQ=255)
                                             gtagtaACAACCAAAAACAAATCAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCa                                                                                                                                                                 >  1:507635/1‑251 (MQ=255)
                                                                   cAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGtt                                                                                                                                            >  1:573292/1‑250 (MQ=255)
                                                                   cAATATCATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGtt                                                                                                                                            >  1:30624/1‑250 (MQ=255)
                                                                         caTTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTg                                                                                                                                      <  1:219089/249‑1 (MQ=255)
                                                                           tttGGCGTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAg                                                                                                                                   <  1:41621/251‑1 (MQ=255)
                                                                                 gTGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTaa                                                                                                                             >  1:594764/1‑251 (MQ=255)
                                                                                  tGATGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTa                                                                                                                                                       <  1:49029/224‑1 (MQ=255)
                                                                                    aTGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAg                                                                                                                          >  1:108104/1‑251 (MQ=255)
                                                                                    aTGCTACGGCAAAGAAAGCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAg                                                                                                                          <  1:111727/251‑1 (MQ=255)
                                                                                                 gaaagCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAAt                                                                                                                            <  1:207813/236‑1 (MQ=255)
                                                                                                  aaagCCGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTAtttttt                                                                                                              <  1:392321/249‑1 (MQ=255)
                                                                                                      ccGCAATGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAg                                                                                                                       <  1:679031/236‑1 (MQ=255)
                                                                                                          aaTGGAAATTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTAc                                                                                                    <  1:668755/251‑1 (MQ=255)
                                                                                                                aaTTCGTACAAAATTGATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATc                                                                                              <  1:299625/251‑1 (MQ=255)
                                                                                                                               gATTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAg                                                                                                                          <  1:405803/208‑1 (MQ=255)
                                                                                                                                aTTAAAAATGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACgg                                                                               <  1:349204/250‑1 (MQ=255)
                                                                                                                                       aTGGTATCGCCTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGt                                                                       >  1:117945/1‑251 (MQ=255)
                                                                                                                                                ccTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGcacgtac                                                               <  1:482277/250‑7 (MQ=255)
                                                                                                                                                ccTCTAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGcacgta                                                                >  1:810929/1‑244 (MQ=255)
                                                                                                                                                   ctAAATCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAg                                                                                                                          <  1:352969/188‑1 (MQ=255)
                                                                                                                                                     aaaTCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTAc                                                         <  1:862273/251‑1 (MQ=255)
                                                                                                                                                     aaaTCCATTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTAc                                                         <  1:244043/251‑1 (MQ=255)
                                                                                                                                                       atccaTTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACtg                                                       <  1:377486/247‑2 (MQ=255)
                                                                                                                                                        tccaTTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACtgt                                                      <  1:453893/248‑3 (MQ=255)
                                                                                                                                                         ccaTTTACTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACt                                                                                                   >  1:653752/3‑205 (MQ=255)
                                                                                                                                                               aCTTAGAGAAAAAAGCACATAAAGCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAg                                                                                                                                                         <  1:230080/145‑1 (MQ=255)
                                                                                                                                                                                      gCAAAAATATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATg                                                                                                           >  1:612077/1‑168 (MQ=255)
                                                                                                                                                                                         aaaaTATATCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGt                     <  1:28804/251‑1 (MQ=255)
                                                                                                                                                                                               tatCCAATCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAg                                                                                                                          >  1:524249/1‑144 (MQ=255)
                                                                                                                                                                                                 tCCAATCTATATTGAAGTAGCTCAAATTGCCCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACaaa             >  1:101749/1‑251 (MQ=255)
                                                                                                                                                                                                    aaTCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTg                                                                                                                                      >  1:507321/1‑127 (MQ=255)
                                                                                                                                                                                                    aaTCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACtgt                                                      <  1:484886/207‑3 (MQ=255)
                                                                                                                                                                                                     aTCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAg                                                                                                                          <  1:554554/138‑1 (MQ=255)
                                                                                                                                                                                                     aTCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACaaact           >  1:659905/1‑247 (MQ=255)
                                                                                                                                                                                                     aTCTATATTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACaaact           >  1:200446/1‑247 (MQ=255)
                                                                                                                                                                                                         atatTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATc                                                                                                                                                    >  1:793134/1‑108 (MQ=255)
                                                                                                                                                                                                         atatTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGtta                                                                                                        <  1:647247/152‑1 (MQ=255)
                                                                                                                                                                                                         atatTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGt     >  1:170678/1‑251 (MQ=255)
                                                                                                                                                                                                          tatTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTa    <  1:286019/251‑1 (MQ=255)
                                                                                                                                                                                                           atTGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTAc   <  1:690801/251‑1 (MQ=255)
                                                                                                                                                                                                            ttGAAGTAGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAGTATCAAGCCGTGATGATATTCAAGGTGAAGTTGCACAAACTCTACGTACt  <  1:649084/251‑1 (MQ=255)
                                                                                                                                                                                                                   aGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATc                                                                                                                                                    >  1:644313/1‑98 (MQ=255)
                                                                                                                                                                                                                   aGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAg                                                                                                                          <  1:544732/124‑1 (MQ=255)
                                                                                                                                                                                                                   aGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTAc                                                         <  1:233191/189‑1 (MQ=255)
                                                                                                                                                                                                                   aGCTCAAATTGCGCAGAGACCGGCGCCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCACGTACTATTACTGTAgtatc                                                >  1:246256/1‑195 (MQ=255)
                                                                                                                                                                                                                                 ||                                                                                                                                                                                                                                    
ATGTCAGTTATAGTAATTTATTTAATTCCGTTATTCGTGATATAGGTAGTAACAACCAAAAGCAAATCAATATAATTTGGCGTGATGCTACGGCAAAGAAAGCCGCAATGAAAATTCGTACAAAATTGATTAAAAATGGTATTACCTCTAAATCTATTTACTTAGAGAAAAAAGCACATAAAGCAAAAACATATCCAATCTATATTGAAGTAGCTCAAATTGCGCCCAAACCGGCACCTTGTGCATATCGTACCGGAGAAGTAGTTTTTTGGAATTCTAATGAGACTTCTTGTGCACTTAAAAGTAATCAGCGTGTTCAATTGAAGTATTAATAGAAGGTATTTTTTATGTTATTACTAGATCAAGAAACAATCACGGCAGACAGTGCTCGCACTATTACTATAGTGTCAAGCCGTGATGATATCCAAGGTGAAGTTGCACAAATTCTACGTACT  >  1/1481567‑1482021

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