Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,479,408 T→C Y73Y (TAT→TAC)  unknown → Flp pilus assembly protein RcpC/CpaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,479,4080TC100.0% 171.9 / NA 60Y73Y (TAT→TAC) unknownFlp pilus assembly protein RcpC/CpaB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (40/20);  total (40/20)

AGGGTAAAAAATGAATTACAGAACGCTCTTGATCATTTCAGTGCTTACGCTAGCGGTTGGCATTGGTGGCATTATGTTTATGCCAACATTTGACAATACTGAAAATCAACAAACAATAATAAATAATCAGCAGGTTCCGGAAGCTGAGAAAAAACCGGAAAAGCTAATTTTAACTGCCGAGTTAAAGAATGACGTAAGCAAAGGTCATTTATTACAGGCAGATGATTATATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTAGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAACGATCCTAGATTCTTTACATCCAGCCTAGATCCGAATCAAGAAGTTGCATATCGCATCTATGTTAAATCTGGTGATGATTACATT  >  1/1479180‑1479654
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
gtaaaaaatgaattacagaacgctcttgatcatttcagtgcttacgctagcggttggcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTc                                                                                                                                                                                                                                   >  1:623529‑M2/115‑250 (MQ=255)
      aatgaattacagaacgctcttgatcatttcagtgcttacgctagcggttggcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTGCAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAg                                                                                                                                                                                                                            <  1:232269‑M2/143‑1 (MQ=255)
       atgaattacagaacgctcttgatcatttcagtgcttacgctagcggttggcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGa                                                                                                                                                                                                                           >  1:497831‑M2/108‑251 (MQ=255)
        tgaattacagaacgctcttgatcatttcagtgcttacgctagcggttggcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAGAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGAc                                                                                                                                                                                                                          <  1:212368‑M2/145‑1 (MQ=255)
                 gaacgctcttgatcatttcagtgcttacgctagcggttggcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCa                                                                                                                                                                                                                    >  1:399151‑M2/98‑248 (MQ=255)
                         ttgatcatttcagtgcttacgctagcggttggcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACaa                                                                                                                                                                                                         <  1:774508‑M2/162‑1 (MQ=255)
                                tttcagtgcttacgctagcggttggcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTa                                                                                                                                                                                                  >  1:10872‑M2/83‑251 (MQ=255)
                                            cgctagcggttggcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACg                                                                                                                                                                                      >  1:178130‑M2/71‑251 (MQ=255)
                                                   ggttggcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGc                                                                                                                                                                                          >  1:674694‑M2/64‑240 (MQ=255)
                                                   ggttggcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTaaaa                                                                                                                                                                               >  1:235054‑M2/64‑251 (MQ=255)
                                                    gttggcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAagag                                                                                                                                                                                              <  1:237051‑M2/173‑1 (MQ=255)
                                                        gcattggtggcattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGc                                                                                                                                                                                                                     >  1:471454‑M2/59‑208 (MQ=255)
                                                                  cattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCAc                                                                                                                                                                                                                   >  1:562842‑M2/49‑200 (MQ=255)
                                                                  cattatgtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGAt                                                                                                                                                                >  1:481173‑M2/49‑251 (MQ=255)
                                                                        gtttatgccagcatttgacaatactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAg                                                                                                                                                          >  1:150762‑M2/43‑251 (MQ=255)
                                                                                            atactgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAAt                                                                                                                                      >  1:873148‑M2/23‑251 (MQ=255)
                                                                                               ctgaaaatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAgg                                                                                                                                                                   >  1:37633‑M2/20‑219 (MQ=255)
                                                                                                    aatcaacaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGc                                                                                                                              <  1:376894‑M2/237‑1 (MQ=255)
                                                                                                         acaaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACCATGACCTAAGAGAGGCTACGAc                                                                                                                                                                                    >  1:433360‑M2/10‑192 (MQ=255)
                                                                                                          caaacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTc                                                                                                                        <  1:828673‑M2/243‑1 (MQ=255)
                                                                                                            aacaatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGAc                                                                                                                                                                                    <  1:502241‑M2/183‑1 (MQ=255)
                                                                                                              caataataatAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTg                                                                                                                                                                                                                               >  1:41512/5‑144 (MQ=255)
                                                                                                              caatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTAc                                                                                                                                                             >  1:504269‑M2/5‑210 (MQ=255)
                                                                                                              caatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGc                                                                                                                                                        >  1:641245‑M2/5‑215 (MQ=255)
                                                                                                              caatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCtttt                                                                                                                    >  1:732558‑M2/5‑251 (MQ=255)
                                                                                                              caatAATAATAAATAATCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCtttt                                                                                                                    >  1:133137‑M2/5‑251 (MQ=255)
                                                                                                                       taaataaTCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAagag                                                                                                                                                                                              <  1:286945/168‑1 (MQ=255)
                                                                                                                       taaataaTCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAagag                                                                                                                                                                                              >  1:159570/1‑168 (MQ=255)
                                                                                                                       taaataaTCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTaa                                                                                                                                                 >  1:106041/1‑213 (MQ=255)
                                                                                                                       taaataaTCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGc                                                                                                                              >  1:242615/1‑232 (MQ=255)
                                                                                                                       taaataaTCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCcaa                                                                                                            >  1:90435/1‑250 (MQ=255)
                                                                                                                       taaataaTCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCcaa                                                                                                            >  1:100097/1‑250 (MQ=255)
                                                                                                                           taatCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAagag                                                                                                                                                                                              <  1:232937/164‑1 (MQ=255)
                                                                                                                           taatCAGTCGGTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAAt                                                                                                       >  1:485190/1‑251 (MQ=255)
                                                                                                                            aatcagtcggTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAAt                                                                                                       >  1:11954/9‑250 (MQ=255)
                                                                                                                                   cggTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTAtc                                                                                               <  1:665572/250‑1 (MQ=255)
                                                                                                                                   cggTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTAtc                                                                                               <  1:279032/250‑1 (MQ=255)
                                                                                                                                     gTTCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTAtc                                                                                               <  1:105153/249‑1 (MQ=255)
                                                                                                                                       tCCGGAAACTGAGAAAAAGCCGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGc                                                                                                                              <  1:347805/216‑1 (MQ=255)
                                                                                                                                                      aaagccGGAAAAGCTAATTTTAACTGCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGAt                                                                            <  1:814244/247‑1 (MQ=255)
                                                                                                                                                                               gCCGAATTAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGGGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGcttgaatccaag                                                   >  1:814100/1‑240 (MQ=255)
                                                                                                                                                                                  gaattAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAAAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGcttgaatccaagtca                                                >  1:426618/4‑237 (MQ=255)
                                                                                                                                                                                    attAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTAtct                                                                                              >  1:563279/2‑203 (MQ=255)
                                                                                                                                                                                    attAAAGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTActtc                                                                  <  1:324211/230‑4 (MQ=255)
                                                                                                                                                                                       aaaGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGc                                                                                                                                                        >  1:407299/1‑142 (MQ=255)
                                                                                                                                                                                       aaaGAATGACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGAt                                                                            <  1:708772/218‑1 (MQ=255)
                                                                                                                                                                                              gACGTGAGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGATATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGAt                                                                            >  1:486708/1‑211 (MQ=255)
                                                                                                                                                                                                cgtgaGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCAc                                                                                                                                                                                                                   <  1:419966/70‑1 (MQ=255)
                                                                                                                                                                                                   gaGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAg                                                                                                                                           <  1:136307/142‑1 (MQ=255)
                                                                                                                                                                                                   gaGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCa                               >  1:681385/2‑251 (MQ=255)
                                                                                                                                                                                                   gaGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCa                               >  1:126383/2‑251 (MQ=255)
                                                                                                                                                                                                    aGCAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCAt                              >  1:711836/1‑251 (MQ=255)
                                                                                                                                                                                                      cAAAGGTCATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATc                             <  1:296670/250‑1 (MQ=255)
                                                                                                                                                                                                             cATTTATTACAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTc                                                                                                                   >  1:578299/1‑157 (MQ=255)
                                                                                                                                                                                                                  attaCAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATcc                                                                                                              >  1:865915/1‑157 (MQ=255)
                                                                                                                                                                                                                  attaCAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTagat                   >  1:761314/1‑245 (MQ=255)
                                                                                                                                                                                                                     aCAGGCAGATGATTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTaaaa                                                                                                                                                                               <  1:790021/89‑1 (MQ=255)
                                                                                                                                                                                                                              tgatTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTActtc                                                                  >  1:110222/1‑186 (MQ=255)
                                                                                                                                                                                                                              tgatTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACa    >  1:736366/1‑251 (MQ=255)
                                                                                                                                                                                                                                atTACATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTGGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACAtt  >  1:348346/1‑251 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
AGGGTAAAAAATGAATTACAGAACGCTCTTGATCATTTCAGTGCTTACGCTAGCGGTTGGCATTGGTGGCATTATGTTTATGCCAACATTTGACAATACTGAAAATCAACAAACAATAATAAATAATCAGCAGGTTCCGGAAGCTGAGAAAAAACCGGAAAAGCTAATTTTAACTGCCGAGTTAAAGAATGACGTAAGCAAAGGTCATTTATTACAGGCAGATGATTATATTTTAAATGAGATTTCTGTCCCTGAAGACAGTGCACTAGTCGACAATGACCTAAGAGAGGCTACGACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAACGATCCTAGATTCTTTACATCCAGCCTAGATCCGAATCAAGAAGTTGCATATCGCATCTATGTTAAATCTGGTGATGATTACATT  >  1/1479180‑1479654

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