Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273277,9850TG40.0% 56.8 / 109.9 75intergenic (‑291/‑4)yagA/yagEIS481‑like element ISEc19 family transposase/2‑keto‑3‑deoxygluconate aldolase
Reads supporting (aligned to +/- strand):  ref base T (31/14);  new base G (19/11);  total (50/25)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACC  >  NZ_CP009273/277899‑278075
                                                                                      |                                                                                           
aCCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCa                                                                                         <  2:644884/91‑1 (MQ=255)
aCCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCa                                                                                         <  2:681864/91‑1 (MQ=255)
     tcCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCg                                                                                    <  1:559343/91‑1 (MQ=255)
          gTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTc                                                                               <  2:425873/91‑1 (MQ=255)
              tGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCcgcg                                                                           >  2:792446/1‑91 (MQ=255)
              tGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCcgcg                                                                           >  2:715074/1‑91 (MQ=255)
              tGCGTCATACTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCACGGAGAGCTCACGCCCAAGTCcgcg                                                                           >  2:406403/1‑91 (MQ=255)
               gCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt                                                                          >  2:772408/1‑91 (MQ=255)
               gCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt                                                                          <  1:530515/91‑1 (MQ=255)
               gCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCgt                                                                          >  2:771998/1‑91 (MQ=255)
               gCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCgt                                                                          >  1:388921/1‑91 (MQ=255)
                cGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTCTAATGTAGAACAAAATGCTTACGCAAGGAGATCTCATGCCGCAGTCCGCgtt                                                                         >  1:788675/1‑91 (MQ=255)
                    aTAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTc                                                                     <  1:339389/91‑1 (MQ=255)
                       gTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACg                                                                 <  2:13992/91‑1 (MQ=255)
                       gTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACg                                                                 <  2:568317/91‑1 (MQ=255)
                       gTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACg                                                                 <  1:472907/91‑1 (MQ=255)
                           gATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGaa                                                             <  2:424620/91‑1 (MQ=255)
                            aTCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGCTCACGGGAAt                                                            >  2:450206/1‑91 (MQ=255)
                             tCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATc                                                           >  2:203375/1‑91 (MQ=255)
                             tCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATc                                                           <  2:342219/91‑1 (MQ=255)
                              cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCg                                                                                    >  1:723284/1‑66 (MQ=255)
                              cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCg                                                                                    <  2:723284/66‑1 (MQ=255)
                               gTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAACTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATGAt                                                         >  2:314470/1‑91 (MQ=255)
                               gTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCAt                                                         >  1:289476/1‑91 (MQ=255)
                                 taaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt                                                                          >  2:572506/1‑73 (MQ=255)
                                 taaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt                                                                          <  1:572506/73‑1 (MQ=255)
                                 taaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTc                                                       >  2:390838/1‑91 (MQ=255)
                                     taaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTccccc                                                   <  1:429820/91‑1 (MQ=255)
                                       aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAgct                                                                                           <  1:777335/50‑1 (MQ=255)
                                       aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGGAATCATTCCCCCt                                                  >  1:369436/1‑91 (MQ=255)
                                       aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGaa                                                             <  2:253971/79‑1 (MQ=255)
                                       aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGaa                                                             >  1:253971/1‑79 (MQ=255)
                                       aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTg                                                 >  2:674836/1‑91 (MQ=255)
                                       aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTg                                                 >  1:786906/1‑91 (MQ=255)
                                        acgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACgg                                                                >  1:213476/1‑75 (MQ=255)
                                        acgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACgg                                                                <  2:213476/75‑1 (MQ=255)
                                           aacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGtctc                                             <  1:287741/91‑1 (MQ=255)
                                            acgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTcc                                            >  1:325863/1‑91 (MQ=255)
                                                 gTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCAt                                       >  2:478257/1‑91 (MQ=255)
                                                         aTGTAGAACAAAATGATTCAGCGAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccg                               >  1:682201/1‑91 (MQ=255)
                                                         aTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccg                               >  1:682295/1‑91 (MQ=255)
                                                         aTGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccg                               >  2:49570/1‑91 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc                             >  2:846158/1‑91 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc                             >  2:43714/1‑91 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc                             >  2:762257/1‑91 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc                             >  2:371224/1‑91 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc                             >  1:746723/1‑91 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc                             >  2:205331/1‑91 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc                             >  1:800702/1‑91 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACAATTTTTAccgcc                             >  1:759685/1‑91 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCgt                                                                          >  1:728020/1‑47 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCgt                                                                          <  2:728020/47‑1 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATc                                                           <  2:806221/61‑1 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATc                                                           >  1:806221/1‑61 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc                             >  2:44781/1‑91 (MQ=255)
                                                           gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc                             >  2:574500/1‑91 (MQ=255)
                                                            tAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTAccgccg                            <  2:786906/91‑1 (MQ=255)
                                                              gAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc                          >  2:729640/1‑91 (MQ=255)
                                                              gAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc                          >  2:521299/1‑91 (MQ=255)
                                                              gAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc                          >  1:380009/1‑91 (MQ=255)
                                                              gAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc                          >  2:516175/1‑91 (MQ=255)
                                                              gAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTc                                                                               <  2:681551/39‑1 (MQ=39)
                                                              gAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTc                                                                               >  1:681551/1‑39 (MQ=39)
                                                              gAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc                          >  2:56380/1‑91 (MQ=255)
                                                              gAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc                          >  2:195356/1‑91 (MQ=255)
                                                              gAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc                          >  1:855802/1‑91 (MQ=255)
                                                                   aaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCa                     >  1:159906/1‑91 (MQ=255)
                                                                   aaaTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTcca                                           >  1:621232/1‑69 (MQ=255)
                                                                   aaaTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTcca                                           <  2:621232/69‑1 (MQ=255)
                                                                   aaaTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCa                     >  1:863815/1‑91 (MQ=255)
                                                                   aaaTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCa                     >  1:718868/1‑91 (MQ=255)
                                                                   aaaTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCa                     >  2:647950/1‑91 (MQ=255)
                                                                     aTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCGCCATTTTTACCGCCGACGGCCAGc                   >  1:320230/1‑91 (MQ=255)
                                                                        aTTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCg                >  1:757476/1‑91 (MQ=255)
                                                                         ttCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGa               <  1:478257/91‑1 (MQ=255)
                                                                           cAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATa             <  1:792446/91‑1 (MQ=255)
                                                                           cAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATa             <  1:583448/91‑1 (MQ=255)
                                                                                aGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCg        >  1:648494/1‑91 (MQ=255)
                                                                                     aTCTAATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCAc   <  2:453255/91‑1 (MQ=255)
                                                                                      tctcATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCAcc  >  1:246941/1‑91 (MQ=255)
                                                                                      |                                                                                           
ACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACC  >  NZ_CP009273/277899‑278075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: