Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,594,1610CT59.5% ‑3.1 / 97.7 84S62F (TCT→TTT) bglJDNA‑binding transcriptional activator BglJ
Reads supporting (aligned to +/- strand):  ref base C (4/30);  new base T (50/0);  total (54/30)
Fisher's exact test for biased strand distribution p-value = 8.57e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-14
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTA  >  NZ_CP009273/4594077‑4594246
                                                                                    |                                                                                     
aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt                                                                                  <  1:683479/90‑1 (MQ=255)
 cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                 <  1:596468/90‑1 (MQ=255)
 cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                 <  1:68232/90‑1 (MQ=255)
 cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                 <  1:1241298/90‑1 (MQ=255)
 cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                 <  2:414804/90‑1 (MQ=255)
 cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                 <  1:618713/90‑1 (MQ=255)
   cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                               <  1:1394081/90‑1 (MQ=255)
     tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAttttttttttttt                                                                                 >  1:740772/1‑85 (MQ=255)
     tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                 <  2:740772/86‑1 (MQ=255)
           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATgaga                                                                       >  2:1251260/1‑90 (MQ=255)
           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATgaga                                                                       >  1:612949/1‑90 (MQ=255)
           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGGCTTTTTTTTTTTTTTGCCAttgga                                                                       >  1:1251641/1‑86 (MQ=255)
           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATgaga                                                                       >  2:448534/1‑90 (MQ=255)
           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATgaga                                                                       >  1:1145314/1‑90 (MQ=255)
           cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACCTTTTTTTTTTTTTTTCCCTgaga                                                                       >  2:673952/1‑90 (MQ=255)
               aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCGTTTTTTTTTTTTTTTCCTTGTGAAggg                                                                   >  1:951944/1‑88 (MQ=255)
                  ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                <  2:846369/90‑1 (MQ=255)
                    tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCCTTAGAAATGGgcgc                                                              >  1:1598760/1‑90 (MQ=255)
                    tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAgcgc                                                              >  1:624293/1‑90 (MQ=255)
                    tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGAgcgc                                                              >  2:972129/1‑90 (MQ=255)
                    tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAgcgc                                                              >  2:374372/1‑90 (MQ=255)
                    tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                              <  1:1244114/90‑1 (MQ=255)
                    tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTGCCCTGAGAAGGGGgcgc                                                              >  1:1450802/1‑90 (MQ=255)
                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCa                                                             >  2:18267/1‑90 (MQ=255)
                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATTAGAAGTGAgggcc                                                             >  1:753012/1‑89 (MQ=255)
                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGGCATTAGGAGTGAgagcc                                                             >  2:801171/1‑89 (MQ=255)
                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCa                                                             >  1:1526392/1‑90 (MQ=255)
                     ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCa                                                             >  1:76217/1‑90 (MQ=255)
                      taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                            <  1:1543748/90‑1 (MQ=255)
                      taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCag                                                            <  1:256462/90‑1 (MQ=255)
                         cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                         <  2:1287484/90‑1 (MQ=255)
                         cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                         <  2:1567800/90‑1 (MQ=255)
                         cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                         <  2:612949/90‑1 (MQ=255)
                         cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                         <  2:972397/90‑1 (MQ=255)
                         cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                         <  1:1185382/90‑1 (MQ=255)
                         cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCTTTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCagggg                                                         >  1:744604/1‑87 (MQ=255)
                          cgTCAGGAGTCATTTAAGGATGGACCGCCGCGGATGTGGTTTTCGGACGGCATTTTTGCCTTGTATGGCATGAGGAGTGACCGCcgagag                                                        <  1:1166846/89‑1 (MQ=255)
                          aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCagggag                                                        >  2:751721/1‑90 (MQ=255)
                          aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCagagag                                                        >  1:5114/1‑90 (MQ=255)
                          aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCagagag                                                        >  2:246194/1‑90 (MQ=255)
                           gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAgg                                                       >  1:1224496/1‑90 (MQ=255)
                           gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTGCATGGGAAGAGAGCGCAGGGAgg                                                       >  1:577460/1‑90 (MQ=255)
                           gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAgggg                                                       >  1:212001/1‑90 (MQ=255)
                           gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGGGGGCAGAGAgg                                                       >  2:256462/1‑90 (MQ=255)
                           gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAgg                                                       >  2:1131318/1‑90 (MQ=255)
                           gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATGGGAAGTGGGGGCAGAGAgg                                                       >  2:1281831/1‑90 (MQ=255)
                            tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCACTAGAAGTGAGCGGAGAGAggg                                                      >  1:372116/1‑90 (MQ=255)
                            tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCAGAGAggg                                                      >  1:1201216/1‑90 (MQ=255)
                            tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCAGAGAggg                                                      >  1:1008529/1‑90 (MQ=255)
                            tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCCTGAGAAGAGAGCGCCGAGAggg                                                      >  2:1069529/1‑90 (MQ=255)
                            tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGgggg                                                      >  1:159845/1‑90 (MQ=255)
                            tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGAGAAGAGAGCGGAGAGAggg                                                      >  1:14712/1‑90 (MQ=255)
                             cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGa                                                     >  2:720844/1‑90 (MQ=255)
                             cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                     <  2:228938/90‑1 (MQ=255)
                             cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                     <  1:448534/90‑1 (MQ=255)
                               ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGGAGGGGGGGGtt                                                   >  1:1031766/1‑90 (MQ=255)
                               ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGGGAAGTGGGGGCAGAGGGGGGtt                                                   >  2:1543748/1‑90 (MQ=255)
                               ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAgggggt                                                   >  1:1320062/1‑87 (MQ=255)
                                  gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc                                                <  1:673952/90‑1 (MQ=255)
                                   tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCt                                               <  1:298920/90‑1 (MQ=255)
                                    cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCAGAGAGGGGTTATCtt                                              >  2:951175/1‑90 (MQ=255)
                                     aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGAGGGGGTATCTTg                                             >  2:489252/1‑90 (MQ=255)
                                      tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                            <  1:18267/90‑1 (MQ=255)
                                      tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                            <  1:533862/90‑1 (MQ=255)
                                         cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctg                                         <  2:211960/90‑1 (MQ=255)
                                              tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTCCGAAGGGAGGGCCGGGAGGGATTTTTTTGCccgactga                                    >  1:744056/1‑90 (MQ=255)
                                               gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTTCATGAGGAGTGAGGGGGGGGGGGGGTTATCTTTCctgtctgta                                   >  1:576429/1‑88 (MQ=255)
                                               gCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGGGAGGGCAGGGGGGGGTTATTTTTCCTTTCTGaa                                   >  1:112901/1‑90 (MQ=255)
                                                ccgtGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATTAGAAGGGGGGGCAGGGGGGGGTTATCTTTCCTGACTGAAc                                  >  2:131118/4‑90 (MQ=255)
                                                ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAAc                                  >  2:94139/1‑90 (MQ=255)
                                                 cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                 <  1:172989/90‑1 (MQ=255)
                                                    gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCGTTAGAAGTGGGCGCCGGGGGGGGGTATTTTGTCTGACTTAACTGGc                              >  1:1572675/1‑90 (MQ=255)
                                                       gcgGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGGGCGGAGGGGGGTTTTTTTGTCTGACTGAACTGGGGAt                           >  1:325396/1‑90 (MQ=255)
                                                         ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGGAGGGAGCGCAGGGGGGGGTTTTTTTGCCTGACTGAACTGGCGATTa                         >  1:836477/1‑90 (MQ=255)
                                                         ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATTTTGCCTGACTGAAATGGCGATTa                         >  2:1500559/1‑90 (MQ=255)
                                                         ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTa                         >  2:515527/1‑90 (MQ=255)
                                                              tcGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGGGCAGAGAGGGGTTATTTTTCCTGTCTGAACTGGCGGTTTAGttt                    >  1:122920/1‑90 (MQ=255)
                                                              tcGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGATTATCTTTCCTGACTGAACTGGCGGTTTAGttt                    >  2:1243465/1‑90 (MQ=255)
                                                              tcGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGGGAGGGCAGAGAGGGGGTTTTTTTCCTGACTGAACTGGGGATTAAGttt                    >  1:1029327/1‑90 (MQ=255)
                                                               cGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGGGCGCAGAGGGGGGTTATTTTTTCTTACTGAACTGGCGGTTAAGtttt                   >  2:534225/1‑89 (MQ=255)
                                                                       gcggTCATTTTTTCTTTTTCTGCCATGAGGAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc           <  2:1139491/90‑1 (MQ=255)
                                                                        cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc          <  2:1008529/90‑1 (MQ=255)
                                                                        cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc          <  1:640115/90‑1 (MQ=255)
                                                                         ggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg         <  1:1273147/90‑1 (MQ=255)
                                                                         ggTCATTTTTTCTTTTTCTGCCATGAGAAGGGGGGGCAGAGAGGGATTATCTTGCCTGTCTGAACAGGGGGTTAAGTTTCCGCGTACCCg         >  2:664651/1‑90 (MQ=255)
                                                                          gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg        <  1:131118/90‑1 (MQ=255)
                                                                          gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg        <  1:94139/90‑1 (MQ=255)
                                                                           tCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGc       <  2:753012/90‑1 (MQ=255)
                                                                             aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt     <  2:569598/90‑1 (MQ=255)
                                                                             aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt     <  2:1593763/90‑1 (MQ=255)
                                                                              ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGtt    <  1:720844/90‑1 (MQ=255)
                                                                              ttttttCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGtt    <  1:657937/90‑1 (MQ=255)
                                                                               tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGttt   <  1:1500559/90‑1 (MQ=255)
                                                                                ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTa  <  2:577460/90‑1 (MQ=255)
                                                                                    |                                                                                     
ACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTA  >  NZ_CP009273/4594077‑4594246

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑