Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,353,1920CA43.8% 21.3 / 10.4 35F106L (TTC→TTA) yrpB2‑nitropropane dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (15/3);  new base A (14/0);  total (32/3)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGATCGCTGCATTCCGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCG‑ACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGC  >  NC_002947/4353046‑4353335
                                                                                                                                                  |                                                                                                                                                
gagaTCGCTGCATTCCGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGATCCTTATTACAGCGAAGTGGGCGCCGATTTc                                                                                                                                                  <  8:139231/147‑1 (MQ=255)
gagaTCGCTGCATTCCGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTc                                                                                                                                                  <  5:103262/147‑1 (MQ=255)
gagaTCGCTGCATTCCGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTc                                                                                                                                                  <  7:164461/147‑1 (MQ=255)
    tCGCTGCATTCCGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTAACgg                                                                                                                                              >  5:335446/1‑147 (MQ=255)
       cTGCATTCCGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGgcgc                                                                                                                                           >  7:272236/1‑147 (MQ=255)
       cTGCATTCCGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGgcgc                                                                                                                                           <  2:324576/147‑1 (MQ=255)
       cTGCATTCCGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGgcg                                                                                                                                            <  4:54721/146‑1 (MQ=255)
           aTTCCGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGAATTAACGGCGCCCAc                                                                                                                                       >  4:84977/1‑147 (MQ=255)
             tCCGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGACCTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGGGCGCCGTTGACGGCGCCCCCgg                                                                                                                                     >  5:4935/1‑146 (MQ=255)
              ccGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGGGGCGCATTAACGGCGCCCACg                                                                                                                                      >  1:302824/1‑145 (MQ=255)
               cgcgCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCCGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTGATTACAGCGAAGTGGGGGCCGCGTTGACGGGGCCCCCCCCg                                                                                                                                   >  5:223350/1‑147 (MQ=255)
                                 cgcCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGGGCCGATTTCACGGCGCCCCCGCCGGTGTCCAACCGTGCGCCt                                                                                                                 >  6:92936/1‑147 (MQ=255)
                                 cgcCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGgcg                                                                                                                                            >  2:287760/1‑120 (MQ=255)
                                 cgcCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGgcg                                                                                                                                            <  1:287760/120‑1 (MQ=255)
                                    ccGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGAATTAAAGGCGCCCAAGCAGGTGTACAACCGTGAGCCttt                                                                                                               >  4:293290/1‑146 (MQ=255)
                                      gCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTAACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTcg‑a                                                                                                           >  8:183045/1‑147 (MQ=255)
                                       ctgaacctgaacTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTAACCGCGCCCACGCCGGTGTCCAACCCTGCGCCTTTcg‑a                                                                                                           >  3:138508/1‑146 (MQ=255)
                                       ctgaacctgaacTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCAATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGGGCCCATTTAACCGCGCCCCCGCCGGTGTGCAACCCTGCGCCTTTag‑ac                                                                                                          >  3:262745/1‑147 (MQ=255)
                                              tgaacTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTAACGGCGCCCAAGCCGGTGTCCAACCGTGCGCCTTTGG‑ACGAGCaga                                                                                                   >  3:225929/1‑147 (MQ=255)
                                                       tctGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCCCGCCGGTGTCCAACCGTGCGCCTTTCG‑ACGAGCAGAGctgcct                                                                                            >  1:63628/1‑145 (MQ=255)
                                                              ccaGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGCATTCACGGGGCCCACGCCGGTGTCCCACCGTGCGCCTTTCG‑ACGCGCCGAGCAGCCGGCAGGTc                                                                                     >  7:353753/1‑145 (MQ=255)
                                                              ccaGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCG‑ACGAGCAGAGCTGCCTGCTGGTc                                                                                     >  1:281478/1‑145 (MQ=255)
                                                                aGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTAAACGCGCCCCCCCCGGTGTCCAACCCGGCGCCGTTAG‑ACGAGCCGAGCGGCCGGGTGGTCGCAg                                                                                 >  4:242741/1‑147 (MQ=255)
                                                                  ccgccAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTGGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTAG‑CCGAGCAGAGCTGCCTGCTGGTCGCAGCg                                                                               >  1:234256/1‑147 (MQ=255)
                                                                             cGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTAACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTAG‑ACGAGCAGAGCTGGCTGCTGGTCGAAGCGTTGCGCCCGGa                                                                    >  3:334672/1‑147 (MQ=255)
                                                                              gATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGCTTTCACGGCGCCCCCGCCGGTGTCCAACCGTGCGCCTTTCG‑ACGAGCAGAGCTGCCTGCTGGTCGCAGCGTTGCGCCCGGAg                                                                   >  7:43950/1‑147 (MQ=255)
                                                                                    gAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGCTTTCACCGCGCCCCCGCCGGTGTGCAACCGTGCGCCTTTCG‑ACCAGCAGAGCTGCCTGCGGGTCGAAGCGTTGCGCCCCGAggggg                                                               >  1:145107/1‑145 (MQ=255)
                                                                                     cgcgcGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGGGTGCAACCGTGCGCCTTTCG‑ACGCGCAGAGCTGCCTGCTGGTCGCAGCGTTGCGCCCGGAGGTGGTGa                                                            >  5:176583/2‑147 (MQ=255)
                                                                                     aGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTTCGCCTTTCG‑ACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGa                                                            >  7:68359/1‑147 (MQ=255)
                                                                                                cgcTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCCCCCAAATAAACCCGCCCAAGCCCCTGTGAAAACGAGCGCCATTAG‑ACCAACCCAGCTGCCAGCTGGTGGCAACCTTGCGCCCCGAAGTGGTGAGCGTCCAatta                                                 >  8:269487/1‑146 (MQ=255)
                                                                                                                       ttattaCAGCGAAGTGGGCGCCGATTTTACGGCGCCCACCCCGGTGTACAACCCTGCGCCTTTAG‑ACGAGCAGAGCTGCCGGCTGGTGGAAGCGGTGCGCCCGGAGGGGGGGAGCTTACATTTAGGCCAGCCGCAGGCCGAAcagct                          >  4:144520/1‑147 (MQ=255)
                                                                                                                       ttattaCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCG‑CCGAGCAGAGCTGCCTGCTGGTCGCAGCGTTGCGCCCGGAGGTGGTGAGCTTACATTTCGGCCTGCCGCAGGCCGAActgct                          >  3:285463/1‑147 (MQ=255)
                                                                                                                               gCGAAGTGGGCGCCGATTTAACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTAG‑ACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCtt                                                        >  7:177687/1‑109 (MQ=255)
                                                                                                                                 gAAGTGGGCGCCGATTTAACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCG‑ACGAGCAGAGCTGCCTGCTGGGCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCCGGCCGAACTGCTGCAGCGGGt                 >  8:168781/1‑146 (MQ=255)
                                                                                                                                 gAAGTGGGCGCCGATTTAACCGCGCCCACGCCGGTGTACAAACGTGCGCCTTTAGAACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCggg                  >  1:105729/1‑146 (MQ=255)
                                                                                                                                   aGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCG‑ACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTaaa              >  4:294968/1‑147 (MQ=255)
                                                                                                                                             gATTTCACGGCGCCCACGCCGGTGTCCTACCGTGCGCCTTTCG‑ACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGt    >  1:216855/1‑147 (MQ=255)
                                                                                                                                                 tCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCG‑ACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGACCTGCTGCAGCGGGTAAAGGCCAGCGGTGc  >  2:178212/1‑145 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GAGATCGCTGCATTCCGCGCCGGCTGCCCGGGGCGCCCGCTGAACCTGAACTTCTTCTGCCACCAGCCGCCAGCGCCCGATGCCGAGCGCGATGCGCGCTGGAAGCAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCG‑ACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGC  >  NC_002947/4353046‑4353335

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