Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273360,7590TG38.5% 117.9 / 18.7 66T329P (ACC→CCC) lacIDNA‑binding transcriptional repressor LacI
Reads supporting (aligned to +/- strand):  ref base T (8/32);  new base G (25/0);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 2.65e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGA  >  NZ_CP009273/360629‑360896
                                                                                                                                      |                                                                                                                                           
gcgcAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGGGGtttttt                                                                                                                                     >  1:510138/1‑138 (MQ=255)
    aGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTttttctttt                                                                                                                                 <  2:523306/139‑1 (MQ=255)
    aGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTttttctttt                                                                                                                                 <  2:338894/139‑1 (MQ=255)
    aGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTttttctttt                                                                                                                                 <  1:602504/139‑1 (MQ=255)
           gTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGT‑TTGG‑GATTGGGGG‑CCCCGGGGGGTTTTTTTTTTTAccgg                                                                                                                           >  2:312132/1‑137 (MQ=255)
           gTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTTCGGGTTTGGGC‑CCAGGGGGGGTTTTTTTTTTCAcacg                                                                                                                           >  2:268167/1‑136 (MQ=255)
           gTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGG‑GTCTTGGGG‑CCCCGGGGGGTTTTTCTTTTTAcccgc                                                                                                                          >  1:306619/1‑136 (MQ=255)
           gTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGG‑ACCAGGGGGGGTTTTCTTTTCACCCGt                                                                                                                          >  2:249394/1‑139 (MQ=255)
              cTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTgag                                                                                                                       <  1:133350/139‑1 (MQ=255)
                 aaTCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                      <  2:105851/137‑1 (MQ=255)
                 aaTCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTgaga                                                                                                                      <  1:250126/137‑1 (MQ=255)
                      cgcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTTC‑GGTTTTGGG‑ACCAGGGGGGGTTTTTTTTTTTCCAGTGGGAGGGGGaa                                                                                                               >  2:309490/1‑139 (MQ=255)
                      cgcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTTGGC‑ACCCGGGGGGGTTTTTTTTTTACCAGGGAGACgggga                                                                                                                >  1:285242/1‑136 (MQ=255)
                      cgcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTTGGC‑ACCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCaa                                                                                                               >  1:26506/1‑139 (MQ=255)
                      cgcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTTTTTTTCCCAGTGAGACGGGCaa                                                                                                               >  1:590659/1‑139 (MQ=255)
                      cgcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCaa                                                                                                               <  2:435692/139‑1 (MQ=255)
                      cgcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGC‑TTGC‑GTATTGGGG‑GCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGcga                                                                                                               >  2:337229/1‑137 (MQ=255)
                      cgcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGGGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGC‑TTTC‑GCGGTGGGCAACCAGGG‑GGTGTTTCTTTTTCCCAACGAGGCggggca                                                                                                               >  2:363082/1‑136 (MQ=255)
                        cgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGG‑CCCAGGGGGGGTTTTCTTTTTACCCGTGAGACGGGcca                                                                                                               >  1:333475/1‑135 (MQ=255)
                        cgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACa                                                                                                             <  2:347147/139‑1 (MQ=255)
                        cgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GGATTTGGC‑GCCAGGGGGGGTTTTTTTTTTACCAGGGGGACGGGGaa                                                                                                               >  2:461647/1‑137 (MQ=255)
                             cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑Ga                                                                                                       <  1:309490/139‑1 (MQ=255)
                                       ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCa                                                                                             <  1:312132/139‑1 (MQ=255)
                                       ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCa                                                                                             <  2:371773/139‑1 (MQ=255)
                                            ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GGCTTTGGC‑GCCAGGGGGGGTTTTTTTTTTACCCGGGAGAAGGGGAACAGCTGGTTTGCCCTTcccccc                                                                                         >  2:67140/1‑137 (MQ=255)
                                              ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTTC‑GGATTGGGC‑ACCAGGGGGGGTTTTTTTTTTACCCGTGCGACGGGGCACAGCT‑GATTGCCCTTCACCCCCTg                                                                                      >  2:299127/1‑139 (MQ=255)
                                              ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCt                                                                                       <  1:299259/138‑1 (MQ=255)
                                                    gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTCTTGC‑GTTTTGGTG‑GCCGGGGT‑GTTTTTTTTTTTTCCAGAGGGACGGGGCACAGCT‑GATTTCCCCTCCAcgccccgccccg                                                                                >  2:448370/1‑137 (MQ=255)
                                                     ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTTC‑GTATTGGGG‑ACCAGGGTGGGTTTTTTTTTCACCAGTGGGACGGGGAAAAGAC‑GGTTTCCCTTCAACCCCCGGCCCTgc                                                                               >  2:244799/1‑138 (MQ=255)
                                                     ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTga                                                                               <  1:337229/139‑1 (MQ=255)
                                                            tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GGATTGGGG‑ACCCGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGGAAAAGCT‑GATTGCCCCTCACCCCCTGGCCCTGagagaggt                                                                        >  2:63240/1‑137 (MQ=255)
                                                                gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GGATTTGGG‑ACCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCT‑GATTTCCCTTTACCGCCCGGCCCCGTGAGAG‑TTgcag                                                                   >  2:25734/1‑139 (MQ=255)
                                                                   ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGG‑AACAGGGGGGGTTTTTTTTTTACCAGGGGGAAGGGCGACAGCT‑GATTGGCCTTCACCGCCTGGCCCCGAGAGAG‑TTgcagca                                                                 >  2:448616/1‑138 (MQ=255)
                                                                    cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAg                                                               <  1:420884/139‑1 (MQ=255)
                                                                    cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑CTGC‑GTATTGGGC‑TCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCACCAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAg                                                               <  1:201670/139‑1 (MQ=255)
                                                                        tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTTC‑GTATTGGGG‑ACCAGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGCAACAGCT‑CATTGCCCTTTACCGCCCGGCCCCGGGAGAG‑TTTCAGCAAGAGGt                                                           >  1:182258/1‑139 (MQ=255)
                                                                          cATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTc                                                          <  1:448616/138‑1 (MQ=255)
                                                                            ttAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTTC‑GGATTGGGC‑ACCAGGGGGGTTTTTCTTTTTACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAACGGTCCAc                                                       >  1:231271/1‑139 (MQ=255)
                                                                              aaTGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACg                                                      <  1:180138/138‑1 (MQ=255)
                                                                                  aaTCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGt                                                 <  1:225774/139‑1 (MQ=255)
                                                                                   aTCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGtt                                                <  1:335872/139‑1 (MQ=255)
                                                                                    tCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTTT‑GTATTGGGG‑ACCAGGGGGGTTTTTTTTTTTACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGttt                                               >  1:519385/1‑139 (MQ=255)
                                                                                    tCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTTGGG‑ACCAGGGGGGGTTTTTTTTTCACCAGTGGGACGGGCAACAGCT‑GGTTGCCCTTCACCCCCTGGCCCTGAGGGAG‑GTGCAGCAAGCGGTCCCCGCTGGgtt                                               >  2:180138/1‑139 (MQ=255)
                                                                                    tCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTTGGC‑ACCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGGAACAGCT‑GATTTCCCTTCCCCCCCTGGCCCCGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGgtt                                               >  1:299506/1‑139 (MQ=255)
                                                                                    tCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTCT‑GGAATGGGC‑ACCAGGGGGGTTTTTCTTTTCACCACGGAGACGGGCAACAACG‑GAATGCCCTTTACCCCCTCGCCCTGAGAGAG‑TTTCAGCAAGCGGGACACGCCGGctt                                               >  1:215543/1‑139 (MQ=255)
                                                                                       gCCAACGCGCGGGGAGAGGC‑GGT‑TTTC‑GGATTTGGG‑ACCAGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGAACAGCC‑GATTGCCCTTCTCCCCCCGGCCCCGAGAGAG‑GTGCAGCAAACGGGCCACGCCGGGTTgtc                                            >  2:361978/1‑137 (MQ=255)
                                                                                       gCCAACGCGCGGGGAGAGGC‑GGT‑TTTC‑GGATTGGGC‑ACCAGGGTGGGTTTTCTTTTCACCCGTGAGAAGGGCAAAAACC‑GATTGCCCCTCACCCCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGGTCACGCCGGTTTGcc                                            >  1:454822/1‑139 (MQ=255)
                                                                                         cAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GGATTGGGC‑ACCAGGGGGGTTTTTTTTTTTACCCAAGAGAAGGGCAACAGGA‑GAATGCCCCCCCCCCCCCGGCCCTGAAAGAGCTTGC‑GCGAACGGTACCCGCGGGTTTTcccc                                          >  2:423772/1‑139 (MQ=255)
                                                                                                  gggAGAGGC‑GGT‑TTGC‑GTATTGGGG‑ACCCGGGGGGTTTTTTTTTTTACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCg                                  >  2:78537/1‑138 (MQ=255)
                                                                                                       agGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATc                            <  2:309339/139‑1 (MQ=255)
                                                                                                                  gC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGa                    <  2:454822/138‑1 (MQ=255)
                                                                                                                   c‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtg                  <  2:69691/139‑1 (MQ=255)
                                                                                                                        ttGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg               <  2:550347/138‑1 (MQ=255)
                                                                                                                               aCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCt                                                                                       <  1:320425/61‑1 (MQ=255)
                                                                                                                               aCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCt                                                                                       >  2:320425/1‑61 (MQ=255)
                                                                                                                               aCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg        <  2:215543/139‑1 (MQ=255)
                                                                                                                                ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg       <  1:551418/139‑1 (MQ=255)
                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAg                                                               >  1:302918/1‑82 (MQ=255)
                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAg                                                               <  2:302918/82‑1 (MQ=255)
                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc      <  1:63240/139‑1 (MQ=255)
                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc      <  1:4141/139‑1 (MQ=255)
                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc      <  1:337080/139‑1 (MQ=255)
                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc      >  2:592178/1‑139 (MQ=255)
                                                                                                                                 cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc      >  2:313548/1‑139 (MQ=255)
                                                                                                                                   gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg    <  2:26506/139‑1 (MQ=255)
                                                                                                                                   gggTGGTTTTTCGTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg    <  1:479429/139‑1 (MQ=255)
                                                                                                                                    ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCt                                                                                 <  2:18592/62‑1 (MQ=255)
                                                                                                                                    ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCt                                                                                 >  1:18592/1‑62 (MQ=255)
                                                                                                                                      tggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa  >  1:12246/1‑138 (MQ=255)
                                                                                                                                      |                                                                                                                                           
GCGCAGATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGT‑TTGC‑GTATTGGGC‑ACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCT‑GATTGCCCTTCACCGCCTGGCCCTGAGAGAG‑TTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGA  >  NZ_CP009273/360629‑360896

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