Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP025268164,4120AC27.9% 230.0 / 32.9 180intergenic (+47/‑149)hrpB/mrcBATP‑dependent helicase HrpB/bifunctional glycosyl transferase/transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (51/78);  new base C (10/40);  total (62/118)
Fisher's exact test for biased strand distribution p-value = 1.42e-02
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 185 total aligned reads displayed.

CAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGG  >  NZ_CP025268/164273‑164551
                                                                                                                                              |                                                                                                                                           
 aaaaaGAGATGAAAGGCCGTTATCCCAAACA‑TGACTGGCCGGACGACCCGGCAAATA‑CTGCACCGCCGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCACCCTGTGGG‑AGAGGGCCGGGGTGAGGGCAt                                                                                                                                            <  2:10654669/140‑1 (MQ=39)
       agaTGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                      >  2:6994369/1‑140 (MQ=255)
        gaTGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                     >  1:8993373/1‑140 (MQ=255)
        gaTGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                     <  2:11447397/140‑1 (MQ=255)
        gaTGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                     >  2:510795/1‑140 (MQ=255)
          tGAAAGGGCGTTATCCCACACA‑TGTCTGGCCGGACGACCCGGAAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGT‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                   <  2:8677417/140‑1 (MQ=255)
          tGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                                   >  1:4635763/1‑140 (MQ=255)
          tGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                   >  2:3228520/1‑140 (MQ=255)
          tGAAAGGGCGTTAACCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                   >  1:9991502/1‑140 (MQ=255)
              tgggCGTTATCCCAAACT‑TGTCGGGCCGGACGACCCGGCAAATC‑CTGCACCGACGCGCCGGACTAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGTGGGCATCAGCGCGCACGtt                                                                                                                               <  2:9339960/139‑1 (MQ=37)
              cgggCGTTATCCCAAACA‑TTTCGGGCCGGACGACCCGCCAAATA‑CTGCACCGCGGCGACGGACTAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCGGGTGGG‑AGTGGGCGGGGGTTGGGGCCTCAGCGCGCACGtt                                                                                                                               <  2:1262534/139‑1 (MQ=255)
              aGGGCGTTCTCCCCACCA‑TGTCTGGCCGGCCGACCCGGCAAAGA‑CTGCACCGACGCGACGGACGAAACCGAATTCGTATGGTTGGTTTTCTCCCGCTCCCTGTGGG‑AGAGGGCCGGGGGGAGGGCATCAGCGCGCACGtt                                                                                                                               <  2:6455485/140‑1 (MQ=25)
              aGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGtg                                                                                                                               >  1:5544432/1‑139 (MQ=255)
                 gCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAACAGTAGTCGTAAGGGTGGTTTTCTCCCTCTCCCTGTGGG‑TGAGGGCCGGGGTGAGGGCATCAGCGCGCCCgggcag                                                                                                                            >  2:9421309/1‑135 (MQ=38)
                  cGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAATTA‑CGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑TGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                           <  2:295629/140‑1 (MQ=255)
                  cGTTATCCCAAACA‑TGTCGGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGTTTTTCTCCCTCACCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                           <  2:1995969/140‑1 (MQ=255)
                  cGTCATCCCAAACA‑TGTCTGGCCGGACGACCCGCCAAATA‑CTGCACCGACGCTACGGACGAAAAAGTATACGTAAGGTTGGTTTTCTCCATCTCCATGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                           <  2:6496176/140‑1 (MQ=38)
                    ttttCCAAAACA‑TGTCTGGCCGTACGACCGGGCAAATA‑CTGCACCGACGCTAGGTAGTAAAATTTTTTGTTATGGTTGTTTTTCTCCCTCTCCCTGTGGT‑AGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                         <  2:8965827/137‑1 (MQ=19)
                    ttATCCCACACA‑TGTCTGGCCGGACGCCCCGGCAAAGA‑CTGCACCGACGCGCCGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcaccct                                                                                                                         <  2:10561028/140‑1 (MQ=37)
                    ttATCCCAAACA‑TTTCTGGCCGGACGACCCGTAAAATC‑CTGCACCTACGCGACGGACTAAAAAGTTTTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                         <  2:5805553/140‑1 (MQ=37)
                    ttATCCCAAACA‑TGTCTGGCCGGACGCCCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAATTTATTCGTACGGTTGGTTTTCGCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                         <  2:1652171/140‑1 (MQ=38)
                    ttATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                         <  2:8739504/140‑1 (MQ=255)
                    ttATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcaccct                                                                                                                         <  2:9529313/140‑1 (MQ=255)
                     taacccAAACATTGTTTGGCCGGACGCCCCGTAAAATA‑CTGCACCGACGCGCCGTACTAAAATGTTTACTTAAGTTTGTTTTTCTCCCACTCCCTGTGGT‑AGTGGGCCGGGGTGAGGCCATCAGCGCGCACGTtcaccct                                                                                                                         <  2:3442456/137‑1 (MQ=255)
                      aTCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                       <  2:10346489/140‑1 (MQ=255)
                         ccACACA‑GGTCTGGCCGGACGCCCCGGCAAATC‑CGGCACCCACGCGACGGACGAAAAAGTCTTCGCAAGGTTGGTTTTCTCCCTCGCCCTGTGGG‑TGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcaccctcaccc                                                                                                                    <  2:10901996/140‑1 (MQ=255)
                         ccAAAAA‑TTTCTGGCCGTACGACCGGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCGGGGGTTAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                    <  2:1544560/140‑1 (MQ=38)
                          cAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTgggcccccacccc                                                                                                                   >  1:5847238/1‑139 (MQ=37)
                          cAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTcaccctcaccct                                                                                                                   >  2:10794839/1‑140 (MQ=255)
                              cA‑TGTGTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                               <  2:4038173/140‑1 (MQ=255)
                               a‑TGTCTGGCCGTACGACCCGGAAAATA‑CTGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑TGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                              <  2:1534945/140‑1 (MQ=255)
                                ttGTCGGGCGGGACGACCCGGAAACCA‑GTGCCCGGACGCGAGGTCCGAAAATGTTTGCTTACGGTTGGTTTACTCCCTCTCCCTGTGGG‑AGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                              <  1:3031067/139‑1 (MQ=255)
                                  gacTGGCCGTACCCCCCGCCAAATA‑CTGCACCTACGCTACGTACTAAAAATTATTCTAACGGTTGTTTTTCTCCCTCCCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                            <  1:2880229/138‑1 (MQ=21)
                                  gTCTGGCCGGCCGACCCGCCAAATA‑CTGCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGTTTTTCTCCCTCTCCCTGGGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                            <  1:10536690/140‑1 (MQ=38)
                                  gTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTTCACCCACACCCAAAcccac                                                                                                            >  2:5581195/1‑138 (MQ=255)
                                  gTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                            <  2:10706072/140‑1 (MQ=255)
                                  gTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAAcccac                                                                                                            >  2:9247136/1‑138 (MQ=255)
                                  gTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAAaccac                                                                                                            >  1:10490356/1‑138 (MQ=255)
                                  gTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAcccaa                                                                                                            >  1:10114715/1‑138 (MQ=255)
                                  gTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCGCCCGCACCCGAAcccac                                                                                                            >  1:7423684/1‑138 (MQ=38)
                                  gTCGGGCCGGACGCCCCGCCACATA‑CTGCACCTACGCGACGGCCGAAAAATTAGTCGAATGTTTGGTTTGCCCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                            <  1:7688529/140‑1 (MQ=25)
                                   tCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCcact                                                                                                           >  2:10344839/1‑140 (MQ=255)
                                   tCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCTCCCACACCCGAAACcact                                                                                                           >  1:3918838/1‑140 (MQ=37)
                                    gtGGCCGGACTCCGCGGCATATA‑GTGCACCGAGGCGACGTACTAAAATGTATTCGTAAGGTTGGTTTTCTCCCGATCCCTGTGGG‑AGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                          <  2:8481113/139‑1 (MQ=25)
                                    cgggCCGGACGACCCGCCACCTC‑CTGCCCCGACGCGCCGGACGAAAATTTATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGG‑AGTGGGCCGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                          <  1:1840354/138‑1 (MQ=25)
                                    cTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTTCACCCACACCCTAACCctcac                                                                                                          >  2:3732201/1‑138 (MQ=255)
                                    cTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCACACCCCAACCCactc                                                                                                          >  2:3732790/1‑140 (MQ=255)
                                    cTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCAAACCCtctc                                                                                                          >  2:11034152/1‑140 (MQ=255)
                                       gCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCGTGCACCCGCACCCGAACCCACAccca                                                                                                       >  2:2762805/1‑139 (MQ=37)
                                       gCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGGGCCCGCACCCAAACCCACAccca                                                                                                       >  1:8916082/1‑139 (MQ=37)
                                         cGGACGACCCCGCAACCA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCATCTCCCTGTGGT‑ATAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                     <  2:6557554/140‑1 (MQ=38)
                                           gacgccCCGGCAAATA‑CGGCACCGCCGCGCCGGACGAAAAAGTATTCGAACGGTTGGTTTTCTCCCTCCCCCTGTGGG‑AGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                   <  2:3021099/140‑1 (MQ=38)
                                             cctcccGGCAAAAACCTGCACCCAGGCCCCGCCCAAAAATTTTTTCGTATGTTTGTTTTACCCCCTCCCCCGGTGGT‑TGTGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                  <  1:713127/137‑1 (MQ=255)
                                               acCCGCAAAATA‑CTGCCCCGACGCGACGGACTAAAATTTATTCGTATGGTGGGTTTTCTCCCTCGCCCTGTGGG‑TGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                               <  1:10312989/140‑1 (MQ=255)
                                                      aaaGA‑CTCCCCCGACGCGACGCACGAAAATGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCCCCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                        <  2:4439760/140‑1 (MQ=37)
                                                          a‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCGCCCACACCCAAACCCACACCCCCACCGGAggggggggcgag                                                                                    >  1:6736191/1‑134 (MQ=21)
                                                               cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCGGTCACCCACACCCGAAACCACACCCACAAGGGGGTGGGGGCCGATCGAg                                                                                >  2:5225062/1‑140 (MQ=25)
                                                               cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCACAACCACAAAGGAGAGGGGGccagtgtca                                                                                >  1:9180942/1‑133 (MQ=25)
                                                               cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCCCACCCCAAACCACACCCACCACGGCGAGGGGGcccagcccg                                                                                >  2:1233084/1‑133 (MQ=21)
                                                                   gacgcgaggGACTAAAATTTATTCTTATGGTGGTTTTTCCCCCTCGCCATGTGGT‑TGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcaca                                                                            <  1:8322017/140‑1 (MQ=25)
                                                                    acgcgtcgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCGCAGCGCGCACGGGCACCCACACCCGAACCCGCTCCCGCAAGGGGGCGGGGCCCGCTCgcggccca                                                                           >  2:9889362/1‑133 (MQ=21)
                                                                    acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCACACCCAAAACCACGCCCACAAGGGAGAGGGGACCGATACAGCACaa                                                                           >  1:4011316/1‑140 (MQ=25)
                                                                    acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCACACCCGAACCCACACCCACAAGGGAGAGGGGCCCGAACGAGcccca                                                                           >  2:8035628/1‑136 (MQ=25)
                                                                     cgcgacgTAAGAAAATTTTTTCTTAGGGTTGTTTTTCTCCCTCTCCCTGTGGG‑AGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                                                          <  1:1530045/140‑1 (MQ=37)
                                                                     cgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCGCACCCTAACCCACTCCCACAAGGGGGAGGGGACCGATGGGGCCCaaa                                                                          >  2:8153386/1‑140 (MQ=25)
                                                                         acggactAAAAAGTATTCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                      <  2:8916082/140‑1 (MQ=255)
                                                                         acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGACCCACACCCAAAACCACACCCACCAGGGAGAGGGGGCCGAGAGAGCCCAAAattt                                                                      >  2:29832/1‑140 (MQ=21)
                                                                              cgAAAATGTATTCGTAAGGTTGTTTTTCTCCATATCCGTGTGGG‑AGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                 <  2:2425812/140‑1 (MQ=39)
                                                                              cgAAAAATTATTCGTATGGTGGGTTTTCTCCCTCTCCCTGTGGG‑AGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                 <  1:2958407/140‑1 (MQ=255)
                                                                              cgAAAAAGTATTCGTAAGTTTGTTTTACTCCCTCTCCCTGTGGG‑TGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                 <  1:2958696/140‑1 (MQ=39)
                                                                                aacaaGTATTCGTAAGGTTGTTTTTCGCCCGCTCCCTGGGGG‑AGTGGGCCGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                               <  1:8390543/137‑1 (MQ=38)
                                                                                  aaaGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                             <  2:1845747/140‑1 (MQ=255)
                                                                                   aatttttCGTCCGGTTGGTTTTCCCCCTCCCCCTGTGGG‑TGTGGCCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                            <  1:1171527/135‑1 (MQ=37)
                                                                                   aaGTATTCGTAAGTTGGGTTTTCTCCCTCTCCCTGTGGG‑AGTGGCCCGGGGTGAGTGCATCAGCGCGCACGTTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGTACCGATCGAGCACAAATTGTGAGAGATAtc                                                            <  2:6677123/140‑1 (MQ=38)
                                                                                     gttttCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGTGGGCCGGGGTGAGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                          <  1:7576052/137‑1 (MQ=255)
                                                                                        ttCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCCAAAACACACCCACCAGGGAGCGGGGACCCAACGAGCACAAAAATTGAGAGataaataaat                                                       >  1:4033458/1‑133 (MQ=21)
                                                                                         tcttctggttggttTTCGCCCTCGCCCTGTGGG‑TTTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                      <  2:4181803/134‑1 (MQ=38)
                                                                                         tCGTATGGGTGTTTTTCCCCCACTCCC‑GTGGGTTGTGGGCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                      <  1:6078579/140‑1 (MQ=25)
                                                                                         tCGTAAGGTGGGTTTACTCCCTCGCCGTGTGGG‑AGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                      <  1:8814548/140‑1 (MQ=39)
                                                                                         gcttttggttggttTTCTCCCTCTCCCTGGGGT‑TGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                      <  1:6472678/134‑1 (MQ=38)
                                                                                              tggtgggttTTCCCCCTCGCCGTGTGGT‑TGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                 <  1:1233084/139‑1 (MQ=37)
                                                                                                 ctggtcTCCTCCCGCTCCCTGTGGG‑CGAGGGCCGGGGTGCGGGCCTCACCGCGCACGTACACCCTCCCCCTACCCCTCTCCCTCAAGGGAGCGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                              <  2:9666206/139‑1 (MQ=255)
                                                                                                      tttCTCCCTCTCCGGGGGGG‑AGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                         <  1:7912655/140‑1 (MQ=39)
                                                                                                      tttCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                         <  1:10693297/140‑1 (MQ=255)
                                                                                                      tttCTCCCTCTCCCTGTAGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                         <  1:2656325/140‑1 (MQ=255)
                                                                                                      tatctccctctcccTGTGGG‑AGTGGGCCGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                         <  1:11011267/138‑1 (MQ=255)
                                                                                                        tctccctctcccTGTGGG‑AGAGGGCCGGGGTGAGGGCATCATCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGTGACCGATCTAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                       <  2:7267455/140‑1 (MQ=255)
                                                                                                                           agagGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCACACCCAAAACCACTCCCACCAGGGCGAGGGGACCCATCGAGCACAAATTTTGAGAGATATCATATTATGTATTTTAAAAGTAGAAAAGAAAAAAGGGcgt                     >  2:10918227/1‑138 (MQ=255)
                                                                                                                            gagGGCGGGGTTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                    <  1:1275572/140‑1 (MQ=255)
                                                                                                                            gagGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                    <  1:941847/140‑1 (MQ=255)
                                                                                                                              gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                  <  2:6050156/140‑1 (MQ=255)
                                                                                                                               ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAAACCACGCCCACAAGGGAGAGGGGCCCGCTCGAGCCCAAAATTTTAGAGAAATAATATTATGTAATGTAAAAGAAAAACAGAAAAAAGGGCtttttc                 >  2:2902334/1‑138 (MQ=255)
                                                                                                                                gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACACCCACCAGGGGGAGGGGGCCGATAGAGCACAAAATTTGTGAGATATCATAATAAGTGTTTTTAAAGAAGAACAGAAAATACGGCGTTTtcg                >  2:978431/1‑140 (MQ=255)
                                                                                                                                 ccGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc               <  1:4947396/140‑1 (MQ=255)
                                                                                                                                 ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc               <  1:8526114/140‑1 (MQ=255)
                                                                                                                                 ccGGGGTGAGGGCATCAGCGCGCACGTTCAACCACACCCTAACCCACTCCCACAAGGGAGAGGGGTCCGATCGAGCCCAAATTTTGTGAGATATCTTCTTATGTGTTGTAACAGAATAAACGAAAAAACGGCTTTTgggc               >  2:6340186/1‑140 (MQ=255)
                                                                                                                                    gggTGAGGGCATCAGCGCGCACGTTGACCCTCACCCGAACCCACTCCCACCAGGGAGAGGGGACCCATAGAGCACAAAATTTTAGAGATATATTAATATTTCTTGTAAAAGAAGAAACCAAAAAAGGGCTTTTTCGCCTg            >  2:1389233/1‑140 (MQ=255)
                                                                                                                                      gTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATTTCTTCTTCTTTCTTGTAACATAAGAACAGAAAATCTGGCTTTTGCGCCTGaa          <  2:9095873/140‑1 (MQ=255)
                                                                                                                                             cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg   <  2:5544432/140‑1 (MQ=255)
                                                                                                                                              aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg  <  2:9450594/140‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATA‑CTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGG‑AGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGG  >  NZ_CP025268/164273‑164551

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: ‑