Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,594,1610CT54.9% 38.7 / 146.8 111S62F (TCT→TTT) bglJDNA‑binding transcriptional activator BglJ
Reads supporting (aligned to +/- strand):  ref base C (10/40);  new base T (61/0);  total (71/40)
Fisher's exact test for biased strand distribution p-value = 4.04e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-17
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 124 total aligned reads displayed.

TAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTT  >  NZ_CP009273/4594075‑4594249
                                                                                      |                                                                                         
tAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttt                                                                                        >  1:253829/1‑90 (MQ=255)
  aCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttcttttt                                                                                      <  2:1386746/90‑1 (MQ=255)
   cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                     <  1:1027888/90‑1 (MQ=255)
   cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                     <  2:218721/90‑1 (MQ=255)
     cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                   <  1:1080102/90‑1 (MQ=255)
     cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                   <  2:1268863/90‑1 (MQ=255)
     cTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                   <  1:400397/90‑1 (MQ=255)
      ttATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAtttttttttttttttc                                                                                  >  1:134636/1‑86 (MQ=255)
       tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGcc                                                                                 <  1:238851/90‑1 (MQ=255)
            gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑Tgag                                                                           <  2:1013078/90‑1 (MQ=255)
             cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTGTTTTTTTTTTTTTTCCA‑Tgaga                                                                          >  2:1479885/1‑90 (MQ=255)
             cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCC‑Tgaga                                                                          >  1:123694/1‑90 (MQ=255)
             cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCA‑Tgaga                                                                          >  2:665685/1‑90 (MQ=255)
             cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑Tgaga                                                                          >  2:212068/1‑90 (MQ=255)
             cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTTCA‑Tgaga                                                                          >  1:1227738/1‑90 (MQ=255)
             cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTGGCA‑Tgagg                                                                          >  2:1049061/1‑89 (MQ=255)
               ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCA‑Tgaggag                                                                        >  2:967184/1‑90 (MQ=255)
               ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCC‑TGAGaaa                                                                        >  1:1516716/1‑89 (MQ=255)
               ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCA‑TGAGAAg                                                                        >  2:115266/1‑90 (MQ=255)
               ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAg                                                                        <  2:1342515/90‑1 (MQ=255)
                cATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCA‑TGAGAAgg                                                                       >  2:746338/1‑89 (MQ=255)
                 aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTGCTTTTTTTTTTTTTTTGCA‑TGAGAAggg                                                                      >  2:1068737/1‑88 (MQ=255)
                 aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTg                                                                      <  2:880563/90‑1 (MQ=255)
                 aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTg                                                                      <  2:851506/90‑1 (MQ=255)
                  tACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGa                                                                     <  2:1207903/90‑1 (MQ=255)
                      tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCAttttttttttttt                                                                                     >  1:655210/1‑70 (MQ=255)
                      tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttt                                                                                     >  2:276040/1‑70 (MQ=255)
                      tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                     <  1:276040/71‑1 (MQ=255)
                      tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                     <  2:655210/71‑1 (MQ=255)
                      tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCA‑TGAGAAGGGAgggc                                                                 >  2:984134/1‑90 (MQ=255)
                      tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAgcgc                                                                 <  1:1284134/90‑1 (MQ=255)
                      tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAgcgc                                                                 <  1:1555368/90‑1 (MQ=255)
                      tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAgcgc                                                                 <  2:1390412/90‑1 (MQ=255)
                       ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGCCA‑TTAGAAGTGAGGGCa                                                                >  2:734957/1‑90 (MQ=255)
                       ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCC‑TGAGAAGGGAGCGCa                                                                >  1:779015/1‑90 (MQ=255)
                       ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCC‑TTTGAAGGGAGGGCa                                                                >  2:1297038/1‑90 (MQ=255)
                       ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCA‑TGAGAAGTGAgcggg                                                                >  1:919917/1‑88 (MQ=255)
                       ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCA‑TGAGAAGTGAgcgga                                                                >  1:329244/1‑88 (MQ=255)
                         aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCaga                                                              <  2:55589/90‑1 (MQ=255)
                          ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCagag                                                             <  2:5395/90‑1 (MQ=255)
                           cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCA‑TGGGAAGTGAGCGCagaga                                                            >  2:1063224/1‑90 (MQ=255)
                           cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCagaga                                                            <  1:1362112/90‑1 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttctttttc                                                                                     <  2:766065/65‑1 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATtttttct                                                                                          <  2:1581296/60‑1 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCC‑TGAGAAGTGAGCGCcgagag                                                           >  1:98673/1‑90 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCC‑TGAGAAGGGAGCGCagagag                                                           >  1:105703/1‑90 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCA‑TGAGAAGAGAGCGCagagag                                                           >  1:625039/1‑90 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCC‑TGAGAAGGGAGGGCaggggg                                                           >  2:1556651/1‑86 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCA‑TGAGAAGTGAGGGCagagag                                                           >  2:159414/1‑90 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCA‑TGAGAAGTGAGCGCagagag                                                           >  2:107114/1‑90 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCA‑TGAGAAGGGGGGGCagagag                                                           >  2:1239084/1‑90 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTg                                                                                   >  1:1001078/1‑67 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTg                                                                                   <  2:1001078/67‑1 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAgcg                                                                  >  2:160101/1‑83 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAgcg                                                                  <  1:160101/83‑1 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCagagag                                                           >  2:91202/1‑90 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTGCCG‑TTAGAAGGGTGCGCagagag                                                           >  1:708153/1‑90 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCAttttttttttttt                                                                                     >  1:766065/1‑64 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCC‑TGAGAAGTGGGGGCagagag                                                           >  1:1282556/1‑90 (MQ=255)
                            aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTCCCACTGA‑AACGGAGCGCagagcg                                                           >  2:766119/1‑88 (MQ=255)
                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCC‑TTAGAAGTGGGCGCAGAGAgg                                                          >  1:1066772/1‑90 (MQ=255)
                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCC‑TGAGAAGTGAGGGCAGAGAgg                                                          >  1:1512141/1‑90 (MQ=255)
                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCA‑TGAGAAGTGAGCGCAGAGAgg                                                          >  1:1436573/1‑90 (MQ=255)
                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCA‑TGAGAAGGGAGCGCAGAGAgg                                                          >  1:630461/1‑90 (MQ=255)
                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCA‑TGAGAAGTGAGCGCagaga                                                            >  1:65129/1‑88 (MQ=255)
                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCA‑TGAGAAGTGAGCGCAGAGAgg                                                          >  2:18385/1‑90 (MQ=255)
                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCagaga                                                            <  2:65129/88‑1 (MQ=255)
                             gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCCTTTTTTTTTTTTTTTCCC‑TGAGAAGAGAGCGCAGAGAgg                                                          >  2:1169649/1‑90 (MQ=255)
                              tCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTCTCCCC‑TGAGAAGTGAGCGCAGAGAggg                                                         >  2:1294405/1‑90 (MQ=255)
                               cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCA‑TGAGAAGGGAGCGCAGAGAgggg                                                        >  2:861908/1‑89 (MQ=255)
                               cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGa                                                        <  1:354608/90‑1 (MQ=255)
                               cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGa                                                        <  1:749916/90‑1 (MQ=255)
                               cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGa                                                        <  1:1625038/90‑1 (MQ=255)
                                 ggAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCC‑TGAGAAGGGAGGGCAGAGGGGGGtt                                                      >  1:1367401/1‑90 (MQ=255)
                                      cATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCACTTTTTTTTTTTTTTCCC‑TGAGAAGGGAGCGCAGAGAGGGGTTATCtt                                                 >  1:439702/1‑90 (MQ=255)
                                        tttCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                               <  1:382791/90‑1 (MQ=255)
                                             ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgac                                          >  2:977750/1‑90 (MQ=255)
                                               aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCA‑TGGGAAGGGAGCGCAGAGGGGGGTTATCTTGCctgactg                                        >  1:473517/1‑90 (MQ=255)
                                                tGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCA‑TGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCctgactga                                       >  2:1396335/1‑90 (MQ=255)
                                                tGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCA‑TGAGAAGGGTGGGCAGGGGGGGGGTATTTTGTctgtctga                                       >  1:720361/1‑90 (MQ=255)
                                                   cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                    <  1:91202/90‑1 (MQ=255)
                                                   cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt                                    <  1:1294405/90‑1 (MQ=255)
                                                    aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCA‑TGAGAAGGGAGGGCAGAGAGGGGTTATCTTGCCTGACTTAACTg                                   >  2:43813/1‑90 (MQ=255)
                                                    aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTg                                   >  2:536579/1‑90 (MQ=255)
                                                     tGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCC‑TTAGAAGGGGGGGCGGAGAGGGGGTTTCTTTCCTGACTGAAATgg                                  >  1:921193/1‑90 (MQ=255)
                                                      gTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCC‑TGAGAAGGGAGGGCAGAGGGGGGTTTTCTTTCCTTGCTGAAATGGc                                 >  1:1541905/1‑90 (MQ=255)
                                                         gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCC‑TTAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAAATGGCGAt                              >  2:1146531/1‑90 (MQ=255)
                                                         gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCA‑TGAGAAGGGAGGGCAGAGAGGGGTTATCTTGCCTGACTGAAATGGCggt                              >  2:1113576/1‑88 (MQ=255)
                                                           ggATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCA‑TGAGAAGGGAGGGCAGAGAGGGGTTTTCTTGCCTGACTGAACTGGCGGTTa                            >  2:1236100/1‑90 (MQ=255)
                                                           ggATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTa                            <  2:625039/90‑1 (MQ=255)
                                                               ctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGtt                        <  1:422267/90‑1 (MQ=255)
                                                                tcGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCC‑TGAGAAGGGAGCGCAGGGAGGGGGTATTTTTCCTGTCTGAAATGGGGGTTTAGttt                       >  1:176460/1‑90 (MQ=255)
                                                                tcGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCA‑TGAGAAGGGAGCGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGCGATTAAgctg                       >  2:24610/1‑87 (MQ=255)
                                                                      gcggcgGGCATTTTTTTTTTTTTTTCCA‑TGAGAAGGGAGCGCAGAGGGGGGTTATTTTTCCTGACTGAACTGGCGATTAAGTTTTcgcgg                 >  1:996702/1‑89 (MQ=255)
                                                                        ggcggTCATTTTTTCTTTTTTTTCCA‑TGAGAAGTGAGCGCCGAGGGGGGTTATTTTTTCTTTCTTAACTGGCGGTTAAAGTTCCCCGTAc               >  1:1421428/1‑90 (MQ=255)
                                                                         gcggTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc              <  2:1180691/90‑1 (MQ=255)
                                                                          cggTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc             <  2:98673/90‑1 (MQ=255)
                                                                           ggTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCg            <  2:1158931/90‑1 (MQ=255)
                                                                               aTTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGt        <  1:1396335/90‑1 (MQ=255)
                                                                                     tctttttctGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAATAAGTTTCCGCGTACCCGGCGTTTAGtt  >  1:1428332/1‑90 (MQ=255)
                                                                                      |                                                                                         
TAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCA‑TGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTT  >  NZ_CP009273/4594075‑4594249

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