Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091359,6270CT20.0% 81.6 / 14.9 41intergenic (+348/+60)yabP/rluApseudogene, pentapeptide repeats‑containing/dual specificity 23S rRNA pseudouridine(746), tRNA pseudouridine(32) synthase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (14/18);  new base T (5/3);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACA  >  NC_000913/59481‑59769
                                                                                                                                                  |                                                                                                                                              
cccGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTc                                                                                                                                              >  6:60343/1‑149 (MQ=255)
       aaaTGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGAt                                                                                                                                        <  3:92329/148‑1 (MQ=255)
           gAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGGAGAAACAAAAAATCCCCGGACTCGGGGATTTATg                                                                                                                                   >  6:146433/1‑149 (MQ=255)
            aaaTATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGAGTCCGGGGAtttttt                                                                                                                                   <  6:51189/148‑4 (MQ=255)
                  aaGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACaa                                                                                                                              >  8:154662/1‑147 (MQ=255)
                             aaCGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTa                                                                                                                 <  2:66756/149‑1 (MQ=255)
                              aCGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAg                                                                                                                <  4:153051/149‑1 (MQ=255)
                                     aagaaAAATTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCGCCTGAATCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAgg                                                                                                         <  7:154662/149‑1 (MQ=255)
                                     aagaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAgg                                                                                                         <  5:60343/149‑1 (MQ=255)
                                     aagaaAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAAATACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAgg                                                                                                         <  4:127966/149‑1 (MQ=255)
                                                cAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGGAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGt                                                                                               <  5:146433/148‑1 (MQ=255)
                                                     aTTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGa                                                                                         <  5:104329/149‑1 (MQ=255)
                                                         gCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAgg                                                                                     <  6:38781/149‑1 (MQ=255)
                                                                   aGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATcc                                                                           <  5:41546/149‑1 (MQ=255)
                                                                       aGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTg                                                                       <  6:144529/149‑1 (MQ=255)
                                                                             ttACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCtt                                                                 >  5:145176/1‑149 (MQ=255)
                                                                                ccGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATc                                                                            <  1:150369/135‑1 (MQ=255)
                                                                                ccGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATc                                                                            >  2:150369/1‑135 (MQ=255)
                                                                                              tAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTg                                                >  4:67833/1‑149 (MQ=255)
                                                                                                 ccACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTg                                                                       >  2:92802/1‑123 (MQ=255)
                                                                                                 ccACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTg                                                                       <  1:92802/123‑1 (MQ=255)
                                                                                                          acaacTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGa                                    >  2:170405/1‑149 (MQ=255)
                                                                                                             acTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATgg                                 >  6:51943/1‑149 (MQ=255)
                                                                                                              cTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATgg                                 <  1:71110/148‑1 (MQ=255)
                                                                                                              cTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATgg                                 >  2:71110/1‑148 (MQ=255)
                                                                                                                 aaaGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGAt                                   <  6:6222/143‑1 (MQ=255)
                                                                                                                 aaaGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGAt                                   >  5:6222/1‑143 (MQ=255)
                                                                                                                       gTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCa                       <  1:170405/149‑1 (MQ=255)
                                                                                                                        tttGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACTTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCaa                      <  3:67833/149‑1 (MQ=255)
                                                                                                                                   acataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGGGTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCAt           >  7:72904/5‑149 (MQ=255)
                                                                                                                                   acataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCAt           >  1:104049/5‑149 (MQ=255)
                                                                                                                                   acataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTGAAAAATCCGCGGGCGCTGTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCAt           >  7:70770/5‑149 (MQ=255)
                                                                                                                                    aaaaaaTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCtt                                                                 >  3:15081/1‑94 (MQ=255)
                                                                                                                                    aaaaaaTCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCtt                                                                 <  4:15081/94‑1 (MQ=255)
                                                                                                                                     ataaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATg          >  6:6504/3‑148 (MQ=255)
                                                                                                                                      taaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATgg                                 <  4:20028/123‑1 (MQ=255)
                                                                                                                                      taaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATgg                                 >  3:20028/2‑124 (MQ=255)
                                                                                                                                      taaaTCCCCGAGTCCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTa        <  5:6504/148‑1 (MQ=255)
                                                                                                                                             ccGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGcaaca  >  8:39845/1‑148 (MQ=255)
                                                                                                                                              cGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGc      <  5:80400/143‑1 (MQ=255)
                                                                                                                                              cGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGc      >  6:80400/1‑143 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CCCGTTAAAATGAAATATAAGAGACGGTCAACGGGTGAAGAAAAAGTTCAAAAATTCGCTGTGGAGCAGGAAGGGAATTACCGAATGGAAAGCGTAGCCACACGCAACAACTGAAAGCAGTTTGGCAGAAACAAAAAATCCCCGGACTCGGGGATTTATGTACAAGAGGCAGCCCTTAGGATGAGGGTATAAACGTACAGGAAAGGTTAAAAATCCGCTGGCGCTTTAAACGTCATACTATTGCCATACGCCGGATGGGTAATCGTCAACATCTCTGCATGTAGCAACA  >  NC_000913/59481‑59769

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