Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,778,6860TG55.6% ‑3.6 / 27.9 27intergenic (‑316/+435)BW25113_RS25705/BW25113_RS13825DUF5507 domain‑containing protein/tRNA‑Ile
Reads supporting (aligned to +/- strand):  ref base T (1/11);  new base G (15/0);  total (16/11)
Fisher's exact test for biased strand distribution p-value = 9.20e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC  >  NZ_CP009273/2778601‑2778774
                                                                                     |                                                                                        
ttAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTc                                                                                      <  2:777304/90‑1 (MQ=255)
              ctcAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGa                                                                        <  2:279772/90‑1 (MQ=255)
                       gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc                                                                                 <  2:602128/72‑1 (MQ=255)
                       gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCa                                                                                <  2:306460/73‑1 (MQ=255)
                       gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTgggc                                                                                 >  1:602128/1‑70 (MQ=255)
                       gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCa                                                                                >  1:306460/1‑73 (MQ=255)
                       gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGGAGGAAAAGttttaag                                                               >  1:530692/1‑86 (MQ=255)
                       gTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAGTAGTAGATGTTTAAgg                                                               >  2:201424/1‑90 (MQ=255)
                         tGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCg                                                             <  2:157436/90‑1 (MQ=255)
                          ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCTTCTGCCTGCGGCGGGGGTTCTGGACAGAAGGAGATGGTTTAGGGGt                                                            >  1:341449/1‑90 (MQ=255)
                          ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTc                                                                                 <  2:249630/69‑1 (MQ=255)
                          ggTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTgggc                                                                                 >  1:249630/1‑67 (MQ=255)
                                    cATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGGAGAGACa                                                  >  1:329904/1‑90 (MQ=255)
                                        gCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTc                                              <  2:652196/90‑1 (MQ=255)
                                             gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGGCAGTAGTAGATGTTTTAAGCGTGGCCGAGACATTTTATCCt                                         >  1:381569/1‑90 (MQ=255)
                                             gCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGGGGTCTTGTCCGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTTATCCt                                         >  2:459240/1‑90 (MQ=255)
                                                        gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTTTGGTCAGTAGGAGATGTTTTAGGCGGGGGCGAGACAATTTATCCTTTCTCTTCGGc                              >  2:251824/1‑90 (MQ=255)
                                                        gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCCCTAATAGATGTTTAAGGCGGGGGAAAGACATTTCATTCTTACTCTACGGc                              >  1:363506/1‑90 (MQ=255)
                                                        gCTCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTTATCCTTACTCTACGGc                              >  1:690739/1‑90 (MQ=255)
                                                          tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt                            <  2:238802/90‑1 (MQ=255)
                                                          tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt                            <  2:329904/90‑1 (MQ=255)
                                                          tCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCAt                            <  2:363506/90‑1 (MQ=255)
                                                          tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTTATCCTTTCTCTACCGCAt                            >  1:642383/1‑90 (MQ=255)
                                                          tCACAAAACTCATCTGCCTGCGGCGGGGGTTCTGGTCAGGAGTAGATGTTTTAGGCGGGGCAGAGACATTTTATCCTTTCTCTTCGGCAt                            >  1:434695/1‑90 (MQ=255)
                                                          tCACAAAACTCATCTGCCTGCGGCGGGGGGTCTGGTCAGTAGTAGATGTTTTAGGCGTGGCAGAGACATTTTATCCTTTCTCTACGGCAt                            >  2:363822/1‑90 (MQ=255)
                                                                       ctgcctgcGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCtacataca               <  1:579610/90‑1 (MQ=255)
                                                                         gcctgcGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGGGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACAtt             >  1:128354/1‑90 (MQ=255)
                                                                                    gtgtTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAc  <  2:100073/90‑1 (MQ=255)
                                                                                     |                                                                                        
TTAGAAATTTATATCTCAATAGCGTTGGTTAATGAGCATAGCCACGCTCCTGTAACGCTCACAAAACTCATCTGCCTGCGGCGGGTGTTCTGGTCAGTAGTAGATGTTTAAGGCGTGGCAGAGACATTTCATCCTTACTCTACGGCATTGTTCTACATACATTGGTTGTGGTAC  >  NZ_CP009273/2778601‑2778774

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