Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 645,748 T→C P488P (CCT→CCC)  unknown → ATP‑dependent helicase HrpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1645,7480TC98.7% 204.4 / NA 77P488P (CCT→CCC) unknownATP‑dependent helicase HrpA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (28/48);  total (28/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGCTATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGGCATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCTTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAATGCTACTAAGTGCGGTCAATTTTTGCAGCGTTTATGAAGT  >  1/645510‑645990
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
gCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGt                                                                                                                                                                                                                                        >  1:702270/1‑251 (MQ=255)
      gTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTc                                                                                                                                                                                                                                   >  1:344591/1‑250 (MQ=255)
             tCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTccccc                                                                                                                                                                                                                               <  1:339053/247‑1 (MQ=255)
                  aTTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAAc                                                                                                                                                                                                                        <  1:563789/249‑1 (MQ=255)
                   ttCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAgg                                                                                                                                                                                                                     <  1:769064/251‑1 (MQ=255)
                     cGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAgggg                                                                                                                                                                                                                   <  1:742171/251‑1 (MQ=255)
                       gAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGCAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAgggggg                                                                                                                                                                                                                 >  1:454942/1‑250 (MQ=255)
                                  caacaaTCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTa                                                                                                                                                                                                      <  1:231214/251‑1 (MQ=255)
                                         cGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATc                                                                                                                                                                                               <  1:181570/251‑1 (MQ=255)
                                          gCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCa                                                                                                                                                                                              >  1:541026/1‑251 (MQ=255)
                                            ccGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCaga                                                                                                                                                                                            <  1:637738/251‑1 (MQ=255)
                                               gAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAAc                                                                                                                                                                                                                        <  1:527926/220‑1 (MQ=255)
                                               gAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAgggg                                                                                                                                                                                                                   <  1:738921/225‑1 (MQ=255)
                                               gAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATc                                                                                                                                                                                               <  1:222402/245‑1 (MQ=255)
                                                   ttACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGa                                                                                                                                                                                     <  1:495214/251‑1 (MQ=255)
                                                   ttACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGa                                                                                                                                                                                     <  1:525827/251‑1 (MQ=255)
                                                   ttACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGa                                                                                                                                                                                     <  1:888539/251‑1 (MQ=255)
                                                            cGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGGAGCCTTCCCATTCGTGGATGCACCGGATAAACGGCATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGGAAATCT‑CCCCCCTCCGCCTTCGGTGCTTCCACCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGctc                                                                                                                                                                            <  1:411094/251‑1 (MQ=255)
                                                                          aCCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAgg                                                                                                                                                              <  1:218417/251‑1 (MQ=255)
                                                                                ttGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATc                                                                                                                                                        <  1:384957/251‑1 (MQ=255)
                                                                                   gCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGTGGGATCTg                                                                                                                                                      >  1:42004/1‑250 (MQ=255)
                                                                                     gTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGcc                                                                                                                                                    >  1:147602/1‑250 (MQ=255)
                                                                                            aTTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGc                                                                                                                                            <  1:461430/251‑1 (MQ=255)
                                                                                            aTTTTGCAAATGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGc                                                                                                                                            <  1:117252/251‑1 (MQ=255)
                                                                                                   aaaTGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAgggg                                                                                                                                     <  1:606023/251‑1 (MQ=255)
                                                                                                    aaTGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAggggg                                                                                                                                    <  1:276967/251‑1 (MQ=255)
                                                                                                    aaTGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAggggg                                                                                                                                    >  1:458839/1‑251 (MQ=255)
                                                                                                      tGACGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAggggggg                                                                                                                                  <  1:783634/251‑1 (MQ=255)
                                                                                                         cGGTGTTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAAt                                                                                                                               <  1:233624/251‑1 (MQ=255)
                                                                                                           gtgtTGGGTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAggggg                                                                                                                                                                                                                  >  1:848044/1‑166 (MQ=255)
                                                                                                                  gTTTGGATGATATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAAc                                                                                                                                                                                                                        >  1:234460/1‑153 (MQ=255)
                                                                                                                           atATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAg                                                                                                                                                               >  1:32710/2‑201 (MQ=255)
                                                                                                                           atATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGaa                                                                                                                                >  1:7202/2‑232 (MQ=255)
                                                                                                                           atATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGaa                                                                                                             >  1:542910/2‑251 (MQ=255)
                                                                                                                             aTTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGaaaa                                                                                                           >  1:693197/1‑251 (MQ=255)
                                                                                                                                   gCCTTCCCGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATAc                                                                                                                             >  1:519697/1‑227 (MQ=255)
                                                                                                                                         ccGTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTcac                                                                                               <  1:341354/251‑1 (MQ=255)
                                                                                                                                           gTTCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACa                                                                                                                                                                >  1:498687/1‑184 (MQ=255)
                                                                                                                                            ttcggtgATGGACCGGATAAACGGCATATTCCGGGTGGGGTGAAATTGTTGGAAGGAGTGGGGGCGTTTGAATGGAAGGCCTA‑TCGCAAATCTGCCCCCTT‑CGCCTTCCGTTAGTCCCCCCCTTAACCGGGGGAGGGAAAGGGTATTTTATCAGGAGGAGGAACTTCGCTCCCCTCGTTTTCAAGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTcacac                                                                                             <  1:444494/244‑1 (MQ=255)
                                                                                                                                            ttCGTGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTtc                                                                                                   >  1:637673/1‑244 (MQ=255)
                                                                                                                                                tGGATGCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGaa                                                                                                                                                                                                            <  1:39958/135‑1 (MQ=255)
                                                                                                                                                     gCACCGGATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTAttt                                                                                                                                                                                                   >  1:223578/1‑139 (MQ=255)
                                                                                                                                                          ggATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAGTGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAAGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCaa                                                                               >  1:445562/1‑250 (MQ=255)
                                                                                                                                                          ggATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGGAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAgg                                                                                                                                       <  1:607284/194‑1 (MQ=255)
                                                                                                                                                          ggATAAACGACATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAgg                                                                                                                                       <  1:340012/194‑1 (MQ=255)
                                                                                                                                                                     atatTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTc                                                                   <  1:601468/251‑1 (MQ=255)
                                                                                                                                                                       atTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAAc                                                                                                                                                                                                                        <  1:180868/100‑1 (MQ=255)
                                                                                                                                                                       atTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGctc                                                                                                                                                                            >  1:745838/1‑144 (MQ=255)
                                                                                                                                                                       atTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAgg                                                                                                                                              <  1:421728/174‑1 (MQ=255)
                                                                                                                                                                       atTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTcc                                                                                                                          <  1:162226/194‑1 (MQ=255)
                                                                                                                                                                       atTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACt                                                                                                                      <  1:496075/198‑1 (MQ=255)
                                                                                                                                                                       atTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCt                                                                                         <  1:509144/227‑1 (MQ=255)
                                                                                                                                                                       atTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTcc                                                                  >  1:264396/1‑250 (MQ=255)
                                                                                                                                                                             gATGGCGTGTAATTTTTGGGAAAATTTGGGGCGTTGGAATTGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGt                                                              <  1:513589/248‑1 (MQ=255)
                                                                                                                                                                             gATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGt                                                              <  1:721227/248‑1 (MQ=255)
                                                                                                                                                                             gATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGt                                                              <  1:737963/248‑1 (MQ=255)
                                                                                                                                                                               tGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATc                                                         <  1:897138/251‑1 (MQ=255)
                                                                                                                                                                                 gCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATc                                                                                                                                                                                               <  1:595172/115‑1 (MQ=255)
                                                                                                                                                                                    tGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGc                                                    <  1:489510/251‑1 (MQ=255)
                                                                                                                                                                                      aaaTTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCt                                                   <  1:656897/250‑1 (MQ=255)
                                                                                                                                                                                      aaaTTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCt                                                   <  1:274565/250‑1 (MQ=255)
                                                                                                                                                                                      aaaTTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCt                                                   <  1:173768/250‑1 (MQ=255)
                                                                                                                                                                                               gaagaaTTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAAt                                         >  1:188644/1‑251 (MQ=255)
                                                                                                                                                                                                  gaaTTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAgg                                                                                                                                       <  1:211685/154‑1 (MQ=255)
                                                                                                                                                                                                              ttGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCt                                                   <  1:800827/226‑1 (MQ=255)
                                                                                                                                                                                                              ttGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAATGCTACTAAGTGCGGt                          >  1:161646/1‑251 (MQ=255)
                                                                                                                                                                                                              ttGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAATGCTACTAAGTGCGGt                          >  1:838312/1‑251 (MQ=255)
                                                                                                                                                                                                                     ggaaGG‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAATGCTACTAAGTGCGGTCAAtttt                   >  1:531446/1‑251 (MQ=255)
                                                                                                                                                                                                                          g‑CTATTCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAATGCTACTAAGTGCGGTCAATTTTTGc                <  1:434640/249‑1 (MQ=255)
                                                                                                                                                                                                                                tCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCacacac                                                                                           >  1:30648/1‑169 (MQ=255)
                                                                                                                                                                                                                                tCGCAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAATGCTACTAAGTGCGGTCAATTTTTGCAGCGttt         >  1:100586/1‑251 (MQ=255)
                                                                                                                                                                                                                                   cAAATCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAATGCTACTAAGTGCGGTCAATTTTTGCAGCGTTTATg      <  1:687521/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aaaTCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAATGCTACTAAGTGCGGTCAATTTTTGCAGCGTTTATGa     <  1:866387/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aaTCT‑CCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAATGCTACTAAGTGCGGTCAATTTTTGCAGCGTTTATGaa    <  1:46619/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tct‑cCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGaa                                                                                                                                >  1:637848/1‑125 (MQ=255)
                                                                                                                                                                                                                                       tct‑cCCCCCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAATGCTACTAAGTGCGGTCAATTTTTGCAGCGTTTATGAAGt  <  1:114523/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ccccCTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCt                                                   >  1:394710/1‑198 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  
GCATTTGTATTCGTCTTTATTCGGAAGACGATTTCAACAATCGCCCGGAATTTACCGATCCGGAAATCCTGCGAACCAATTTGGCGTCGGTGATTTTGCAAATGACGGTGTTGGGTTTGGATGCTATTGAAGCCTTCCCGTTCGTGGATGCACCGGATAAACGGCATATTCAGGATGGCGTGAAATTGTTGGAAGAATTGGGGGCGTTGGAATGGAAGG‑CTATTCGCAAATCT‑CCCCCTTCCGCCTTCGGTACTTCCCCCCGTAAACAGGGGGAGGGAAAGGGTATTTTATCAGAAGAAGGATCTTCGCTCCCTCCGTTTACAGGGGGATCTGCCCGTAGGGCTGAGGGGGGGAATACTCCAACTTTCAAAGAAAAACGTGTTCTCACACACCTCGGTCGCCAACTCGCCCAGCTCCCAGTCGATCCATGCTTGGCAAAAATGCTACTAAGTGCGGTCAATTTTTGCAGCGTTTATGAAGT  >  1/645510‑645990

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