Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP025268164,4120AC26.1% 208.3 / 15.3 161intergenic (+47/‑149)hrpB/mrcBATP‑dependent helicase HrpB/bifunctional glycosyl transferase/transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (55/64);  new base C (5/37);  total (60/101)
Fisher's exact test for biased strand distribution p-value = 6.92e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 162 total aligned reads displayed.

CAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCG  >  NZ_CP025268/164273‑164550
                                                                                                                                            |                                                                                                                                            
    aaGAGATTAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGTACTAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGATAGGGCCGGGGTGAGGGCATCAg                                                                                                                                          <  2:7748592/140‑1 (MQ=38)
        ttttaatgggAGTTATCCCAAACT‑TGTCTGGCCGGACGACCCGCAAAATACTGCACCGACGCGACGGACTAAAAATTTTTCGTAAGGTTGGTTTTCTCCCTCACCCTGTGGGAGTGGGCCGGGTTGCGGGCCTCAgcgcg                                                                                                                                      <  2:6215647/133‑1 (MQ=255)
        gaTTAATGGGCGTTATCCAAAACT‑TGTCTGGCCGGACGCCCCGCAAATTACTGCACCGACGCGACGTACTAAAAAGTATTCTAATGTTTGTTTTTCTCCCACTCCCTGTGGGTGAGGGCGGGGGTGAGGGCATCAgcgcg                                                                                                                                      <  2:9151232/140‑1 (MQ=19)
        gaTGACAGGGCTTTATCCCAAACT‑TTGGTGGCCGGACGCCCCGGCAAATCCTGCACCGACGCGAGGGACGAAAATTTATTCTTATGGTTGTTTTTCTCCCTCACCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAgcgcg                                                                                                                                      <  1:7986928/140‑1 (MQ=19)
        gaTGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAgcgcg                                                                                                                                      >  1:8689448/1‑140 (MQ=255)
          ttaacgggCGTTATCCCAAACT‑TGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCa                                                                                                                                    <  2:4427081/135‑1 (MQ=37)
                ggCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgga                                                                                                                              >  2:3458503/1‑137 (MQ=255)
                  cgatATCCAAAACT‑TGTCTGGCGTCACGCCCCCTCAATTACAGCACCGACGCGACTGTGTAAAAAGTATTCTTACTGTTGTTTTTATCCATCTCCCTTGGGGATAGGGCAGGGTTGCGGGCATCAGCGCGCACGTtc‑acc                                                                                                                           <  2:6484766/137‑1 (MQ=255)
                  cGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggacc                                                                                                                            >  1:5991954/1‑135 (MQ=255)
                    ttATCCCACCCA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGGGGGCATCAGCGCGCACGGtc‑accct                                                                                                                         <  2:8885572/140‑1 (MQ=38)
                    ttATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc‑accct                                                                                                                         <  2:1988100/140‑1 (MQ=255)
                    ttATCCCAAACA‑GGTCGGGCCGGACTCCCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCGTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc‑accct                                                                                                                         <  2:8681920/140‑1 (MQ=38)
                    ttATCCAAAACT‑GGTGTGGCCGGACGACCCGGAAAATCCTGCACCGACGCGAGGGAGTAAAATTTTTTCGTTTGGTGGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTtc‑accct                                                                                                                         <  2:5669151/140‑1 (MQ=19)
                    ttATCCAAAAAT‑TGTCTGGCCGTACGCCCCGGAAAATACTGCACCGACGCGAGGGACTAAAATTTTTTCGTATGGTTGGTTTTATCCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTtc‑accct                                                                                                                         <  2:4430914/140‑1 (MQ=255)
                    ttACCCCACACA‑TGTCTGCCCGGACGACCCGGCAAACCCTGCCCCGCCGCGACGGCCGAACACGTCTACGTACGGTTGGTTGTCGCCCCCCCCCTGTGGGCGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTtc‑accct                                                                                                                         <  2:3852172/140‑1 (MQ=255)
                    ttACCCCAAACA‑TGTCTGGCCGGACGACCCGTAAAAACCGGCACCGACGCGCCGGACGAAAAAGTTTTCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc‑accct                                                                                                                         <  2:9046909/140‑1 (MQ=25)
                    ttACCCCAAAAT‑TTTGTGGCCGGCCGACCGGGCAAATCCGCCACCGACGCTACGGCCTAAAATGTATTCTTATGTTTGTTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGTTGTGGGCATCAGCGCGCACGTtc‑accct                                                                                                                         <  2:6232146/140‑1 (MQ=19)
                     ttatccAAACATTGTCGGGCCGTACTCCCCGGAAAATACGTCACCGACGCTCCGTACTAAAATTTTTTTTTATGGTTGGTTTTCTCCCACTCCGTGGGGGAGTGGGCGGGGTTGTGTCCATCAGCGCGCACGTtc‑accct                                                                                                                         <  2:5738241/136‑1 (MQ=255)
                         ccAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgc‑acccgcaccc                                                                                                                    >  2:6831874/1‑140 (MQ=255)
                           aaaCA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgc‑ccccacacccga                                                                                                                  >  2:8216131/1‑138 (MQ=38)
                            aaCA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTaa                                                                                                                 <  2:9819184/140‑1 (MQ=255)
                            aaCA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCGCACCCCaa                                                                                                                 >  2:5854084/1‑140 (MQ=255)
                              cA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAAcc                                                                                                               <  2:6090748/140‑1 (MQ=255)
                                ttGTCTGGCGGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAATTTTTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAAccc                                                                                                              <  2:4306969/139‑1 (MQ=37)
                                  gtggggCCGCACGCCCCGCCAAACACGCCCCCGCCGCGCCGGACCAAAAATTATTCTTACGGTGGTTTTGCCCCCACCCCCTGTGGGTGTGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCctc                                                                                                            <  1:9519113/136‑1 (MQ=255)
                                  gTCTGGCGGGACGACCCGGCAAATCATGCACCGACGCGCCGGAGTAAACAGTATTCGTCAGGTTGGTTGTCCCCCACTCCCTGTGGGATTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCctc                                                                                                            <  2:7227038/140‑1 (MQ=37)
                                  gTCTGGCCGTACGACCCGGCAAATACTCCACCGACGCGACGGACTAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCctc                                                                                                            <  2:7262575/140‑1 (MQ=38)
                                  gTCTGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCctc                                                                                                            <  2:804061/140‑1 (MQ=255)
                                  gTCTGGCCGGACGACCCGGCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTC‑ACCCACCCCCCAAccccc                                                                                                            >  2:3744645/1‑138 (MQ=37)
                                  gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGC‑ACCCTCACCCGAAcccac                                                                                                            >  2:205779/1‑138 (MQ=38)
                                  gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGG‑GCCCACACCCAAAcccac                                                                                                            >  1:9255361/1‑138 (MQ=37)
                                  gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGG‑ACCCGCCCCCAAAaccac                                                                                                            >  1:8792183/1‑138 (MQ=37)
                                  gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGG‑ACCCCACCCCCAAcccac                                                                                                            >  1:3762079/1‑138 (MQ=37)
                                  gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGC‑ACCCTCACCCAAAcccac                                                                                                            >  2:2979968/1‑138 (MQ=255)
                                  gTCTGGACTGACGCCCCCTCAATTACTGACCCTCCGCGAAGGACGAATACGTTTTCGTAAGGTTGGTTTTCTCCATCTCCCTGTGGGAGTGGTCCGTGGTGAGGTCATCAGCGCGCACGGAC‑ACACTCACCCTTACCctc                                                                                                            <  2:7622711/140‑1 (MQ=255)
                                  gTCGGGCCGGACGACCCTGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACTGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCctc                                                                                                            <  2:957723/140‑1 (MQ=255)
                                    cgggccggccGACCCGGAAAATACTGCACCGCCGCGAGGGACCAAAAATTATTCGTATGGTTGTTTTTCTCCCTCTCCCGGGGGGAGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCtctc                                                                                                          <  2:3506463/138‑1 (MQ=255)
                                    cTGGCCGGACGACCCGGCAAATCCTGCACCGCGGCGCCGGACGAAAAAGAATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCtctc                                                                                                          <  2:8442219/140‑1 (MQ=255)
                                    cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCtctc                                                                                                          <  2:6432091/140‑1 (MQ=255)
                                    cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCGCACCCTAACCCtctc                                                                                                          >  1:9872567/1‑140 (MQ=255)
                                     tGGCCGGACGCCCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTcc                                                                                                         <  2:5025367/140‑1 (MQ=255)
                                       gCCGGACGACCCGGCAATTCCGGCACCGACGCGACGGACGAAAAAGTATTCGTCCGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTACCCCTCTCCCt                                                                                                       <  2:3982306/140‑1 (MQ=37)
                                       gCCGGACGACCCGGCAAATACTGCACCGACGCGAGGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCACCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCt                                                                                                       <  2:9706606/140‑1 (MQ=255)
                                       gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCACACCCAAACCCTCTCACt                                                                                                       >  2:184886/1‑140 (MQ=255)
                                        ccGGACGACCCGCAAAATACTCCACCGACGCGACGGACTAAAAATTTTACGTATGGTTGTTTTTCACCCTCCCCCTGTGGCAGTGGGCCGGGGTGAGGCCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTc                                                                                                      <  1:4068903/140‑1 (MQ=255)
                                          ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCGGCGCGCACGTGC‑ACCCTCAACCAAACCCGATcccgggt                                                                                                    >  1:2012625/1‑136 (MQ=37)
                                             cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTG‑ACCCTCACCCCAAACCTCGCCCCCCAggg                                                                                                 >  1:7671698/1‑140 (MQ=37)
                                              gacCCGGCAAATCCTCCACCGCCGCGCGGGACGAAAATTTTTGCTTATGGTTGTTTTCCCCCCTCACCCTGTGGTTGAGGGCGGGGGTGGGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                <  1:2069578/140‑1 (MQ=255)
                                                cccGGCAAATACTGCCCCGACGCGACGGACTAAAATGTATTCGTAAGGTTGGATTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                              <  2:8677505/140‑1 (MQ=255)
                                                cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                              <  2:7739127/140‑1 (MQ=255)
                                                cccGGAAACTCCTGCCCCGCCGCGAGGGAGTAAAATGTATTCGTATGTTTGGTTTTCTCCCTCACCCTGGGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                              <  2:4804321/140‑1 (MQ=255)
                                                      aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc‑accc                                                                                                                          >  1:7979726/1‑106 (MQ=255)
                                                          acggCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGATGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCCCCCTCAAGGGAGAGGTGACCGAt                                                                                    <  2:2719223/137‑1 (MQ=39)
                                                          aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTA‑ACCCACACCCTAACCCACACCCGCAAGGGTGGGGGGACCGAt                                                                                    >  1:8983968/1‑140 (MQ=37)
                                                          aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGC‑ACCCCCACCCAAAACCACACCCACAACGGGGAGGGGAcccct                                                                                    >  1:9826964/1‑137 (MQ=25)
                                                          aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGC‑ACCCACACCCTAAACCACTACCACCAGGGAGTGGGGACCGAt                                                                                    >  2:1860887/1‑140 (MQ=37)
                                                          aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGGGC‑ACCCTCACCCTAACCCTCTCCCGCCCGGGCGCGGGGAcccag                                                                                    >  2:8497389/1‑137 (MQ=37)
                                                           cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                 <  2:1018509/95‑1 (MQ=255)
                                                           cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                 >  1:1018509/1‑95 (MQ=255)
                                                            ggCACCGACGCTACGTACTAAAAAGTATTCTTAAGGTTGTTTTTCTCCCTCTCCGTGTGGGTTTGGGCCGGGGTGAGGTCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                  <  2:4717917/139‑1 (MQ=37)
                                                              cACCGACGCGACGGACTAAAATGTTTCCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                <  1:1684207/140‑1 (MQ=38)
                                                              cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTG‑GCCCCCACCCAAAACCACACCCACCAAGGAGAGGGGcccccacgac                                                                                >  1:559785/1‑133 (MQ=21)
                                                                 cgccgcgacgGACTAAAACTTATTCGTTAGTTTGGTTTTCTCCCTCTCCCGGGGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca‑c                                                                            <  2:2283293/137‑1 (MQ=37)
                                                                  gccccgacgTACGAAACCGTATTCCTAAGGTGGGTCAGCACCCTCCCCCTCTGGGCCCGGGCCGGGCTGCGGGCATCCCCGCCCACGTAC‑ACCCTCACCCTAACCCTCACCCTCAAGGGAGAGGGGACCCATCGAGca‑ca                                                                           <  2:8756106/136‑1 (MQ=255)
                                                                   acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCCAGGGAGAGGGGACCGATCGAGCA‑Caa                                                                          >  2:8894203/1‑140 (MQ=255)
                                                                   acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGTC‑CCCCGCACCCAAACCCTCTCCCTCAAGGGAGAGGGGTCCGCTCGAGCA‑Caa                                                                          >  2:9486995/1‑140 (MQ=37)
                                                                   acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTG‑CCCCTCACCCAAAACCACACCCTAAAGGGGGAGGGGACCGATCGAGAA‑Caa                                                                          >  1:3099143/1‑140 (MQ=25)
                                                                   acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCCCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑Caa                                                                          >  1:4027997/1‑140 (MQ=255)
                                                                    cgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGG‑GACCTCACCCAAACCCACAACCACAAAGGGGGGGGGGCCGCGCGGGca‑caca                                                                         >  1:1443662/1‑138 (MQ=18)
                                                                          ggactAAAATGTATTCTTAAGTTTGGTTTTCTCCCTCGCCCTGTGGTTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTga                                                                   <  1:1785952/140‑1 (MQ=38)
                                                                           gactAAAAAGTATCCGTAAGGTTGGTTTACTCCCTCGCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTgag                                                                  <  1:7590620/140‑1 (MQ=38)
                                                                             ctaaaaaGTTTTCTTATGGTTGGTTTTCTCCCTCTCCGTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTgagag                                                                <  2:5953666/138‑1 (MQ=39)
                                                                             cgAAAAATTTTTCGTATGTTTGGTTTTCTCCCTCTCCGTGTGGTTGTGGGCGGGGTTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTgagag                                                                <  2:8150590/140‑1 (MQ=38)
                                                                                aactgTATGCGTACGGTTGGTTATCGCCCTCCCCCGGTGGCTGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGata                                                             <  1:9293788/136‑1 (MQ=37)
                                                                                        tggTATGGTTGTTTTTCACCCTCGCCCTGTGGGTGTGGCCGGGGGTGAGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATAtcttcttc                                                     <  1:8236271/138‑1 (MQ=38)
                                                                                        tctttaGGTTGGTTTACTCCCTCTCCCTGGGGGTGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAANCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATAtcttcttc                                                     <  1:6319021/135‑1 (MQ=255)
                                                                                        tCGTACGGTTGGTTTTCTCCCTGGCCCGGGGGGTGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAACTTTTGAGAGATAtcttcttc                                                     <  2:4255496/140‑1 (MQ=38)
                                                                                        tCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATAtcttcttc                                                     <  2:2689078/140‑1 (MQ=255)
                                                                                        ccgaacggttggttTGCCCCCACGCCCGGGGGGAGAGGCCCGGGGTGAGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATAtcttcttc                                                     <  1:6779534/134‑1 (MQ=255)
                                                                                          gTATGGTGGGTTTACTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTg                                                   <  1:5097008/140‑1 (MQ=255)
                                                                                          gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTAACCCGAACCCTCTCCCTCCAGGGTGAGGGGACCGATCGAGCA‑CAAAATTTGAGAGATATCTTATTCTg                                                   >  2:7643672/1‑140 (MQ=37)
                                                                                           gacggttgtttTTCCCCCCCCCCGGGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGt                                                  <  1:1620734/137‑1 (MQ=255)
                                                                                             tggtgggttTTCGCCCACTCCCGGGGGGAGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCt                                                <  1:8468441/139‑1 (MQ=38)
                                                                                               gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCACAAGGGAGAGGGGACCGATCGAGCC‑CAAAATTTGAGAGATATATTCTTATGTCTTg                                              >  2:7753365/1‑140 (MQ=38)
                                                                                                tgggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTCAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                             <  2:6690392/138‑1 (MQ=255)
                                                                                                tgggttTTCGCCCTCGCCGTGTGGGAGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                             <  1:4932715/138‑1 (MQ=39)
                                                                                                      ttCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                       <  2:1443662/140‑1 (MQ=39)
                                                                                                       gctccctctccctGTGGGAGAGGCCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                      <  1:8785780/139‑1 (MQ=255)
                                                                                                          ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGG‑GCCCCCACCCCAACCCCCTCCCACACGGGGGAGTGGGGCGAGCGGGCACCAAATTTTG‑GAGGTATCTTATTAAGTCTTGTGACAgaagaa                                   >  2:2325907/1‑140 (MQ=12)
                                                                                                                          gagGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                   <  1:7040902/140‑1 (MQ=255)
                                                                                                                           agGGCCGGGGTGAGGGCATCAGCGCGCACGGTC‑ACCCGCACCCAAACCCGCACCCTCAAGGGAGAGGGGGCCGGTAGAGCA‑CAAAATTTGTGAGATATATTATTATTTGTTTTAAAACAAGAACAAAAAAAAGGGCtttt                  >  2:3095277/1‑140 (MQ=255)
                                                                                                                            ggcccGGGGGGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                 <  1:9724615/137‑1 (MQ=255)
                                                                                                                             ggCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAAACCTCACCCTCAAGGGTGAGGGGGCCCATCGAGCA‑CAAAATTTGTTTGATATATTATTATGTGTTTTAACAGAAGAACAGAAAATCgggggttttg                >  2:8442763/1‑133 (MQ=255)
                                                                                                                             ggCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCACCCCCTAAACCAAACCCTCAAGGGTGAGGGGTCCCATCGCGCA‑CAAATTTTGAGGTTTATAATAATATTTCGTGTTAAAGAAGAACAAAAAATAGGGCtttttg                >  2:7016745/1‑138 (MQ=255)
                                                                                                                             ggCCGGGGTGAGGGCATCAGCGCGCACGTGCAACCCCCACCGAAACCC‑CACCCACAACGGAGAGGGGCCCGATCCGCCC‑CAAAAATTGATAGATATCTTAATATTGCTTGTAACACAATAAAAAAAAAAAGGGCtttttg                >  2:1291621/1‑138 (MQ=255)
                                                                                                                               ccGGGGTGAGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc              <  1:6870629/140‑1 (MQ=255)
                                                                                                                               ccGGGGTGAGGGCCTCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc              <  2:5155686/140‑1 (MQ=255)
                                                                                                                                 ggggTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt            <  1:414571/140‑1 (MQ=255)
                                                                                                                                   ggTGAGGGCATCAGCGCGCACGTTC‑ACCCACACCCGAAACCACACCCACAAGGGAGAGGGGCCCGCTCGAGCA‑CAAATTTTGTGAGATATAATATTATGTCTTGTAACAGAAGAACAGAAAATAGGGCTGTTGCGCCGGa          >  2:6564818/1‑140 (MQ=255)
                                                                                                                                           cATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg  <  1:2979968/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                            
CAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTC‑ACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCA‑CAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCG  >  NZ_CP025268/164273‑164550

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