Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,814,861	C→T	10.9%	noncoding (60/78 nt)	micA →	sRNA antisense regulator of ompA, lamB, ompX, and phoP; Hfq‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,814,861	0	C	T	10.9%	164.0 / 7.6	55	noncoding (60/78 nt)	micA	sRNA antisense regulator of ompA, lamB, ompX, and phoP; Hfq‑dependent
Reads supporting (aligned to +/- strand):  ref base C (26/23);  new base T (4/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGC  >  NC_000913/2814713‑2815006                                                                                                                                                     |                                                                                                                                                  gAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                                                                                                                   <  2:686362/149‑1 (MQ=255)          tAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCt                                                                                                                                          >  4:426667/1‑149 (MQ=255)          tAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCGGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCt                                                                                                                                          >  2:90027/1‑149 (MQ=255)             caacGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCGTGAGTGGccaaaa                                                                                                                                       <  2:250525/149‑5 (MQ=255)                cGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttct                                                                                                                                    <  2:548332/149‑1 (MQ=255)                            aCCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGc                                                                                                                        >  2:246777/1‑149 (MQ=255)                             ccTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCg                                                                                                                       >  2:668861/1‑149 (MQ=255)                               tCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgtg                                                                                                                     >  2:553025/1‑149 (MQ=255)                                       tttttttAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCa                                                                                                             <  3:574081/149‑1 (MQ=255)                                                    tttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaaa                                                                                                <  1:668861/149‑1 (MQ=255)                                                          aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTc                                                                                          >  3:586102/1‑149 (MQ=255)                                                           aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCg                                                                                         <  3:420412/149‑1 (MQ=255)                                                            ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                                                                                                      <  3:563483/103‑1 (MQ=255)                                                            ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCgg                                                                                        <  1:157955/149‑1 (MQ=255)                                                            ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTGATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                                                                                                                      >  4:563483/1‑103 (MQ=255)                                                               tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttt                                                                                                                                  >  3:189571/1‑104 (MQ=255)                                                               tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttt                                                                                                                                  <  4:189571/104‑1 (MQ=255)                                                                                       tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTcc                                                                    >  1:680249/1‑142 (MQ=255)                                                                                       tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTcc                                                                    <  2:680249/142‑1 (MQ=255)                                                                                         tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCtt                                                                            >  3:253666/1‑132 (MQ=255)                                                                                         tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCtt                                                                            <  4:253666/132‑1 (MQ=255)                                                                                         tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACgg                                                           >  2:689205/1‑149 (MQ=255)                                                                                          gAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGc                                                          >  1:345598/1‑149 (MQ=255)                                                                                          gAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGc                                                          >  3:20427/1‑149 (MQ=255)                                                                                           aaaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATGTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCg                                                         >  4:299134/1‑149 (MQ=255)                                                                                            aaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCTCTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAACCACACCGCAAACGGTTCGGTATCAGCGGCTTTCATTTCCTGCTGACGGCGt                                                        <  4:586102/149‑1 (MQ=255)                                                                                                 gcgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGa                                                   <  3:426667/149‑1 (MQ=255)                                                                                                 gcgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGa                                                   >  1:484463/1‑149 (MQ=255)                                                                                                 gcgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGa                                                   <  1:406730/149‑1 (MQ=255)                                                                                                 gcgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGa                                                   <  1:90027/149‑1 (MQ=255)                                                                                                  cgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGNNNGCGTTCAGAc                                                  <  4:280751/149‑1 (MQ=255)                                                                                                  cgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGAc                                                  <  1:689205/149‑1 (MQ=255)                                                                                                      tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaaa                                                                                                <  2:639387/99‑1 (MQ=255)                                                                                                      tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCaaa                                                                                                >  1:639387/1‑99 (MQ=255)                                                                                                                tcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttt                                                                                                                  <  1:295633/71‑1 (MQ=255)                                                                                                                tcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttt                                                                                                                  >  2:295633/1‑71 (MQ=255)                                                                                                                tcatcCATGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGc                                                >  2:581607/1‑137 (MQ=255)                                                                                                                tcatcCATGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGc                                                <  1:581607/137‑1 (MQ=255)                                                                                                                      cTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACaaa                              >  1:121275/1‑149 (MQ=255)                                                                                                                                  aTGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGAc                                                             <  2:349640/106‑1 (MQ=255)                                                                                                                                  aTGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGAc                                                             >  1:349640/1‑106 (MQ=255)                                                                                                                                      aaTTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGtt                                                       <  2:458930/108‑1 (MQ=255)                                                                                                                                      aaTTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGtt                                                       >  1:458930/1‑108 (MQ=255)                                                                                                                                       aTTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGaa             <  1:553025/149‑1 (MQ=255)                                                                                                                                         aaaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAAtt           >  1:547466/4‑149 (MQ=255)                                                                                                                                          aaggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAAttt          >  2:642367/3‑149 (MQ=255)                                                                                                                                           aggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTc         >  2:561630/2‑149 (MQ=255)                                                                                                                                            ggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCgcgc                 <  3:478499/140‑1 (MQ=255)                                                                                                                                            ggCCACTCGTGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCgcgc                 >  4:478499/1‑140 (MQ=255)                                                                                                                                             gCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCa       <  4:20427/149‑1 (MQ=255)                                                                                                                                             gCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCa       <  1:246777/149‑1 (MQ=255)                                                                                                                                               cACTCACGAGTGGCCTTTTTCTTTTCTGTGAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGc     >  4:38422/1‑149 (MQ=255)                                                                                                                                               cACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGc     >  1:412316/1‑149 (MQ=255)                                                                                                                                               cACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGc     >  4:60909/1‑149 (MQ=255)                                                                                                                                                aCTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCg    >  2:352086/1‑149 (MQ=255)                                                                                                                                                  tCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGc  >  3:479214/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GAAGCTATCTAACAACGGCATTTAGCCACCTCCGGTAATTTTTTTAAAAATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGC  >  NC_000913/2814713‑2815006

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