Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP025268164,4120AC28.6% 109.9 / 18.1 105intergenic (+47/‑149)hrpB/mrcBATP‑dependent helicase HrpB/bifunctional glycosyl transferase/transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (33/42);  new base C (9/21);  total (42/63)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
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 105 total aligned reads displayed.

AAAGAGATGAAAGGGCGTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCG  >  NZ_CP025268/164276‑164550
                                                                                                                                            |                                                                                                                                          
aaaGAGATGAAAGGGCGTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCa                                                                                                                                         >  1:1376971/1‑140 (MQ=255)
     gaTGAAAGGGCGTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGGGGGCATCAgcgcg                                                                                                                                    <  2:6433607/140‑1 (MQ=255)
     gaTGAAAGGGCGTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCCTCAgcgcg                                                                                                                                    >  1:3036797/1‑140 (MQ=255)
     gaTGAAAGGGCGTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAgcgcg                                                                                                                                    >  1:5931891/1‑140 (MQ=255)
     gaTGAAAGGGCGTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAgcgcg                                                                                                                                    >  2:1505976/1‑140 (MQ=255)
     gaTGAAAGGGCGTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAgcgcg                                                                                                                                    >  1:2742065/1‑140 (MQ=255)
     gaTGAAAGGGCGTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAgcgcg                                                                                                                                    >  1:2742869/1‑140 (MQ=255)
           tgggCTTTATCCAAAA‑ATTGTGTGGCGGGACCACCGTCAAAATACTCCACCGACGCGACGTACTAAAATTTTTTGTTATGTTTGTTTTA‑CTCCCTCTCCG‑TGTGGGTGTGGGCG‑GGGGTGCGGGCATCAGCGCGCACGtt                                                                                                                              <  1:6094327/139‑1 (MQ=255)
           aGGGCGTTAACCAAAA‑CATGTCTGGCGGGACGACCCGGAAAATACTGCACCTACGCGACGGACTAAAAAGTTTTCGTACGTTTGTTTTT‑CTCCCTCTCCC‑TGTGGGTGTGGGCC‑GGGGTGAGGGCATCAGCGCGCACGtt                                                                                                                              <  2:5600126/140‑1 (MQ=25)
               cGTTATCCCCAA‑CATGTATGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCCTCAGCGCGCACGTtcacc                                                                                                                          >  1:1151028/1‑140 (MQ=255)
               cGTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCACCA‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                          <  2:4646792/140‑1 (MQ=255)
                gTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                         >  2:3631410/1‑140 (MQ=255)
                gTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                         >  2:3631138/1‑140 (MQ=255)
                 ttttCCCAAA‑CTTTTCTGGACGGCCGACCCGGCAAAAAGTGCACCTTCGCGACGGACTAAAAAGTATTAGTAAGGTTGTATTA‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GTGGTGCTGGCATCAGCGCGCACGTtcaccct                                                                                                                        <  2:5792134/137‑1 (MQ=25)
                 ttATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCCGCGCGCACGTtcgcccg                                                                                                                        >  1:5005950/1‑139 (MQ=255)
                 ttATCCCAAA‑CATGTCTGGCCGGACGACCCGGAAAATACTGCACCGAGGCGACGGACGAAAATTTATTCTTAATTTTGGTTTT‑ATCCCTCTCCA‑TGTGGGTGAGGGCC‑GGGGTGTGGGCATCAGCGCGCACGTtcaccct                                                                                                                        <  2:2734140/140‑1 (MQ=255)
                 ttACCCCAAC‑CAGCTCGGGCCGGACGACCCGCCAAATCCTCCACCCACGCGCCGGACGAAAAAGTATTCGTAAGTTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                        <  2:5792120/140‑1 (MQ=25)
                  ttccccAAACCATGTCGGGCCGCACGACCCGGCAACCACTGCACCGCCGCGCCGCACTAAAACGTATGCTTAAGGTGGGTTTT‑CTCCCACACCC‑TGTGGGTGAGGGCC‑GGGGTGCGGCCCTCAGCGCGCACGTtcaccct                                                                                                                        <  2:1369945/137‑1 (MQ=255)
                   accccAAA‑CTTGTGGGCCCGGCCCCCCGGGAAAATTCGTCACCGCCGCGACGGACTAAAATTTATTC‑TTATGTTGGTTTTACTCCCCCTCAG‑TGTGGTTGCGGGCG‑GGGGTGAGGCCCTCAGCGCGCACGTtcaccctca                                                                                                                      <  2:1151028/138‑1 (MQ=255)
                   aTCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTggacccaca                                                                                                                      >  1:251706/1‑140 (MQ=255)
                      ccAAA‑CTTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGTGTGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                   <  2:1437596/140‑1 (MQ=255)
                        aaa‑CTTGTCTGGCGGGACGACCCGTCAAATACTGCACCGACGCTACGTACTAAAAAGTTTTCTTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGTGGGCC‑GGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTa                                                                                                                 <  2:18585/140‑1 (MQ=37)
                          accaggTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGTTTTT‑CTCCCACACCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                <  2:5255506/135‑1 (MQ=38)
                             atgacgggCCGGACGCCCCGGAAAATCCTCCACCGACGCGACGGACGAAAAAGTATACGAATGTTTGTTTTA‑CGCCCTCACCC‑TGTGGGAGTGGGCG‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                             <  1:1372771/134‑1 (MQ=25)
                             aTGTCGGGCGGGACGACCGGTCACATCCTGCCCCGCCGCGCCGGACGAAAATGTTTTCGTATGTTTGTTTTC‑CTCCCTCGCCC‑TGTGGGTTCGGGCG‑GGGGGGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAccc                                                                                                             <  2:3855941/140‑1 (MQ=255)
                               ggcgggccggccCCCCCGGAAAATCCTGCCCCGCCGCGACGGACGAAAAAGTATTCTTACGTTTGGTTTT‑CTCCCTCTCCC‑TGGGGGTGTGGGCC‑GGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                           <  2:2051317/136‑1 (MQ=25)
                               gTCTGGCCGGACGCCCCGGCAAATACTGCACCGCCGCGACGGACGAAAAAGTATGCGGAAGGTGGGTGTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCGCCCTAACCctc                                                                                                           >  2:6380714/1‑140 (MQ=25)
                               gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCtaaacctaaaccaa                                                                                                           >  2:3220076/1‑138 (MQ=255)
                               gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAaccaa                                                                                                           >  2:3170253/1‑138 (MQ=255)
                               gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCCCCCCAAcccac                                                                                                           >  1:1588345/1‑138 (MQ=37)
                               gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCGAAacccc                                                                                                           >  1:6338235/1‑138 (MQ=38)
                                tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCCTCCGCGCGCACGGTCACCCACACCCTAACCctct                                                                                                          >  2:1583363/1‑140 (MQ=255)
                                tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                          >  2:1583147/1‑140 (MQ=255)
                                tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAAaccaca                                                                                                          >  2:1260722/1‑137 (MQ=37)
                                 gtGGCCGGACTCCCCGGAAAATCCTCCACCGACGCGAGGGACTAAAATTTTTTCTTAAGGTTGGTTTT‑CTCCCTCTCCC‑TTTGGGTGTGGGCC‑GGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                         <  1:2765023/139‑1 (MQ=255)
                                 cTGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                         <  2:5434392/140‑1 (MQ=255)
                                 cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAAaccacac                                                                                                         >  1:4627953/1‑136 (MQ=37)
                                        gcccacCCGCCAAATACTGCACCGACGCGCCGGACGAAAATGTATTCTTAAGGTTGGTTTT‑CCCCCTCTCCC‑TGGGGGTGTGGGCG‑GGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                  <  1:380422/136‑1 (MQ=255)
                                        gacgacCCGGAAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGAAAGGTTGGTTTT‑CTCCCTCGCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                  <  2:3706845/140‑1 (MQ=255)
                                         acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTCAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGTGGGCCTCCGCGCGCACGTTGACCCACACCCAAACCCTCTCCCTCAAgg                                                                                                 >  1:4302110/1‑140 (MQ=38)
                                         acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCCAAACCACACCCTCACgg                                                                                                 >  1:1567270/1‑140 (MQ=25)
                                            ccccgcaaaaTACTGCCCCGCCGCGACGTACTAAAATTTATTCTTTTGTTTGTTTTA‑CACCCTCTCCT‑TGGGTGTTTGGGCG‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                              <  1:804647/133‑1 (MQ=255)
                                            ccccGGCAAATCCCGCACCCCCGCGACGTCCGAAACACTATTCGTCCGGTTGGTTTA‑CTCCCTCGCCC‑TGTGGGTGAGGGCC‑GGGTTGCGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                              <  2:2913652/139‑1 (MQ=255)
                                            ccccGGCAAATACTGCCCCCTCGCGTCGCACTAAAATGTTTTCGTATGTTTGTTTTT‑CCCCCCCTCCC‑TGGGGGAGCGGGCC‑GGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                              <  2:1997629/139‑1 (MQ=255)
                                            acCCGGCAAATCCTGCACCGCCGCGACGGACGAAAATGTATGCTTAAGGTGGGTTTC‑CCCCCCCTCCC‑TGGGGGCGCGGGCC‑GGGGGGCGGGCCTCAGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                              <  2:2372550/140‑1 (MQ=255)
                                            acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                              <  2:5373497/140‑1 (MQ=255)
                                            acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                              <  2:3688881/140‑1 (MQ=255)
                                                 gCAAATACTCCACCGACGCGCCGGACTAAAAAGTATTCGTATGGTGGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggt                                                                                         <  2:576707/140‑2 (MQ=37)
                                                  cAAATACTGCACCTACGCGACGGACGAAAAATTAATCGTAAGGTTGGTTTT‑CTCACTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGAACGTTAACCCACAAACAAAACCAAACCCAAAAGGGAGAGGGGa                                                                                        >  2:4298661/1‑140 (MQ=21)
                                                       aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCCTCAGCGCGCACGTGCACCCAAACCCGAAACCACTCCCGCAAGGGAGGGGGGGCCGAt                                                                                   >  2:1250594/1‑140 (MQ=25)
                                                          gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCCTCAGCGCGCACGTTGACCCGCACCCAAACCCACACCCACCAGGGAGAGGGGGcccaaaga                                                                                >  2:2536573/1‑134 (MQ=25)
                                                           cACCGACGCGTCGGACGAAGACGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGGGAGGGCCTCAGCGCGCACGTTCTCCCGCCCCCTCACCCTCTGCCTCTAGGGCGTGGGGACCGCTCGAg                                                                               >  2:1493562/1‑140 (MQ=25)
                                                           cACCGACGCGCGGTACGAAAATTTTTTCTTACGGTTGTTTTC‑ACCCCACTCCC‑GGTGGTAGCGGGCG‑GGGGTGTGGCACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                               <  1:5166179/140‑1 (MQ=19)
                                                           cACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCACACCCACCAGGGGGTGGGGACCGATCGAg                                                                               >  2:4693014/1‑140 (MQ=37)
                                                            aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGGGGACCCGCACCCCAACCCACTCCCGCCACGGTGTGGGGGCCGCTGGAGc                                                                              >  1:5215570/1‑140 (MQ=21)
                                                             ccgccgCGACTGCCGAATAAGTACTCGACTGGTGGGAGCT‑CTCCCTCCCCC‑TGTGGGATTGGGCC‑GGGGTGTGGCCCTCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                             <  1:5751001/140‑1 (MQ=19)
                                                             ccGACGCGAGGGAATAAAAATTATTCGTATGGTTGGTTTT‑CTCCCTCACCC‑TGTGGGTGTGGGCG‑GGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                             <  1:5121406/140‑1 (MQ=38)
                                                                acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCCTCAGCGCGCACGTGGCACCGCACCCAAAACCACTACCACAAGGGGGGGGGGGCCGATCGGGCACaa                                                                          >  2:2880752/1‑140 (MQ=21)
                                                                acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCTCGTTGACCCACCCCCACACCCACACCCGAAAGGGAGAGGGGTCCCATCGAGCACaa                                                                          >  2:3300828/1‑140 (MQ=25)
                                                                acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGGGGCCCCGCCCCCCAACCCACTCCCACCCGGGGGAGGGGCCCGGTCGCGCCCaa                                                                          >  1:1520705/1‑140 (MQ=18)
                                                                       gtacgAAAATTTTTTCGTTTGGTTGTTTTC‑CCCCCTCGCGG‑GGGGGGTGGGGGCG‑GGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                   <  1:2536573/138‑1 (MQ=19)
                                                                       gcccgAAAACGAATGCGTAAGGTTGGTGTC‑CCCCCACCCCG‑GGGGGGTGCGGGCC‑GGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGTGCACAAATTTTga                                                                   <  1:5654214/137‑1 (MQ=19)
                                                                          ctaaaaaGTTTGCGTAGGTTTGGTTTT‑CTCCCTCTCCC‑TGTGGTTGTGGGCG‑GGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                <  1:3500122/138‑1 (MQ=37)
                                                                          ctaaaaaGTATTCTTATGTTTGGTTTT‑C‑CCCCGGCCC‑CGTGGGTCTTGGCCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGTGGGGACCGATCGAGCACAAATTTTgagag                                                                <  1:351722/138‑1 (MQ=255)
                                                                          cgAAAATTTATTCGTAAGGTTGTTTTT‑CTCCCTCTCCC‑TGGGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                <  2:5015756/140‑1 (MQ=255)
                                                                          cgAAAAAGTATTCTTACGGTTGGTTTT‑CTCCCTCGCCC‑TGTGGGAGTGGGCC‑GGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                <  2:3110687/140‑1 (MQ=39)
                                                                          cgAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACGAATTTTgagag                                                                <  2:2742065/140‑1 (MQ=255)
                                                                          cgAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑GGGGGGAGAGGGCC‑GGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                <  2:2742869/140‑1 (MQ=255)
                                                                            aaaaaGTATTCGTTAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAACGAGCACAAATTTTGAGAGat                                                              <  2:6227612/140‑1 (MQ=255)
                                                                            aaaaaGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                              <  2:6228042/140‑1 (MQ=255)
                                                                            aaaaaGTAGTCGTACGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGTGTGGGCC‑GGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                              <  1:3631410/140‑1 (MQ=38)
                                                                              aaaGTAGTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                            <  2:1389712/140‑1 (MQ=255)
                                                                               atgttttCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGTGAGGGCC‑GGGGTGAGGGCATCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                           <  2:5265601/135‑1 (MQ=255)
                                                                               aaGTATTCGTAAGGTTGGTTTA‑CTCCCTCACCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                           <  2:5170238/140‑1 (MQ=255)
                                                                                   aTTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTGACCCACACCCAAACCCGCTCCCCAAAGGGAGGGGGGTCCGATCGCGCCCAAAATTTTTGTGAGAtcttat                                                       >  2:152389/1‑138 (MQ=21)
                                                                                   aTTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCTAACCCACAACCTCAAGGGTGAGGGGACCGGTCGAGCCCCAATTTTGAGAGATAtattat                                                       >  2:3488256/1‑138 (MQ=25)
                                                                                   aTTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCTAACCCACACCCACAAGGGAGAGGGGACCGAACGAGCACAAAATTTGAGAGataaattat                                                       >  2:4491739/1‑134 (MQ=25)
                                                                                     tCGTATGTTTGGTTTC‑CCCCCTCGCCA‑GGGGGGTGTGGGCC‑GGGGTTTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                     <  1:1461329/140‑1 (MQ=37)
                                                                                     tCGTATGGTTGGTTTT‑CTCCCTCTCCC‑TGGGGGAGTGGGCC‑GGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                     <  2:2019530/140‑1 (MQ=38)
                                                                                     tCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGTGGGCG‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                     <  2:3284524/140‑1 (MQ=255)
                                                                                     gcGTACGGTTGGTTTT‑CTCCCTCACCC‑GGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                     <  1:6289647/139‑1 (MQ=255)
                                                                                      tctttttttggttTT‑CCCCCCCCCCCGTGGGGGTGGGGGCG‑GGGTTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                     <  1:2876194/133‑1 (MQ=25)
                                                                                          aGGTTGGTTTT‑CTCCCTCTCCC‑TGGGGGAGAGGTCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTACTGTCt                                                <  2:6057198/140‑1 (MQ=255)
                                                                                               gttttt‑CCCCCTCTCCC‑TGTGGTAGAGGGCC‑GGGTTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                           <  1:4071105/138‑1 (MQ=39)
                                                                                                     actccctctccc‑tGTGGGAGTGGGCC‑GGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                      <  1:1831163/139‑1 (MQ=255)
                                                                                                          ctctcca‑gGGGGGAGAGGGCC‑GGGGGGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                 <  1:1683829/140‑1 (MQ=39)
                                                                                                           tctccc‑tGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTGACCCTCACCCTAACCCACTCCCGCAAGGGAGAGGGGGCCGATCGAGCACAAATTTTGTGAGATATATTGTTATGTGTTGTAACAGAAAAAAAg                                >  2:4250577/1‑140 (MQ=37)
                                                                                                                  tgtgGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTGCCCCTCACCCGAACCCTCTCCCTCAAGGGAGTGGGGACCGATCGAGCACAAAATTTGAGAGATATATTCTTATGTCTTGTTACAGAAGTAAAGAAAAta                          >  2:1684551/1‑139 (MQ=38)
                                                                                                                        tgtgggCG‑GGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                    <  1:5834996/137‑1 (MQ=255)
                                                                                                                         gcgggCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                   <  1:4616570/138‑1 (MQ=255)
                                                                                                                         gagGGCC‑GGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                   <  1:311921/140‑1 (MQ=255)
                                                                                                                           gggCG‑GGGGTGAGGGCCTCAGCGCGCCCGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                 <  1:5583587/140‑1 (MQ=255)
                                                                                                                           gggCC‑GGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                 <  2:2151506/140‑1 (MQ=255)
                                                                                                                            ggCC‑GGGGTGAGGGCATCAGCGCGCACGTTGACCCTCGCCCTATACCTCGCCCTCAAGGGGGTGGGGTCCGTTAGAGCACCAATTTTGTGAGATATATTAGTAGGTCTTGTTACAGAAGGACAGAAAATCGGGCTTTTgg                >  2:1866758/1‑139 (MQ=255)
                                                                                                                            ggCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACTCCCACAAGGGAGAGGGGACCCAACGAGCACAAATTTTGAGAGAAATCATCATATTTAATGTAACAGAAGAACAGAAAAACGGGCtttttc                >  2:6237987/1‑138 (MQ=255)
                                                                                                                             gCC‑GGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg               <  1:902608/140‑1 (MQ=255)
                                                                                                                              cc‑GGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCGTACCCACGCCCACCAGGGAGTGGGGGCCGCTCGAGCCCAAATTTTGAGATATATATTCTTATGTCTTGGAACAGAAGAACCGGAAATAGGGCGTTTgcgc              >  2:4640404/1‑140 (MQ=255)
                                                                                                                              cc‑GGGGAGAGGCCCTCAGCGCGCACGTTCACCCTCACCCGCCCCCTCACCCTCAAGGGAGAGGGGACCGATCGAGCACACATTTCGAGAGCTATCTTCTTCGGTCTGGGAACAGACGAACAGACACTCGGGCTTTTgcgc              <  2:1881204/140‑1 (MQ=255)
                                                                                                                                 ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt            <  2:1003691/140‑1 (MQ=255)
                                                                                                                                           cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGCCCACAAATTATGAGAGATAACTTCTACTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg  <  2:1423755/140‑1 (MQ=39)
                                                                                                                                            |                                                                                                                                          
AAAGAGATGAAAGGGCGTTATCCCAAA‑CATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCC‑TGTGGGAGAGGGCC‑GGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCG  >  NZ_CP025268/164276‑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: ‑