Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,778,6860TG52.5% 72.3 / 90.9 59intergenic (‑316/+435)BW25113_RS25705/BW25113_RS13825DUF5507 domain‑containing protein/tRNA‑Ile
Reads supporting (aligned to +/- strand):  ref base T (8/20);  new base G (31/0);  total (39/20)
Fisher's exact test for biased strand distribution p-value = 1.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-11
Rejected as consensus: Frequency below/above cutoff threshold.

AATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC  >  NZ_CP009273/2778606‑2778774
                                                                                |                                                                                        
aaTTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc                                                                                 <  1:1167519/90‑1 (MQ=255)
aaTTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc                                                                                 <  2:367641/90‑1 (MQ=255)
aaTTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc                                                                                 <  1:372649/90‑1 (MQ=255)
        cagcAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGtagtag                                                                         <  1:806577/87‑1 (MQ=255)
                  gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTTGGCAGTAGTAGATGTTTAAgg                                                               >  2:248576/1‑90 (MQ=255)
                  gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGGCAGTAGTAGATGttttaag                                                               >  1:1225395/1‑86 (MQ=255)
                  gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTGGGCAGTAGTAGATGTTTAAgg                                                               >  2:54417/1‑90 (MQ=255)
                  gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTCTGGTCAGTAGTAGATGTTTTAgg                                                               >  1:245156/1‑90 (MQ=255)
                     ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGGCAGTAGTAGATGTTTTAGGCgg                                                            >  2:81599/1‑89 (MQ=255)
                     ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTTAGGCgg                                                            >  1:1355716/1‑89 (MQ=255)
                     ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTCTGGGCAGTAGGAGATGTTTAAGGCgg                                                            >  1:808636/1‑89 (MQ=255)
                         caggAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCa                                                        <  1:219080/87‑1 (MQ=255)
                               cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGAt                                                                       <  2:204877/69‑1 (MQ=255)
                               cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGAt                                                                       >  1:204877/1‑69 (MQ=255)
                                      aCGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTTTGGTCAGTAGGAGATGTTTTAGGGGTGGCAGAGACATTTTATc                                           >  1:674935/1‑90 (MQ=255)
                                      aCGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc                                                                                 <  2:219258/52‑1 (MQ=255)
                                      aCGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc                                                                                 >  1:219258/1‑52 (MQ=255)
                                      aCGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACAtt                                                >  1:534955/1‑85 (MQ=255)
                                      aCGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACAtt                                                <  2:534955/85‑1 (MQ=255)
                                      aCGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGGCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATc                                           >  1:1315314/1‑90 (MQ=255)
                                      aCGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTTGTCAGTAGTAGATGTTTAAAGCGGGGGCGAGACAatttctg                                           >  2:219080/1‑86 (MQ=255)
                                      aCGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTCTGGTCAGTAGTAGATGTTTTAGGGGGGGCAGAGACATTTTATc                                           >  2:421908/1‑90 (MQ=255)
                                        gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCCGTAGTAGATGTTTTAGGCGTGGCAGAGACATTTCATCCt                                         >  1:199379/1‑90 (MQ=255)
                                        gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGTGGAAGGAGATGTTTAAAGGCTGGGAGAGACATTTCATCCt                                         >  1:1307535/1‑90 (MQ=255)
                                        gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTGGTCGGTAGTAGGAGTTTTAAGCGTGGCAGAGACATTTTATCCt                                         >  1:1082157/1‑90 (MQ=255)
                                        gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTTTGGTCAGGAGTAGAGGTTTAAAGCGTGGCAGAGACATTTCATCCt                                         >  2:1373797/1‑90 (MQ=255)
                                         cTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCtt                                        <  1:845139/90‑1 (MQ=255)
                                           ccTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTAc                                      <  2:873189/90‑1 (MQ=255)
                                           ccTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTAc                                      <  1:627131/90‑1 (MQ=255)
                                           ccTGTAACGCTCACAAAACTCATCTGCCCGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTAc                                      <  2:1050982/90‑1 (MQ=255)
                                                 aCGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTTATCCCTTCTCTACg                                >  2:578380/1‑90 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa                                                                                >  1:710810/1‑40 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa                                                                                <  2:710810/40‑1 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTTGTCAGTAGTAGATGTTTTAGGCGGGGGAGAGACATTTCATCCTTTATCTTCGGc                              >  2:1075754/1‑90 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGTAGATGTTTAAGGGGGGGCGGAGACATTTTATCCTTTCTCTACGGc                              >  2:1297883/1‑90 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGGCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTTATCCTTACTCTACGGc                              >  2:607393/1‑90 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTTAGGCGTGGGAGAAACATTTCATCCTTACTCTTCGGc                              >  1:670191/1‑90 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTTAGGCGGGGGCGAGACATTTTTTCCTTTATCTACGGc                              >  1:1417633/1‑90 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGc                              >  1:636107/1‑90 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGc                              >  1:309909/1‑90 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTCATCCTTACTCTACGGc                              >  1:1077101/1‑90 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAAGCGGGGGAGAGACATTTTATCCTTTCTCTACGGc                              >  2:627131/1‑90 (MQ=255)
                                                   gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAGTAGTAGATGTTTAAGGCGGGGGAGAGAAATTTCATCCTTTCTCTACGGc                              >  1:1088626/1‑90 (MQ=255)
                                                     tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa                                                                                <  1:1010511/38‑1 (MQ=255)
                                                     tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa                                                                                >  2:1010511/1‑38 (MQ=255)
                                                     tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTTAGGCGGGGGAGAGACATTTTATCCTTTCTCTACGGCAt                            >  1:1297731/1‑90 (MQ=255)
                                                     tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt                            >  1:884147/1‑90 (MQ=255)
                                                     tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt                            <  1:421908/90‑1 (MQ=255)
                                                     tCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGTAGATGTTTAAAGGGGGGGAGAGACATTTTATCCTTTCTCTACGGCAt                            >  1:1273616/1‑90 (MQ=255)
                                                     tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTTGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTTATCCTTTCTCTACGGCAt                            >  2:536519/1‑90 (MQ=255)
                                                     tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTTAAGCGGGGCAGAGACATTTTATCCTTTCTCTACGGCAt                            >  2:167620/1‑90 (MQ=255)
                                                     tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAGGAGTAGATGTTTAAGGCGGGGCAGAGACATTTCATCCTTTCTCTACGGCAt                            >  1:759447/1‑90 (MQ=255)
                                                      cacaAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTCATCCTTACTCTACGGCAtt                           >  1:1339488/1‑90 (MQ=255)
                                                              tCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCtaca                   <  2:757667/90‑1 (MQ=255)
                                                              tCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTTATCCTTTCTCTACGGCATTGTTCtaca                   >  2:423360/1‑90 (MQ=255)
                                                               cATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCtacat                  <  1:536519/90‑1 (MQ=255)
                                                                   tgcctgcGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTacatacat              <  2:884147/90‑1 (MQ=255)
                                                                       tgcGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGt          <  2:144285/90‑1 (MQ=255)
                                                                         cggcggGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTtg        <  1:914861/90‑1 (MQ=255)
                                                                               gtgtTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAc  <  1:167620/90‑1 (MQ=255)
                                                                               gtgtTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAc  <  1:1209486/90‑1 (MQ=255)
                                                                               gtgtTCCGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGGCATTTCATCCTTTCTCTACGGCATTGTTCTACATACATTGGTTGTGGttc  <  2:821712/90‑3 (MQ=255)
                                                                                |                                                                                        
AATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC  >  NZ_CP009273/2778606‑2778774

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