Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091359,6280TC21.7% 83.6 / 21.4 46intergenic (+349/+59)yabP/rluApseudogene, pentapeptide repeats‑containing/dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (17/19);  new base C (5/5);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.76e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAA  >  NC_000913/59489‑59771
                                                                                                                                           |                                                                                                                                               
aaTGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGAttt                                                                                                                                        <  2:54671/149‑1 (MQ=255)
  tGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTAt                                                                                                                                      >  1:114608/1‑149 (MQ=255)
  tGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTAt                                                                                                                                      <  2:114608/149‑1 (MQ=255)
    aaaTATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGt                                                                                                                                    <  3:31488/149‑1 (MQ=255)
     aaTATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAttttttgt                                                                                                                                   <  6:53999/149‑6 (MQ=255)
                     aaCGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAATGGAATTACCGAATGGAAAGAGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTa                                                                                                                   <  3:18162/149‑1 (MQ=255)
                      aCGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAg                                                                                                                  <  5:1764/149‑1 (MQ=255)
                             aagaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAgg                                                                                                           <  4:28048/149‑1 (MQ=255)
                               gaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTa                                                                                                                   <  3:12275/139‑1 (MQ=255)
                               gaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTa                                                                                                                   >  4:12275/1‑139 (MQ=255)
                               gaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAggg                                                                                                          >  3:11418/1‑148 (MQ=255)
                               gaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAtttttt                                                                                                                                     >  3:17191/1‑118 (MQ=255)
                               gaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAtttttt                                                                                                                                     >  3:15204/1‑118 (MQ=255)
                               gaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAtttttt                                                                                                                                     <  4:15204/121‑4 (MQ=255)
                               gaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAtttttt                                                                                                                                     <  4:17191/121‑4 (MQ=255)
                                aaaaaGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTAt                                                                                                                                      >  3:65296/1‑119 (MQ=255)
                                aaaaaGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTAt                                                                                                                                      <  4:65296/119‑1 (MQ=255)
                                           aaaTTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAgg                                                                                                                            >  1:113827/1‑118 (MQ=255)
                                           aaaTTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAgg                                                                                                                            <  2:113827/118‑1 (MQ=255)
                                                   tgtgGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGt                                                                                      >  8:72299/1‑148 (MQ=255)
                                                     tgGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTaa                                                                                   <  4:11418/149‑1 (MQ=255)
                                                               aGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTg                                                                         >  4:120606/1‑149 (MQ=255)
                                                                   aaTTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATcc                                                                             <  8:133687/141‑1 (MQ=255)
                                                                   aaTTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATcc                                                                             >  7:133687/1‑141 (MQ=255)
                                                                           aaTGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACg                                                             >  2:142036/1‑149 (MQ=255)
                                                                                    cGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACTAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGTAAGGTTAAAAATCCGCTGGCGCTTTATACGTCATACTAt                                                    <  3:83272/149‑1 (MQ=255)
                                                                                       aGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGc                                                 >  6:92033/1‑149 (MQ=255)
                                                                                             acGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCtt                                                                   <  6:97539/125‑1 (MQ=255)
                                                                                             acGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCtt                                                                   >  5:97539/1‑125 (MQ=255)
                                                                                                   caacTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGAt                                     <  5:28316/149‑1 (MQ=255)
                                                                                                       tGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATggg                                  >  4:138610/1‑148 (MQ=255)
                                                                                                       tGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATggg                                  >  2:31807/1‑148 (MQ=255)
                                                                                                                     cAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACAtct                   >  1:142008/1‑149 (MQ=255)
                                                                                                                         aaCAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGc               <  1:142036/149‑1 (MQ=255)
                                                                                                                             ataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTAc                                                                                               <  5:33915/63‑1 (MQ=25)
                                                                                                                             ataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTAc                                                                                               >  6:33915/3‑65 (MQ=25)
                                                                                                                             ataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGAt                                     <  1:85826/121‑1 (MQ=255)
                                                                                                                             ataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGAt                                     >  2:85826/3‑123 (MQ=255)
                                                                                                                             ataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATg            >  2:57981/3‑148 (MQ=255)
                                                                                                                             aaaaaTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGt           <  3:120606/149‑1 (MQ=255)
                                                                                                                               aaaTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAg         <  3:138610/149‑1 (MQ=255)
                                                                                                                                   ccccGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATggg                                  <  4:137352/120‑1 (MQ=255)
                                                                                                                                   ccccGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATggg                                  >  3:137352/1‑120 (MQ=255)
                                                                                                                                      cGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACaaa  >  1:65740/1‑149 (MQ=255)
                                                                                                                                       ggACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCaa                        >  5:3776/1‑126 (MQ=255)
                                                                                                                                       ggACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCaa                        <  6:3776/126‑1 (MQ=255)
                                                                                                                                           tCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCa                                                          >  7:4660/1‑88 (MQ=255)
                                                                                                                                           tCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCa                                                          <  8:4660/88‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                               
AATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACAAA  >  NC_000913/59489‑59771

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