Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,262,1710GC21.7% 35.5 / 11.4 23S196T (AGC→ACC) dusAtRNA‑dihydrouridine synthase A
Reads supporting (aligned to +/- strand):  ref base G (11/7);  new base C (3/2);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTG  >  NC_000913_3_pae_tpiA/4262023‑4262319
                                                                                                                                                    |                                                                                                                                                    
atgATTCTGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAg                                                                                                                                                      >  2:261171/3‑149 (MQ=255)
   agtgcggtgacggtgaAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGccc                                                                                                                                                   <  2:581347/145‑1 (MQ=255)
              ggtgaAAGCGCGTATTGGCATCGACGACCAGGACAGGTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGTTTTGTTTATCATCGACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCACGAAAGGAAAcc                                                                                                                                        >  2:250597/1‑149 (MQ=255)
                       ggtaagtGGTATCGACGGCAGGGACGGCGATGAAGTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATcc                                                                                                                               <  2:590003/143‑1 (MQ=255)
                                       gacCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATcc                                                                                                               >  1:540399/1‑149 (MQ=255)
                                                   tATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGACTATCCGCATGTCTATCa                                                                                                   >  2:429520/1‑149 (MQ=255)
                                                               tGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGa                                                                                       >  2:507189/1‑149 (MQ=255)
                                                                                      cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaa                                                                                                                                               >  2:618344/1‑70 (MQ=39)
                                                                                      cGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGCTTAACCCCgaaa                                                                                                                                               <  1:618344/70‑1 (MQ=39)
                                                                                              gCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGCCTTTCCGCATCTCACAATGTCGATTAACggt                                                        >  2:35548/1‑149 (MQ=255)
                                                                                                        gaTGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAACCGTGACTTTCCGCATCTGACAATGTCGATTCAAGTGGGTATCCAGt                                              >  2:275439/1‑149 (MQ=255)
                                                                                                            ttCATCATCCACGCACGTAAAGCCTGGCTTTCGGTGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCt                                          >  2:686447/1‑149 (MQ=255)
                                                                                                                    ccACGCACGTAAAGCCGGGGTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTACCGGTGGTATCAAGTCGCTGGAAGAgg                                  <  2:394821/149‑1 (MQ=255)
                                                                                                                     cacgcaggTAAAGCCTGGGTTTCGGGGTTAAGCACGAAAGAAAACGGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCTACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGc                                 <  2:616648/149‑1 (MQ=255)
                                                                                                                             tAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCACGTCGCTGGAAGAGGCCAAAGcac                         >  2:92344/1‑149 (MQ=255)
                                                                                                                              aaaGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGAGTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCCACTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGcaca                        <  2:42860/149‑1 (MQ=255)
                                                                                                                                 gCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCt                     >  2:533251/1‑149 (MQ=255)
                                                                                                                                   cTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGCGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGCCCAAAGCACACCTGc                   >  1:657343/1‑149 (MQ=255)
                                                                                                                                     ttcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCaa                 >  2:386246/3‑149 (MQ=255)
                                                                                                                                      tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAgg                                  <  2:609343/130‑1 (MQ=255)
                                                                                                                                      tcTTTCGGGCTTAACCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAgg                                  >  1:609343/2‑131 (MQ=255)
                                                                                                                                             gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGTAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATatgga         <  1:616075/149‑1 (MQ=255)
                                                                                                                                             gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCt                                                                           <  2:43234/83‑1 (MQ=255)
                                                                                                                                             gggTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCt                                                                           >  1:43234/1‑83 (MQ=255)
                                                                                                                                                    gCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGATGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTg  >  2:477891/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGATTCCGGTGACGGTGAAAACGCGTATTGGCATCGACGACCAGGACAGCTATGAATTTCTCTGCGATTTCATCAATACCGTTTCCGGCAAAGGCGAGTGTGAGATGTTCATCATCCACGCACGTAAAGCCTGGCTTTCGGGGTTAAGCCCGAAAGAAAACCGTGAAATCCCGCCGCTCGATTATCCGCGTGTGTATCAACTGAAGCGTGACTTTCCGCATCTGACAATGTCGATTAACGGTGGTATCAAGTCGCTGGAAGAGGCCAAAGCACACCTGCAACATATGGATGGCGTG  >  NC_000913_3_pae_tpiA/4262023‑4262319

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