Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	888,014	0	T	G	39.5%	45.6 / 16.9	38	noncoding (44/81 nt)	rybB	small regulatory RNA RybB
Reads supporting (aligned to +/- strand):  ref base T (7/16);  new base G (15/0);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 1.96e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATA  >  NC_000913/887877‑888153                                                                                                                                               |                                                                                                                                             ggATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGTTAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTg                                                                                                                                            >  1:498261/1‑140 (MQ=255) ggATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTg                                                                                                                                            >  1:66746/1‑140 (MQ=255)    tGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATg                                                                                                                                         >  1:169069/1‑140 (MQ=255)          tGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTcc                                                                                                                                   <  2:500054/140‑1 (MQ=255)          tGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTcc                                                                                                                                   >  2:215763/1‑140 (MQ=255)           gCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCca                                                                                                                                  <  2:66746/140‑1 (MQ=255)            cATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCcac                                                                                                                                 >  2:629022/1‑140 (MQ=255)                    cACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATggg                                                                                                                        <  1:415316/140‑1 (MQ=255)                         accgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACAt                                                                                                                   <  1:360303/140‑1 (MQ=255)                           cgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATggg                                                                                                                        <  1:146235/133‑1 (MQ=255)                           cgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAATggg                                                                                                                        >  2:146235/1‑133 (MQ=255)                                    cTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAATGGGGACATCaaaaaaaaac                                                                                                        >  1:197176/1‑138 (MQ=255)                                    cTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAATGGGGACATCAAAGaaaaac                                                                                                        >  2:83530/1‑138 (MQ=255)                                    cTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAATGGGGACATCAAAAAAAAga                                                                                                        >  2:734490/1‑139 (MQ=255)                                     tATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCa                                                                                                       >  2:342976/1‑140 (MQ=255)                                      aTCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGTTAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAg                                                                                                      <  2:498261/140‑1 (MQ=255)                                          cgcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAATGGGGGCATCAAAGAAAAACAGTggg                                                                                                  >  2:760843/1‑139 (MQ=255)                                           gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCa                                                                                                 <  1:215763/140‑1 (MQ=255)                                           gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCa                                                                                                 <  2:163096/140‑1 (MQ=255)                                           gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCAC‑‑AAATGGGGACATCAAAGAAAAGCAGTGGCa                                                                                                 <  2:709532/139‑1 (MQ=255)                                           gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCAC‑‑AAATGGGGACATCAAAGAAAAGCAGTGGCa                                                                                                 >  1:709532/1‑139 (MQ=255)                                                                          tGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACTGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAAATGGGACATCAAAAAAAAGCAGGGGCACTAATTAAAACTGATGCCCTGCGGaa                                                                       >  1:613562/1‑140 (MQ=255)                                                                          tGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAAAGGGGGCATCAAAGAAAAACAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGtt                                                                  >  1:181102/1‑140 (MQ=255)                                                                                       aTAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCa                                                     >  1:602122/1‑140 (MQ=255)                                                                                         aaGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAAAAGGGACAACAAAAAAAAACAGTGGCACTAATTAAAACTGATGCCCTGCGGAAAAATTCTTCGGTTGTGCaaa                                                   >  2:426890/1‑140 (MQ=255)                                                                                           gTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGAGGCGCTG‑TCTCTTATACACAAATGAGGCTGCCGACGAAAAGGATAGGGGATAATTAAGACTGATGCCCTGCGGGAAAAATCTGCGGGTGTGCaaaaa                                                 >  2:677179/1‑140 (MQ=255)                                                                                            ttttAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCaaaaaa                                                <  1:629022/140‑1 (MQ=255)                                                                                                 gATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTc                                           >  1:796233/1‑140 (MQ=255)                                                                                                 gATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTc                                           <  2:602122/140‑1 (MQ=255)                                                                                                      cAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGGCATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGGAAAGTTCTGCGGTTGTGCAAAAAAATTTCAtttt                                      >  2:152556/1‑140 (MQ=255)                                                                                                      cAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACCGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCAtttt                                      >  1:246915/1‑140 (MQ=255)                                                                                                      cAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCAtttt                                      >  1:479165/1‑140 (MQ=255)                                                                                                             cccATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCa                               >  2:655576/1‑140 (MQ=255)                                                                                                                aTCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACt                            >  1:267780/1‑140 (MQ=255)                                                                                                                aTCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGGTGATGGGCTCCACA‑AAAAGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGGTGTGCAAAAAAATTTCATTTTCAGGGCAACt                            >  1:794328/1‑140 (MQ=255)                                                                                                                       t‑‑‑‑‑tGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGttt                     <  1:342976/140‑1 (MQ=255)                                                                                                                                    aaaTGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTaa             >  1:21887/1‑140 (MQ=255)                                                                                                                                     aaTGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAAt            <  1:655576/140‑1 (MQ=255)                                                                                                                                     aaTGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAAt            <  2:21887/140‑1 (MQ=255)                                                                                                                                      aTGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATc           <  2:796233/140‑1 (MQ=255)                                                                                                                                       tGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACCGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGAGCAACTTCAGTTTTATCCTAATcc          <  2:246915/140‑1 (MQ=255)                                                                                                                                               ttGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATa  <  1:575943/140‑1 (MQ=255)                                                                                                                                               |                                                                                                                                             GGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCT‑‑‑‑‑TGAACCGAAATGGCGGGGTTGATGGGCTCCACA‑AAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCATA  >  NC_000913/887877‑888153

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