Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP025268164,4320TA23.8% 312.4 / 11.8 135intergenic (+67/‑129)hrpB/mrcBATP‑dependent helicase HrpB/bifunctional glycosyl transferase/transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (4/89);  new base A (29/0);  total (46/89)
Fisher's exact test for biased strand distribution p-value = 4.24e-24
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 160 total aligned reads displayed.

TTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGA  >  NZ_CP025268/164293‑164570
                                                                                                                                             |                                                                                                                                           
ttATCCCACACATGTCTGGCCGGACGCCCCGGCAAAGACTGCACCGACGCGCCGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcaccct                                                                                                                                             <  2:10561028/140‑1 (MQ=37)
ttATCCCAAACTTGTCTGGCCGGACGACCGGGAAAATACGGCACGGACGCTACGTACTAAAAAGTATTCGTATGTTTGTTTTTCTCCCTCTCAATGTGGTTTAGGGCGGGGTTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                             <  2:7130950/140‑1 (MQ=255)
ttATCCCAAACATTTCTGGCCGGACGACCCGTAAAATCCTGCACCTACGCGACGGACTAAAAAGTTTTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                             <  2:5805553/140‑1 (MQ=37)
ttATCCCAAACATGTCTGGCCGGACGCCCCGGCAAATACTGCACCGACGCGACGGACGAAAATTTATTCGTACGGTTGGTTTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                             <  2:1652171/140‑1 (MQ=38)
ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                             <  2:8739504/140‑1 (MQ=255)
ttATCCCAAACATGTCTGGCCGGACGACCCGGAAAATCCGGCACCGACGCGACGGACGAAAATGTATTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGTTGAGGGCGGGGGTGTGGGCCTCAGCGCGCACGTtcaccct                                                                                                                                             <  2:3551524/140‑1 (MQ=37)
cgagcccAATCCTGTCTGCCCGGCCGACCCGGACAATCCTGCACCGACGCGACGCCCGATAACGTATTCGACTGGTGGGTTGTCTCCCTCTCCCTCTGGGAGAGGGCCGGGTTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                             <  2:9448417/136‑1 (MQ=16)
 taacccAAACATTGTTTGGCCGGACGCCCCGTAAAATACTGCACCGACGCGCCGTACTAAAATGTTTACTTAAGTTTGTTTTTCTCCCACTCCCTGTGGTAGTGGGCCGGGGTGAGGCCATCAGCGCGCACGTtcaccct                                                                                                                                             <  2:3442456/137‑1 (MQ=255)
     ccAAAAATTTCTGGCCGTACGACCGGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTTAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                                        <  2:1544560/140‑1 (MQ=38)
      cAAACATGTCTGGCGGGCCGCCCCGGCAAATACTCCACCGCCGCGAGGGACTAAAATGTTTTCGAACGGTTGGTTTTCACCCACGCCGTGGGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTcaccctcaccct                                                                                                                                       <  2:4320594/140‑1 (MQ=255)
      cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTgggcccccacccc                                                                                                                                       >  1:5847238/1‑139 (MQ=37)
      cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTcaccctcaccct                                                                                                                                       >  2:10794839/1‑140 (MQ=255)
        aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccccacccaca                                                                                                                                     >  1:8875112/1‑137 (MQ=255)
         aaattGTCGGGCCGTACTACCCGGAAAATACTGAACCTCCGCTTCGTACTAAAATGTTTTCGTAGGTTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                                     <  1:7534278/136‑1 (MQ=16)
         aCATGTCTGGCCGGCCTACCCGGCAAATCCTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                                    <  2:2820280/140‑1 (MQ=38)
          cATGTGTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                                   <  2:4038173/140‑1 (MQ=255)
           aTGTCTGGCCGTACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                                                  <  2:1534945/140‑1 (MQ=255)
            ttGTCGGGCGGGACGACCCGGAAACCAGTGCCCGGACGCGAGGTCCGAAAATGTTTGCTTACGGTTGGTTTACTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                                                  <  1:3031067/139‑1 (MQ=255)
              gacTGGCCGTACCCCCCGCCAAATACTGCACCTACGCTACGTACTAAAAATTATTCTAACGGTTGTTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                <  1:2880229/138‑1 (MQ=21)
              gTCTGGCGGGACGACCCGTCAAATACTGCACCGACGCGACGGACGAAAATTTATTCTTATGGTTGTTTTACTCCCTCACCCTGTGGGCGTGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                <  1:9557765/140‑1 (MQ=255)
              gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                <  2:10706072/140‑1 (MQ=255)
              gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAAcccac                                                                                                                                >  2:9247136/1‑138 (MQ=255)
              gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAAaccac                                                                                                                                >  1:10490356/1‑138 (MQ=255)
              gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCCaaaacaa                                                                                                                                >  1:6550746/1‑135 (MQ=37)
              gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAcccaa                                                                                                                                >  1:10114715/1‑138 (MQ=255)
              gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCGCCCGCACCCGAAcccac                                                                                                                                >  1:7423684/1‑138 (MQ=38)
              gTCTCGCCGGCCGACCCGGCAAATACTGCACCGACGCGACGTCCGAAAAATTATTCTAATGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                <  1:10536375/140‑1 (MQ=37)
              gTCGGGCCGGACGCCCCGCCACATACTGCACCTACGCGACGGCCGAAAAATTAGTCGAATGTTTGGTTTGCCCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                <  1:7688529/140‑1 (MQ=25)
               tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCcact                                                                                                                               >  2:10344839/1‑140 (MQ=255)
               tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCTCCCACACCCGAAACcact                                                                                                                               >  1:3918838/1‑140 (MQ=37)
               tCTGCCCGGACGCCCCGCCAAATACTGCACCGACGCGACGCACGAAAAAGTATTCGTAAGGTTGGTTTGCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCctct                                                                                                                               <  2:6379155/140‑1 (MQ=255)
                cgggCCGGACGACCCGCCACCTCCTGCCCCGACGCGCCGGACGAAAATTTATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                              <  1:1840354/138‑1 (MQ=25)
                cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCACACCCCAACCCactc                                                                                                                              >  2:3732790/1‑140 (MQ=255)
                 tGGCCGGAGGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                             <  2:9591777/140‑1 (MQ=255)
                   gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCGTGCACCCGCACCCGAACCCACAccca                                                                                                                           >  2:2762805/1‑139 (MQ=37)
                       gccgacCCTGAAAATCCTCCACCGCGGCGAGGTACTAAAATTTATTCTTTTGGTTGGTTTTCTCCCTCTCCCGGTGGTAGTGGGCGGGGGTGTGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                       <  1:5544691/138‑1 (MQ=255)
                       gacgccCCGGCAAATACGGCACCGCCGCGCCGGACGAAAAAGTATTCGAACGGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                       <  2:3021099/140‑1 (MQ=38)
                       gaccacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATCCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                       <  1:7592115/140‑1 (MQ=255)
                           ccccGGCAAACACCCCCCCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCCCCCGCGCCCGGGGGGCGCGGGCCGGGGTGCGGGCAGCAGCGCGCACGTTCACCCTCACCCTCACCCTCTCCCTCAAGGgag                                                                                                                   <  2:8923856/139‑1 (MQ=255)
                           acCCGCAAAATACTGCCCCGACGCGACGGACTAAAATTTATTCGTATGGTGGGTTTTCTCCCTCGCCCTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                                                   <  1:10312989/140‑1 (MQ=255)
                            cccgcaaaaTCCTCCACCGCCGCGCCGGACCAAAATTTATTCGTCTGGTTGTTTTCCCCCCGCCCCGTGTGGGTGGGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                  <  1:10794839/134‑1 (MQ=255)
                             ccGGCAAATACTGGACCGACGCGACGGACTAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                 <  2:1852130/140‑1 (MQ=255)
                             ccGGCAAATACTCCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                 <  2:9066907/140‑1 (MQ=255)
                                  aaaTCCTGCACCGCCGCGCCGGACCAAAATGTATGCTTAAGGTTGGTTTCCCCCCACTCCCTGTGGCCGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                            <  1:9402301/140‑1 (MQ=255)
                                  aaaGACTCCCCCGACGCGACGCACGAAAATGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCCCCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                            <  2:4439760/140‑1 (MQ=37)
                                      aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCGCCCACACCCAAACCCACACCCCCACCGGAggggggggcgag                                                                                                        >  1:6736191/1‑134 (MQ=21)
                                      aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGGCCGCACCCAAACCCCAGCCCCCAAggggggggggggcgcg                                                                                                        >  1:7958699/1‑134 (MQ=21)
                                          cACCGACGCGACGTACTAAAATTTATTCGTATGGTTGTTTTTAACCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                    <  1:1481904/140‑1 (MQ=37)
                                          cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAACCACAACCACAAAGGAGAGGGGGccagtgtca                                                                                                    >  1:9180942/1‑133 (MQ=25)
                                          cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCCCACCCCAAACCACACCCACCACGGCGAGGGGGcccagcccg                                                                                                    >  2:1233084/1‑133 (MQ=21)
                                              gacgcgtcgGACGACAAAGTATTCGTAAGGTTGTTTTACTCCCTCTCCCTGTGGGAGAGCGCCGGGTTGAGGGCATCACCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcaca                                                                                               <  2:10455978/140‑1 (MQ=38)
                                              gacgcgaggGACTAAAATTTATTCTTATGGTGGTTTTTCCCCCTCGCCATGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcaca                                                                                               <  1:8322017/140‑1 (MQ=25)
                                               acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCACACCCAAAACCACGCCCACAAGGGAGAGGGGACCGATACAGCACaa                                                                                              >  1:4011316/1‑140 (MQ=25)
                                               acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCACACCCAAACCCACTACCTCAAGGGAGTGGGGACCGATAGAGCACaa                                                                                              >  2:5778719/1‑140 (MQ=37)
                                               acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCACACCCGAACCCACACCCACAAGGGAGAGGGGCCCGAACGAGcccca                                                                                              >  2:8035628/1‑136 (MQ=25)
                                                cgcgacgTAAGAAAATTTTTTCTTAGGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                                                                             <  1:1530045/140‑1 (MQ=37)
                                                cgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCGCACCCTAACCCACTCCCACAAGGGGGAGGGGACCGATGGGGCCCaaa                                                                                             >  2:8153386/1‑140 (MQ=25)
                                                   gacggacgAAAACTTATGCGTACGGTGGGTTTACTCCCTCTCCTTGTGGGAGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCCCCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                          <  2:2412283/140‑1 (MQ=37)
                                                    acggactAAAAAGTATTCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                         <  2:8916082/140‑1 (MQ=255)
                                                    acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGACCCACACCCAAAACCACACCCACCAGGGAGAGGGGGCCGAGAGAGCCCAAAattt                                                                                         >  2:29832/1‑140 (MQ=21)
                                                    acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGGTCACCCGCACCCGAACCCACTCCCGGGAGGGGGGGGGGGGCGGTGGGGCACAAAtttg                                                                                         >  1:10590127/1‑139 (MQ=21)
                                                      ggacgAAAATGTATTGGTAAGGTGGGTTTTCGCCCTCTCCCTGGGGTTGCGGGCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                       <  2:4087631/140‑1 (MQ=37)
                                                       gactAAAATTTATTCTTATGTTTGTTTTTCTCCATCGCCCTGTGGGTGTGGGCCGGGGTCTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                                      <  1:8422513/140‑1 (MQ=37)
                                                         cgAAAATGTATTCGTAAGGTTGTTTTTCTCCATATCCGTGTGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                                    <  2:2425812/140‑1 (MQ=39)
                                                         cgAAAAATTATTCGTATGGTGGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                                    <  1:2958407/140‑1 (MQ=255)
                                                         cgAAAAAGTATTCGTAAGTTTGTTTTACTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                                    <  1:2958696/140‑1 (MQ=39)
                                                              aatttttCGTCCGGTTGGTTTTCCCCCTCCCCCTGTGGGTGTGGCCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                               <  1:1171527/135‑1 (MQ=37)
                                                              aaGTATTCGTAAGTTGGGTTTTCTCCCTCTCCCTGTGGGAGTGGCCCGGGGTGAGTGCATCAGCGCGCACGTTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGTACCGATCGAGCACAAATTGTGAGAGATAtc                                                                               <  2:6677123/140‑1 (MQ=38)
                                                                  aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACCCCCCAACCCACCCCCCAAAGGGAGAGGGGCCCGCAAGCGCCCAAAATTTTTGAGatataatat                                                                           >  2:5285518/1‑135 (MQ=12)
                                                                   ttCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCCAAAACACACCCACCAGGGAGCGGGGACCCAACGAGCACAAAAATTGAGAGataaataaat                                                                          >  1:4033458/1‑133 (MQ=21)
                                                                    tCGTATGGGTGTTTTTCCCCCACTCCCGTGGGTTGTGGGCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                                         <  1:6078579/140‑1 (MQ=25)
                                                                      gaaaGGTTGGTTATCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTg                                                                       <  2:11031402/138‑1 (MQ=38)
                                                                         tggtgggttTTCCCCCTCACCGTGGGGGAGGGGGGGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                                    <  1:1233230/139‑1 (MQ=37)
                                                                            ttggttTTCTCCCACTCCCTGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                                 <  1:1868059/140‑1 (MQ=255)
                                                                            ctggtcTCCTCCCGCTCCCTGTGGGCGAGGGCCGGGGTGCGGGCCTCACCGCGCACGTACACCCTCCCCCTACCCCTCTCCCTCAAGGGAGCGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                                 <  2:9666206/139‑1 (MQ=255)
                                                                                 tttCTCCCTCTCCGTGGGGGTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                            <  1:4131939/140‑1 (MQ=38)
                                                                                 tatctccctctcccTGTGGGAGTGGGCCGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                            <  1:11011267/138‑1 (MQ=255)
                                                                                   tctccctctcccTGTGGGAGAGGGCCGGGGTGAGGGCATCATCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGTGACCGATCTAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                          <  2:7267455/140‑1 (MQ=255)
                                                                                   tctccctcgcccGGGGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                          <  2:5151742/140‑1 (MQ=255)
                                                                                                 tgGGAGAGGGCCGGGGTGAGGGCATCAGCGGGCACGTTGACCCGAAACCCAAACCAAAACCACCATGGAGTGGGGTCCCATCGAGCCACAAATTTTTGAGTTATATTATTATTTGTTGTAAAAGAATAAAAGAAAATAgg                                            >  2:3175787/1‑140 (MQ=255)
                                                                                                      agagGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCACACCCAAAACCACTCCCACCAGGGCGAGGGGACCCATCGAGCACAAATTTTGAGAGATATCATATTATGTATTTTAAAAGTAGAAAAGAAAAAAGGGcgt                                        >  2:10918227/1‑138 (MQ=255)
                                                                                                       gagGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                                       <  1:941847/140‑1 (MQ=255)
                                                                                                         gggCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                                     <  1:6323005/140‑1 (MQ=255)
                                                                                                          ggCCGGGGTGAGGGCATCAGCGCGCACGTTCGCCCCCAACCTAAACCACACCCTCAAGGGCGAGGGGACCGGTCCAGCACAAAATTTGAGTGATATATTATTGAGTAGTGTGAAAGAAGAACAGAACAGAGGGCTTTTgg                                    >  2:6538715/1‑139 (MQ=255)
                                                                                                           gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACACCCACCAGGGGGAGGGGGCCGATAGAGCACAAAATTTGTGAGATATCATAATAAGTGTTTTTAAAGAAGAACAGAAAATACGGCGTTTtcg                                   >  2:978431/1‑140 (MQ=255)
                                                                                                            ccGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                  <  1:4947396/140‑1 (MQ=255)
                                                                                                            ccGGGGTGAGGGCATCAGCGCGCAGGTGGACCCACACCCTAACCCACACCCACAAGGGAGTGGGGACCCATCGAGCCCCACATTTGAGAGCTATATTAATCTGTCTTTGAACCGAAGGACAGACAAAACGGCGTTTgcgc                                  >  2:11179659/1‑140 (MQ=255)
                                                                                                            ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                  <  1:8526114/140‑1 (MQ=255)
                                                                                                            ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACACCCCCAAGGGAGAGGGGACCCATAGCGCACAAATTTTGAGAGAAAAAATATTATGTCATGTAACAAAACCACACCAAAAAGGGCGTTTgcgc                                  >  2:10683271/1‑140 (MQ=255)
                                                                                                            ccGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCTAAACCACAACCACAAGGGAGTGGGGGCCCATCGAGCCCAAATTTTGTGTGATATATAATTCTGTCATTTTAAAGAAGAACAGCAAATCGGGCGTTTgggc                                  >  2:1815268/1‑140 (MQ=255)
                                                                                                               gggTGAGGGCATCAGCGCGCACGTTGACCCTCACCCGAACCCACTCCCACCAGGGAGAGGGGACCCATAGAGCACAAAATTTTAGAGATATATTAATATTTCTTGTAAAAGAAGAAACCAAAAAAGGGCTTTTTCGCCTg                               >  2:1389233/1‑140 (MQ=255)
                                                                                                                 gTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATTTCTTCTTCTTTCTTGTAACATAAGAACAGAAAATCTGGCTTTTGCGCCTGaa                             <  2:9095873/140‑1 (MQ=255)
                                                                                                                      gccctccgcgcgcACGCCCCCCCTCACACACACCCTCGCCCTCCAGCGAGAGGCGCCCGATCGCGCACAAATTTTGAGAGCTATCTTCTTCTGTCATGTAACAGAAGAACAGACCATCGGGCTTTTGCGCCTGAATATTg                        <  2:6959803/133‑1 (MQ=255)
                                                                                                                        cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                      <  2:5544432/140‑1 (MQ=255)
                                                                                                                         ctCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTACTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg                     <  2:6119866/139‑1 (MQ=255)
                                                                                                                         aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg                     <  2:9450594/140‑1 (MQ=255)
                                                                                                                                  cACGTTCACCCTCACCCTACCCCTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACCGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc            <  1:6990319/140‑1 (MQ=255)
                                                                                                                                    cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt          <  1:11391189/140‑1 (MQ=255)
                                                                                                                                      ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATgg        <  2:6736191/140‑1 (MQ=255)
                                                                                                                                        caccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGAGATCTTCTTCTGTCTTGTAACCGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGcc      <  1:7908938/140‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                           
TTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGA  >  NZ_CP025268/164293‑164570

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: