Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	888,014	0	T	G	34.0%	75.9 / 13.4	53	noncoding (44/81 nt)	rybB	small regulatory RNA RybB
Reads supporting (aligned to +/- strand):  ref base T (9/26);  new base G (18/0);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 9.67e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGC  >  NC_000913/887877‑888149                                                                                                                                          |                                                                                                                                        ggATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTg                                                                                                                                       <  1:756861/140‑1 (MQ=255)     gACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATgg                                                                                                                                   <  1:4626/140‑1 (MQ=255)     gACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATgg                                                                                                                                   >  1:1169034/1‑140 (MQ=255)      aCGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATggg                                                                                                                                  <  2:1395716/140‑1 (MQ=255)       cGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGc                                                                                                                                 <  1:923909/140‑1 (MQ=255)                    cACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATggg                                                                                                                    <  1:1004609/140‑1 (MQ=255)                    cACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATggg                                                                                                                    <  1:1002472/140‑1 (MQ=255)                    cACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATggg                                                                                                                    <  1:824209/140‑1 (MQ=255)                         accgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACAt                                                                                                               <  1:120752/140‑1 (MQ=255)                         accgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGGTGGGCTCCACAAAATGGGGGCAt                                                                                                               >  2:120752/1‑140 (MQ=255)                              aaGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAg                                                                                                          <  1:1352352/140‑1 (MQ=255)                                    cTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGc                                                                                                    >  1:579061/1‑140 (MQ=255)                                    cTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGGCATCAAAGAAAAGc                                                                                                    >  1:893254/1‑140 (MQ=255)                                     tATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCa                                                                                                   >  1:226057/1‑140 (MQ=255)                                        cgcgcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGACATCCAaaaaaaaaaggg                                                                                                >  2:98389/1‑134 (MQ=255)                                           gcgAGGGGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                                                                                             <  2:458219/140‑1 (MQ=255)                                           gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                                                                                             <  2:580465/140‑1 (MQ=255)                                           gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                                                                                             <  2:179147/140‑1 (MQ=255)                                           gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                                                                                             <  2:1288586/140‑1 (MQ=255)                                           gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                                                                                             <  2:1169034/140‑1 (MQ=255)                                                   aTTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTaa                                                                                     <  2:672716/140‑1 (MQ=255)                                                         aGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCCCAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGa                                                                               >  1:1251279/1‑140 (MQ=255)                                                         aGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAAAAGGGGCATCAAAGAAAAGCAGTGGCACTAATTAAGACTGa                                                                               >  2:1121634/1‑140 (MQ=255)                                                             tGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGcc                                                                           <  1:235534/140‑1 (MQ=255)                                                              ggTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGccc                                                                          >  1:1192806/1‑140 (MQ=255)                                                                     agGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGaa                                                                   <  1:457665/140‑1 (MQ=255)                                                                     agGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGaa                                                                   <  1:1121634/140‑1 (MQ=255)                                                                         ttGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGt                                                               >  2:786734/1‑140 (MQ=255)                                                                          tGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTCCACAAAATGGGGGCATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGtt                                                              >  1:178092/1‑140 (MQ=255)                                                                          tGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCCCAAAAAGGGGGCATCAAAGAAAAACAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGtt                                                              >  2:1328695/1‑140 (MQ=255)                                                                                 tagtagATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTCCACAAAATGGGGGCATCAAAAAAAAGCAGGGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTTTGCggg                                                       >  2:1015274/1‑139 (MQ=255)                                                                                 tagtagATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGGAAAGTTCTGCggg                                                       >  2:210231/1‑139 (MQ=255)                                                                                   gtagATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTCCACAAAATGGGGACATCAAAGAAAAACAAGGGCACTAATTAAGACTGATGCCCTGCGGAAAAGATCTGCgggtt                                                     >  2:1327661/1‑137 (MQ=255)                                                                                         aaGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCaaa                                               >  2:106376/1‑140 (MQ=255)                                                                                         aaGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCCCAAAATGGGGGCATCAAAGAAAAACAGGGGGACTAATTAAAACTTATGCCCTGCGGAAAAGGTCTTCGGGTGTGCaaa                                               >  2:659500/1‑140 (MQ=255)                                                                                           gTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCaaaaa                                             >  1:291896/1‑140 (MQ=255)                                                                                           gTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACCGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCaaaaa                                             >  2:405981/1‑140 (MQ=255)                                                                                                aGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAAACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAAttt                                        >  2:1223013/1‑140 (MQ=255)                                                                                                 gATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCCCCACAAAATGGGGACATCAAAAAAAAGCAGGGGCACTAATTAAAACTGATGCCCTGCGGAAAAGTTTTTCGGTTGTGCAAAAAAAtttt                                       >  2:233775/1‑139 (MQ=255)                                                                                                      cAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCAtttt                                  >  1:66338/1‑140 (MQ=255)                                                                                                      cAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCAtttt                                  >  1:852380/1‑140 (MQ=255)                                                                                                      cAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCCCCACAAAAAGGGGACATCAAAAAAAAGCAGTGGCACCAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGGTGTGCAAAAAAAATTCAtttt                                  >  2:1073824/1‑140 (MQ=255)                                                                                                                aTCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACt                        >  1:231916/1‑140 (MQ=255)                                                                                                                 tCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAAAAAAACAGGGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCGACtt                       >  2:527214/1‑140 (MQ=255)                                                                                                                        tGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAAAGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAAAATGATTCCCCGCGGAAAAAATCTGCGGGTGTGCAAAAAAATTTCATTTTTAAGGCAACATCAGtttt                >  2:715081/1‑140 (MQ=255)                                                                                                                         gAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAACAGGGGCACTAATTAAGACTGATGCCCCGCGGAAAAATTCTTCGGGTGTGCAAAAAAAATTCATTTTCAGGGCAACTTCAGttttt               >  2:351550/1‑139 (MQ=255)                                                                                                                           aCCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATc             <  1:405981/140‑1 (MQ=255)                                                                                                                            ccGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATcc            <  1:715081/140‑1 (MQ=255)                                                                                                                                aaTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAAt        <  2:579061/140‑1 (MQ=255)                                                                                                                                aaTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAAt        <  1:351550/140‑1 (MQ=255)                                                                                                                                aaTGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAAt        >  1:1197115/1‑140 (MQ=255)                                                                                                                                 aTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATc       <  2:226057/140‑1 (MQ=255)                                                                                                                                     cGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTgg   <  2:912109/140‑1 (MQ=255)                                                                                                                                     cGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTgg   <  1:98389/140‑1 (MQ=255)                                                                                                                                      ggggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGc  <  2:1304238/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                        GGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGC  >  NC_000913/887877‑888149

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