Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,594,1610CT52.9% 71.0 / 196.9 137S62F (TCT→TTT) bglJDNA‑binding transcriptional activator BglJ
Reads supporting (aligned to +/- strand):  ref base C (10/54);  new base T (72/0);  total (83/54)
Fisher's exact test for biased strand distribution p-value = 4.54e-28
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-27
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 147 total aligned reads displayed.

AACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGT  >  NZ_CP009273/4594076‑4594248
                                                                                     |                                                                                       
aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttttttt                                                                                     >  1:2833742/1‑90 (MQ=255)
  cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                   <  2:1514577/90‑1 (MQ=255)
    cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCAttttttttttttttt                                                                                 >  2:1203731/1‑87 (MQ=255)
      tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                               >  2:405885/1‑90 (MQ=255)
         aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATg                                                                            <  1:1055382/90‑1 (MQ=255)
           gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgag                                                                          <  2:1662687/90‑1 (MQ=255)
            cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATgaga                                                                         >  2:1059453/1‑90 (MQ=255)
            cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATgaga                                                                         >  1:28328/1‑90 (MQ=255)
            cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATgaga                                                                         >  1:2358627/1‑90 (MQ=255)
            cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATggga                                                                         >  1:358058/1‑90 (MQ=255)
            cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATgaga                                                                         >  2:672212/1‑90 (MQ=255)
            cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgaga                                                                         >  2:487836/1‑90 (MQ=255)
            cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTgaga                                                                         >  2:2482435/1‑90 (MQ=255)
            cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATgaga                                                                         >  1:348701/1‑90 (MQ=255)
             tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGaa                                                                        >  1:2676441/1‑90 (MQ=255)
             tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGaa                                                                        >  2:1833885/1‑90 (MQ=255)
              ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAg                                                                       >  1:933670/1‑90 (MQ=255)
              ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATgaggag                                                                       >  1:1342225/1‑90 (MQ=255)
               cATACCTATACCAATCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATgaggagg                                                                      >  2:2744916/1‑89 (MQ=255)
                aTACCTATCCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTg                                                                     <  2:2842418/90‑1 (MQ=255)
                aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTg                                                                     <  1:1106908/90‑1 (MQ=255)
                   ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgc                                                                  <  2:1729613/90‑1 (MQ=255)
                    cTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcg                                                                 <  1:825228/90‑1 (MQ=255)
                     tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGTGAgcgc                                                                >  1:21678/1‑90 (MQ=255)
                     tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTTAGAAGTGAgggc                                                                >  1:2906145/1‑90 (MQ=255)
                      ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAgcgga                                                               >  1:1825788/1‑88 (MQ=255)
                      ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCa                                                               >  1:2161195/1‑90 (MQ=255)
                      ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCa                                                               >  1:1166348/1‑90 (MQ=255)
                      ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCa                                                               >  1:2921870/1‑90 (MQ=255)
                      ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGGGAGCGCa                                                               >  2:6517/1‑90 (MQ=255)
                       taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGGGCag                                                              >  2:2654244/1‑90 (MQ=255)
                       taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTCTCTGCCATGAGAAGTGAGCGCag                                                              <  1:2654244/90‑1 (MQ=255)
                        aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga                                                             <  2:2238521/90‑1 (MQ=255)
                          cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                           <  1:1064096/90‑1 (MQ=255)
                          cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                           <  1:487836/90‑1 (MQ=255)
                          cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                           <  1:2125832/90‑1 (MQ=255)
                          cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                           <  2:2968132/90‑1 (MQ=255)
                          cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga                                                           <  1:1570270/90‑1 (MQ=255)
                          cAGCCAGGAGTCATTTCAGGATGCCATGCCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTGAGAAGTGAGCGCagaga                                                           <  2:2669773/90‑1 (MQ=255)
                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGGGGGCACagagag                                                          >  1:2008290/1‑90 (MQ=255)
                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGGGAAGGGAGGGCcgggag                                                          >  2:1676536/1‑90 (MQ=255)
                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCagagag                                                          >  2:442978/1‑90 (MQ=255)
                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCagagag                                                          >  1:245260/1‑90 (MQ=255)
                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCagagag                                                          >  2:1296077/1‑90 (MQ=255)
                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCagagag                                                          >  2:1212055/1‑90 (MQ=255)
                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgag                                                                          <  1:1087802/74‑1 (MQ=255)
                           aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATGAGGAGGGGGAGCcgagag                                                          >  2:2714009/1‑90 (MQ=255)
                           aGTCAGGACTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                <  1:386639/84‑1 (MQ=255)
                           aGTCAGGACTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                >  2:386639/1‑84 (MQ=255)
                            gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAttttttttttttt                                                                                   >  1:853365/1‑63 (MQ=255)
                            gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                   <  2:853365/64‑1 (MQ=255)
                            gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                   <  1:1901307/64‑1 (MQ=255)
                            gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAATGGGCGCAGAGAgg                                                         >  1:1610076/1‑90 (MQ=255)
                            gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGGGAAGTGAGCGCAGAGAgg                                                         >  2:367470/1‑90 (MQ=255)
                            gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGGGCAGAGAgg                                                         >  2:1048896/1‑90 (MQ=255)
                            gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCATTTTCTCCCGTTAGAAGTGAGCGCGGAGAgg                                                         >  2:1720747/1‑90 (MQ=255)
                             tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttct                                                                                        <  1:2961862/58‑1 (MQ=255)
                             tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGGGAAGTGAGGGCAGAGAggg                                                        >  1:269522/1‑90 (MQ=255)
                             tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAGAggg                                                        >  1:171930/1‑90 (MQ=255)
                             tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATgag                                                                          <  2:1089677/72‑1 (MQ=255)
                             tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTTCCTTGAGAAGGGAGGGCAGAGgggg                                                        >  1:2237399/1‑90 (MQ=255)
                              cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                       <  1:3020851/90‑1 (MQ=255)
                              cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                       <  2:1095956/90‑1 (MQ=255)
                              cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa                                                       <  1:282562/90‑1 (MQ=255)
                                ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGTGAAGGGAGCGCCGAGAGGGGtt                                                     >  2:1526846/1‑90 (MQ=255)
                                ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTTCCCTGAGGAGTGAGCGCagggagggggt                                                     >  1:2329957/1‑87 (MQ=255)
                                ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGtt                                                     >  1:1845142/1‑90 (MQ=255)
                                   gTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGATTATc                                                  >  2:350009/1‑90 (MQ=255)
                                     cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                <  1:740937/74‑1 (MQ=255)
                                       tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                              <  2:245260/90‑1 (MQ=255)
                                            ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                         <  2:1845142/90‑1 (MQ=255)
                                            ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                         <  1:449364/90‑1 (MQ=255)
                                            ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                         <  2:2934309/90‑1 (MQ=255)
                                               tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGGGGGGGTATCTTTCctggctga                                      >  1:1034299/1‑90 (MQ=255)
                                               tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCCGAGAGGGGTTATCTTGGctgactga                                      >  2:2028866/1‑90 (MQ=255)
                                               tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCagaga                                                           >  1:949127/1‑69 (MQ=255)
                                                 ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGAGGGGGTTTTTTTCCTGACTGAAc                                    >  1:400950/1‑90 (MQ=255)
                                                  cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                   <  1:747472/90‑1 (MQ=255)
                                                  cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                   <  1:1019327/90‑1 (MQ=255)
                                                  cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                   <  2:2928278/90‑1 (MQ=255)
                                                  cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                   <  2:2971142/90‑1 (MQ=255)
                                                  cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                   <  2:2342210/90‑1 (MQ=255)
                                                    tGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGGGGGGGCAGAGGGGGGGTATTTTTCCTGactggactgg                                 >  1:2601651/1‑90 (MQ=255)
                                                    tGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTgg                                 <  2:47485/90‑1 (MQ=255)
                                                      tCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTCCATGAGAAGGGAGCGCAGAGGGGGATTATCTTGCCTGACTGAAATGGCg                               >  2:762139/1‑90 (MQ=255)
                                                        gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGGGCGCAGAGGGGGGTTTTCTTTCCTGACTGAAATGGCggt                             >  1:581109/1‑88 (MQ=255)
                                                        gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAGGGGGGTTATCTTGCCTGACTGAACTGGCggt                             >  1:152423/1‑88 (MQ=255)
                                                        gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCAGAGAGGGGGTATCTTTCCTGACTGAAATGGCggt                             >  1:2083442/1‑88 (MQ=255)
                                                        gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGGGCGCAGAGAGGGGGTATTTTTCCTGACTGAACTGGCGAt                             >  1:1667647/1‑90 (MQ=255)
                                                        gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTTTTTTGCCTGGCTGGACTGGCggt                             >  2:2295935/1‑88 (MQ=255)
                                                        gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCggt                             >  2:356768/1‑88 (MQ=255)
                                                          ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTa                           <  2:2163727/90‑1 (MQ=255)
                                                             tctcGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGGGGGGTTATTTTGCCTGACTGAACTGGCGGTTAAGt                        >  1:879236/1‑90 (MQ=255)
                                                              ctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGtt                       <  1:2482209/90‑1 (MQ=255)
                                                                   ttGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCcgc                  >  2:1369256/1‑90 (MQ=255)
                                                                   ttGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGGGGGGTTATCTTGCCTGACTGAACTGGGGATTAAGTTTCcgc                  >  2:1948584/1‑90 (MQ=255)
                                                                        gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc             <  2:1601899/90‑1 (MQ=255)
                                                                         cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc            <  1:1493876/90‑1 (MQ=255)
                                                                           gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg          <  1:656784/90‑1 (MQ=255)
                                                                              aTTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt       <  1:140636/90‑1 (MQ=255)
                                                                                     |                                                                                       
AACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGT  >  NZ_CP009273/4594076‑4594248

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