Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	164,581	A→C	19.0%	intergenic (+47/‑149)	hrpB → / → mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	164,581	0	A	C	19.0%	74.0 / 2.1	47	intergenic (+47/‑149)	hrpB/mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (16/21);  new base C (2/8);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGA  >  NC_000913/164432‑164721                                                                                                                                                      |                                                                                                                                             ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                              <  2:557411/150‑1 (MQ=255)  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                             <  2:1214270/150‑1 (MQ=255)  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                             <  2:489338/150‑1 (MQ=255)  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                             <  2:489211/150‑1 (MQ=255)   gTGAGGTGCAAAAAGAGATGAAAGGGTGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                            <  2:1932329/150‑1 (MQ=255)         tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                      >  2:2033607/1‑150 (MQ=255)         tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                      >  2:1753708/1‑150 (MQ=255)                ttataTGAAAGGGCGTTATCCCAAACATTTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGAGGGACGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTtc                                                                                                                               <  2:1265292/146‑1 (MQ=38)                   gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggacc                                                                                                                            >  2:256246/1‑145 (MQ=255)                      gAAATGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccgc                                                                                                                         >  2:506923/1‑148 (MQ=255)                       aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTgcacccaca                                                                                                                        >  2:365159/1‑150 (MQ=255)                         aGGGCGTTATCCCAAACATGTCGGGCCGGACGACCCGCAAAATACTGCACCGACGCGACGGACTAAAAATTTTTGTTATGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctcacc                                                                                                                      <  2:179917/150‑1 (MQ=37)                            gCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgaccctgaaccat                                                                                                                   >  1:790969/1‑143 (MQ=255)                                    caaaaCAGGTCGGGCCGCACGACCCGCCAAATCCTCCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGTTTTCCCCCCCCGCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                           <  2:1868152/148‑1 (MQ=255)                                      aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCGGCGCGCACGTGGACCCGCCCCCCACCCCACAccc                                                                                                         >  2:2326825/1‑150 (MQ=37)                                      aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACACCCCAAaccacaacc                                                                                                         >  1:1631327/1‑144 (MQ=37)                                      aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCccacccccaccc                                                                                                         >  2:32298/1‑150 (MQ=37)                                      aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCACCCCAAaccacaacc                                                                                                         >  2:1917399/1‑144 (MQ=37)                                             tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGGCCCCCACCCCCACCCACCCCCCCaaacg                                                                                                  >  1:2201168/1‑147 (MQ=37)                                                 gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACACCCTCAACGgagtg                                                                                              >  2:1201210/1‑148 (MQ=37)                                                         acCCGGCAAATACTGCACCGACGCGCCGCCCGAAAACTTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttaggg                                                                                      <  2:1713833/150‑7 (MQ=37)                                                         acCCGCCAAATACTGCACCGACGCGCCGGACTAAAAAGTATTCGTATGGTTGGTTTTCGCCCTCCCCCTGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgggttaggg                                                                                      <  2:1713900/150‑7 (MQ=255)                                                            cgccAAACACTGCCCCGCCGCGACGGACCAAACAGTTTTCTTAGGGTTGTTTTTCTCCCCCCCCCGGGGGGAGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                   <  2:474353/147‑1 (MQ=255)                                                            cGGCACATACTGCACCGACGCGACGGACGAAAAAGTCTTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGCCCCCCCCCCCCAACCCCCCCCCGCCAGGGGGCGGGGGCCgcgcg                                                                                   >  2:974224/1‑146 (MQ=21)                                                            cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCACCCCCCCCCGAACCCACGCCCCCCACGGCGCGGGGCCCgcgcg                                                                                   >  1:1963967/1‑146 (MQ=25)                                                              gCAAATACTGCACCGACGCGCCGGCCGAAAAAGTATTCGTAAGGTTGGTTTTCCCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                 <  2:1355360/150‑1 (MQ=255)                                                              gCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                 <  2:589959/150‑1 (MQ=255)                                                                          ccgccgGGACGGCCGAAAATGTTTTCGTACGTTTGTTTTTAGCCCTCTCCGTGTGGTAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                     <  1:335966/150‑1 (MQ=255)                                                                             acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCACTCAAGGGAGAGGGGACCGATCGAGCACACATTTTgagag                                                                  >  2:1538618/1‑150 (MQ=255)                                                                                 gacgtacgAAAATTTTTTCGTATGTTTGGTTTTCTCCCTCTCCCTGTGGTAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                              <  2:2171796/150‑1 (MQ=255)                                                                                       cgAACATGTATGCTTATGGTTGGTTTACTCCCCCTCCCTGTGGGAGTGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGTGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                        <  1:1873300/150‑1 (MQ=255)                                                                                        taaaatGTATTCGTACGTTTGGTTTTCCCCCTCTCCGTGGGGGTGCGGGCGGGGTTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                       <  2:1899074/149‑1 (MQ=255)                                                                                        gAAAATTTATTGGTATGTTTGTTTTACCCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                       <  2:1898929/150‑1 (MQ=255)                                                                                           aatgttttCGTATGGTTGGTTTTCGCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                    <  2:350855/144‑1 (MQ=38)                                                                                           aaaGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                    <  2:1554494/150‑1 (MQ=255)                                                                                            atgTATTCGTACGTTTGTTTCTCACCCCCGCCGGGTGGGAGGGGGCGGGGGGGCGGCCATCAGCGCGCACGTTCCCCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                   <  1:208588/148‑1 (MQ=25)                                                                                            aaGTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                   <  1:1768236/150‑1 (MQ=255)                                                                                              ttATTCGTACGGTTGGTTTTCCCCCCCGCCCGGTGGTAGCGGGCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                 <  1:1843432/149‑1 (MQ=255)                                                                                              gTATTCGTCAGGTTGCTGTCCGCCCGCGCCCGGGGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTCCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                 <  1:1670457/150‑1 (MQ=255)                                                                                                  tcttctggttgtttTCCCCCCTCGCCGGGGGGGAGAGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                             <  2:1043141/144‑1 (MQ=255)                                                                                                  tCGTAAGGGTGGTTTTCGCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                             <  2:1528716/150‑1 (MQ=255)                                                                                                                       tctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACAACCACAAGGGGGAGGGGGCCCATAGCGCCCAAAATTTGGGGGAAATAATCTTCTGTCTTGTAAAAGAAGAAAAGAAAAACGGGc                        >  2:1570903/1‑150 (MQ=16)                                                                                                                         tccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCGAACCCACACCCACAAGGGAGAGGGGACCCGTCGAGCACAAATTTTGGGAGATATCTTCTTCTGTCTTGGGACAGAAGAACAGAAAAGCGGGcgt                      >  2:2182507/1‑148 (MQ=38)                                                                                                                         tccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCAAACCCACAACCTCAAGGGGGTGGGGGCCGCTCGAGCACCAATATTGAGAGGTATATTCTTATGTGTTTTTACACAACAACCGACAATCGGGcgt                      >  2:2140538/1‑148 (MQ=19)                                                                                                                           cctGGGCGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                    <  2:2230011/150‑1 (MQ=255)                                                                                                                           ccggtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                    <  2:1779002/147‑1 (MQ=255)                                                                                                                                           ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg    <  1:755323/150‑1 (MQ=255)                                                                                                                                             ggTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga  <  1:23925/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                             CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGA  >  NC_000913/164432‑164721

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