Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0349433,245,0500CA36.8% 46.5 / 58.3 57intergenic (‑72/+121)EO946_RS16785/EO946_RS16790UvrD‑helicase domain‑containing protein/trimeric intracellular cation channel family protein
Reads supporting (aligned to +/- strand):  ref base C (21/15);  new base A (6/15);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 5.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTCCTGCAGCACCTCATCTATTCTTGACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGCTGTTAGTTATCTAAACGAAA  >  NZ_CP034943/3244915‑3245188
                                                                                                                                       |                                                                                                                                          
gCTCCTGCAGCACCTCATCGACTCTTGACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAata                                                                                                                                        <  2:1172229/140‑1 (MQ=255)
      gcagcaCCTCATCTATTCTTGACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGAt                                                                                                                                  >  2:1202987/1‑140 (MQ=255)
            ccTCATCTATTCTTGACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGc                                                                                                                            >  2:664718/1‑140 (MQ=255)
                aTCTATTCTTGACTGCTCGTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGc                                                                                                                        <  1:1304809/140‑1 (MQ=255)
                   tcttCGTGACTGCTCTTCCTTACATTCCTTATCCTGCTGATTCATTAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACGGTAAAATATAATTGGATACAGATAATCAATATATGATTGAAGCTGGCaaa                                                                                                                     <  2:1074809/138‑1 (MQ=255)
                   tATTCTTGACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCaaa                                                                                                                     <  2:1477866/140‑1 (MQ=255)
                    aTTCTTGACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAAt                                                                                                                    >  2:1975834/1‑140 (MQ=255)
                         tGACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAAt                                                                                                               <  2:796245/140‑1 (MQ=255)
                         tGACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAAt                                                                                                               <  2:714640/140‑1 (MQ=255)
                          gACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAAtt                                                                                                              >  2:813752/1‑140 (MQ=255)
                          gACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAAtt                                                                                                              >  1:94421/1‑140 (MQ=255)
                          gACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCaaa                                                                                                                     >  2:880011/1‑133 (MQ=255)
                          gACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCaaa                                                                                                                     <  1:880011/133‑1 (MQ=255)
                          gACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAAtt                                                                                                              >  1:1652498/1‑140 (MQ=255)
                          gACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAAtt                                                                                                              >  1:2052927/1‑140 (MQ=255)
                             ttctctTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTAt                                                                                                           <  2:1515457/138‑1 (MQ=255)
                             tGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATa                                                                                                                   <  1:411310/132‑1 (MQ=255)
                             tGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATa                                                                                                                   >  2:411310/1‑132 (MQ=255)
                                  ttccttccATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAg                                                                                                      <  1:1541158/140‑1 (MQ=255)
                                  ttccttccATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAg                                                                                                      <  1:286027/140‑1 (MQ=255)
                                  ttccttccATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAg                                                                                                      <  1:286063/140‑1 (MQ=255)
                                    ccttccATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATTATTGAAGCTGGCAAATACAATTTATCTTATca                                                                                                    >  2:1103184/1‑140 (MQ=255)
                                    ccttccATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAgaa                                                                                                    >  2:1856115/1‑138 (MQ=255)
                                     cttccATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTGGcac                                                                                                   <  1:1202987/140‑1 (MQ=255)
                                      ttccATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTTTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGcaca                                                                                                  <  2:2004336/140‑1 (MQ=255)
                                      ttccATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGcaca                                                                                                  <  2:350317/140‑1 (MQ=255)
                                      ttccATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGcaca                                                                                                  <  2:350370/140‑1 (MQ=255)
                                      ttccATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGcaca                                                                                                  <  2:1360435/140‑1 (MQ=255)
                                      ttccATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGcaca                                                                                                  <  2:1981773/140‑1 (MQ=255)
                                                     tgctgATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCtcaa                                                                                   >  1:955992/1‑140 (MQ=255)
                                                      gctgATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAAATTATCTTAGCACATTCACCCACTCtcaat                                                                                  >  2:1371794/1‑140 (MQ=255)
                                                       ctgATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTTTAAAATATAATTGGATACATATAATCAAAATATTATTGAATCTGGCAAATAAAAATTATCTTAGCACATTCACCCAATAtcaatc                                                                                 >  1:2060445/1‑140 (MQ=255)
                                                       ctgATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCtcaatc                                                                                 >  2:227084/1‑140 (MQ=255)
                                                       ctgATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCtcaatc                                                                                 >  1:1667/1‑140 (MQ=255)
                                                       ctgATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCtcaatc                                                                                 >  2:888378/1‑140 (MQ=255)
                                                       ctgATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACAGATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCtcaatc                                                                                 >  1:2012177/1‑140 (MQ=255)
                                                       ctgATTAATGAACGGCACCTCCAGAAAATATTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATTATTGAAGCTGGAAAATAAAATTTATCTTAGCACATTCACCCAATCtcaatc                                                                                 <  1:2047886/140‑1 (MQ=255)
                                                            tCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGt                                                                            <  1:1324647/140‑1 (MQ=255)
                                                            tCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATAAATGATTGAAGCTGGCAAATACAATTTATCTTATCACATTCACCCAATCTCAAACAAGGt                                                                            >  1:956484/1‑140 (MQ=255)
                                                            tCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGt                                                                            <  2:1782075/140‑1 (MQ=255)
                                                            tCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGt                                                                            <  1:1975834/140‑1 (MQ=255)
                                                             cATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGtt                                                                           <  1:664718/140‑1 (MQ=255)
                                                                 aaCGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTggg                                                                       <  2:1972928/140‑1 (MQ=255)
                                                                 aaCGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTggg                                                                       <  2:874895/140‑1 (MQ=255)
                                                                                                       aaTACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATc                                 <  2:1738105/140‑1 (MQ=255)
                                                                                                                  atatAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGc                      <  2:94421/140‑1 (MQ=255)
                                                                                                                           gATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATAGGTTTTTTTGCTgttagtta             >  1:453253/1‑140 (MQ=255)
                                                                                                                             tACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGCt                     <  2:280365/130‑1 (MQ=255)
                                                                                                                             tACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGCt                     >  1:280365/1‑130 (MQ=255)
                                                                                                                             tACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGCTGTTAGTTATc           >  1:1209851/1‑140 (MQ=255)
                                                                                                                             tACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGc                      <  2:1699480/129‑1 (MQ=255)
                                                                                                                             tACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGc                      >  1:1699480/1‑129 (MQ=255)
                                                                                                                             tACATATAATAAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGCTGTTAGTTATc           >  2:1828244/1‑140 (MQ=255)
                                                                                                                                atatAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCACTCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAATTCGGTTTTTTTGCTGTTAGTTATCTaa        >  1:575448/1‑140 (MQ=255)
                                                                                                                                     aTCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGCTGTTAGTTATCTAAACGaa   <  2:1667/140‑1 (MQ=255)
                                                                                                                                      tCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGCTGTTAGTTATCTAAACGaaa  >  2:1625933/1‑140 (MQ=255)
                                                                                                                                      tCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGCTGGTATTTATCTAAACGaaa  >  2:1151432/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                          
GCTCCTGCAGCACCTCATCTATTCTTGACTGCTCTTCCTTCCATTCCTTATCCTGCTGATTCATGAACGGCACCTCCAGAAAATAGTTGACAAACACAGATGAAATACTGTAAAATATAATTGGATACATATAATCAATATATGATTGAAGCTGGCAAATACAATTTATCTTAGCACATTCACCCAATCTCAATCAAGGTTTGGGTTATTTATATGCAAAAAAACCGATTTCAATTTGAAATCGGTTTTTTTGCTGTTAGTTATCTAAACGAAA  >  NZ_CP034943/3244915‑3245188

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: