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

GGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCA  >  NC_000913/887876‑888151                                                                                                                                           |                                                                                                                                          gggATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGtt                                                                                                                                          >  2:11122/1‑140 (MQ=255)    aTGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt                                                                                                                                       <  2:908633/140‑1 (MQ=255)    aTGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt                                                                                                                                       <  1:471464/140‑1 (MQ=255)    aTGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGAt                                                                                                                                       >  1:1390922/1‑140 (MQ=255)        cGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGc                                                                                                                                   <  1:1069292/140‑1 (MQ=255)           tGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTcc                                                                                                                                <  2:1406290/140‑1 (MQ=255)           tGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTcc                                                                                                                                <  1:994823/140‑1 (MQ=255)           tGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTcc                                                                                                                                >  2:861630/1‑140 (MQ=255)                          accgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATggg                                                                                                                      <  2:1502467/135‑1 (MQ=255)                          accgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGGAGGGCGCCCCAAAATGGGGACAt                                                                                                                 >  2:936032/1‑140 (MQ=255)                          accgTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATggg                                                                                                                      >  1:1502467/1‑135 (MQ=255)                                       aTCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATTGGCTCCACAAAAAAGGGACATCAAAAAAAAACaa                                                                                                    >  2:773080/1‑139 (MQ=255)                                        tCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCCCAAAAAGGGGGCATCAAAGAAAAGCAgg                                                                                                   >  2:305093/1‑139 (MQ=255)                                            gcgAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCa                                                                                               <  1:209965/140‑1 (MQ=255)                                               aggagATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTa                                                                                            >  1:1052800/1‑140 (MQ=255)                                                  agaTTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAAtt                                                                                         <  1:11122/140‑1 (MQ=255)                                                    aTTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTaa                                                                                       <  2:827342/140‑1 (MQ=255)                                                            gATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGATGATGGGCTCCACAAAATGGGGGGAACGAAAAAAAGCAATGCAAGTAAATAAGACTGAta                                                                               >  2:109761/1‑139 (MQ=255)                                                              tGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGcc                                                                             <  1:245563/140‑1 (MQ=255)                                                                         gTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAg                                                                  >  2:72743/1‑140 (MQ=255)                                                                          ttGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGGCATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGGAAAGt                                                                 >  2:684932/1‑140 (MQ=255)                                                                           tGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGtt                                                                >  1:578936/1‑140 (MQ=255)                                                                                  tagtagATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGt                                                         >  1:920136/1‑140 (MQ=255)                                                                                  tagtagATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAACAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGt                                                         >  2:1340713/1‑140 (MQ=255)                                                                                   agtagATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGtt                                                        >  1:465138/1‑140 (MQ=255)                                                                                       gATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGc                                                    <  2:465138/140‑1 (MQ=255)                                                                                          aaGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCaaa                                                 >  2:901236/1‑140 (MQ=255)                                                                                          aaGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCaaa                                                 >  1:1268697/1‑140 (MQ=255)                                                                                          aaGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCaaa                                                 >  1:930466/1‑140 (MQ=255)                                                                                            gTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGACATCAAAAAAAAGCAGGGGCACTAATTAAAACTGATGCCCTGCGGGAAAGTTCTGCGGGTGTGCaaaaa                                               >  2:140451/1‑140 (MQ=255)                                                                                               ttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAAt                                            >  2:318508/1‑140 (MQ=255)                                                                                                 aGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAAttt                                          >  1:1021286/1‑140 (MQ=255)                                                                                                            aaCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAggg                               <  1:936032/140‑1 (MQ=255)                                                                                                             aCCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGGTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGc                              <  2:578936/140‑1 (MQ=255)                                                                                                               ccATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGGTGTGCAAAAAAATTTCATTTTCAGGGCaa                            >  1:355981/1‑140 (MQ=255)                                                                                                                 aTCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACt                          >  1:1229816/1‑140 (MQ=255)                                                                                                                  tCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACtt                         >  1:218472/1‑140 (MQ=255)                                                                                                                  tCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACtt                         >  1:1309316/1‑140 (MQ=255)                                                                                                                  tCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACtt                         >  1:1229425/1‑140 (MQ=255)                                                                                                                  tCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGGCATCAAAGAAAAGCAGTGGGCCTAATTAAGACTGATGCCCTGCGGGAAAGTTCTGAGGGTGTGCAAAAAAAATTCATTTTTAGGGCAACtt                         >  2:1484848/1‑140 (MQ=255)                                                                                                                      ccTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGt                     <  1:305093/140‑1 (MQ=255)                                                                                                                           aaCCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTAt                >  1:103767/1‑140 (MQ=255)                                                                                                                             ccGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATcc              <  1:72743/140‑1 (MQ=255)                                                                                                                             ccGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATcc              <  1:140451/140‑1 (MQ=255)                                                                                                                                aaaTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGttt                   <  1:540969/132‑1 (MQ=255)                                                                                                                                aaaTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTaa           >  1:81815/1‑140 (MQ=255)                                                                                                                                aaaTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAACAGGGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGttt                   >  2:540969/1‑132 (MQ=255)                                                                                                                                aaaTGGCGGGGGTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTaa           >  1:512996/1‑140 (MQ=255)                                                                                                                                aaaTGGCGGGGGTGATGGGCTCCACAAAAAGGGGGCATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTaa           >  2:671799/1‑140 (MQ=255)                                                                                                                                aaaTGGCGGGGGTGATGGGCTCCACAAAAAGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGGTGTGCAAAAAAATTTTATTTTTAGGGCAACTTCAGTTTTTTCCTaa           >  2:972800/1‑140 (MQ=255)                                                                                                                                  aTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATc         <  2:1229425/140‑1 (MQ=255)                                                                                                                                   tGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATcc        <  2:103767/140‑1 (MQ=255)                                                                                                                                    ggCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCt       <  2:1390922/140‑1 (MQ=255)                                                                                                                                      cGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTgg     <  2:218472/140‑1 (MQ=255)                                                                                                                                       ggggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGc    <  1:901236/140‑1 (MQ=255)                                                                                                                                        gggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGcc   <  2:930466/140‑1 (MQ=255)                                                                                                                                         ggTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCa  <  2:1268697/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                          GGGATGACGCTGCATTTTGTCACCGACCGTAAGCCGCTATCGCGCGAGGAGATTTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGTTGTGCAAAAAAATTTCATTTTCAGGGCAACTTCAGTTTTATCCTAATCCTGGCCA  >  NC_000913/887876‑888151

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: ‑