Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,812,185	A→G	7.5%	intergenic (+105/‑288)	PP_4238 → / → PP_5665	hypothetical protein/McbA family microcin B17‑like processing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,812,185	0	A	G	7.5%	146.1 / 2.8	53	intergenic (+105/‑288)	PP_4238/PP_5665	hypothetical protein/McbA family microcin B17‑like processing protein
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base G (2/2);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

CACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATC  >  NC_002947/4812038‑4812330                                                                                                                                                     |                                                                                                                                                  cacaTTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                                  <  8:266215/149‑1 (MQ=255)  acaTTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑GTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCaaa                                                                                                                                                  >  2:186448/1‑148 (MQ=255)    aTTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACa                                                                                                                                                >  1:334774/1‑148 (MQ=255)     ttGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACAt                                                                                                                                               >  4:61084/1‑148 (MQ=255)       ggTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                             <  4:278526/148‑1 (MQ=255)       ggTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                             <  3:61084/148‑1 (MQ=255)        gTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAg                                                                                                                                            >  1:20622/1‑148 (MQ=255)               gACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTcc                                                                                                                                    >  2:374041/1‑148 (MQ=255)               gACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTc                                                                                                                                     <  6:212879/148‑1 (MQ=255)                  gCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTccgaa                                                                                                                                 <  1:374041/149‑5 (MQ=255)                         gTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCttat                                                                                                                          <  5:228928/149‑1 (MQ=255)                          ttCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCttatt                                                                                                                         <  4:316859/149‑1 (MQ=255)                               tCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTaa                                                                                                                       <  4:266703/146‑1 (MQ=255)                                 gCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCa                                                                                                                                          <  1:70199/125‑1 (MQ=255)                                 gCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCa                                                                                                                                          >  2:70199/1‑125 (MQ=255)                                  cAAGCAGTAAATCTCCCGAAATCAATCTCACGCCCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACtt                                                                                                                  >  2:262399/1‑148 (MQ=255)                                   aaGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTc                                                                                                                 >  2:338753/1‑148 (MQ=255)                                         tAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                             <  6:354202/114‑1 (MQ=255)                                         tAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCa                                                                                                                                             >  5:354202/1‑114 (MQ=255)                                                  cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATTTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaa                                                                                                 <  1:290011/149‑1 (MQ=255)                                                   gAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa                                                                                                <  2:20622/149‑1 (MQ=255)                                                      atcattcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGCAACAATTCTTCGTAGCATGCGGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGaa                                                                                             <  6:366970/149‑1 (MQ=255)                                                      atcaatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTACAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGaa                                                                                             <  6:245467/149‑1 (MQ=255)                                                        caatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCAt                                                                                                                                         <  8:116103/103‑1 (MQ=255)                                                        caatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCAt                                                                                                                                         >  7:116103/1‑103 (MQ=255)                                                             tcACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTca                                                                                       >  6:332280/1‑148 (MQ=255)                                                               aCGACCCGGCAAATCTATCCCCACGCAGGTCTAACCGCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTc                                                                                        >  1:253562/1‑146 (MQ=255)                                                               aCGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTc                                                                                        <  2:253562/145‑1 (MQ=255)                                                                  aCCCGGCAATCCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAAGAACTAGCTACAATTCTTCTGTCCATGCTGATGTCAAACCTCAGCAAGCTCGCGCCCCTTTTTAAAACTACACCCCGGGAATAACAAAGAACACACACACGCt                                                                                 >  7:237730/1‑149 (MQ=255)                                                                   cccGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACa                                                                                                  >  8:158170/1‑131 (MQ=255)                                                                   cccGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACa                                                                                                  <  7:158170/131‑1 (MQ=255)                                                                   cccGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCATCGGAGCATGCTGATGTCAAACATCGGCATGCTCGCGCCCCGTATTAAAACTTCACACCGCGCATGACAAAGAACAGTCACACGCTc                                                                                >  8:57775/1‑149 (MQ=255)                                                                        cAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcacca                                                                            >  6:178919/1‑148 (MQ=255)                                                                           aTCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACaga                                                                        <  2:334774/149‑1 (MQ=255)                                                                           aTCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACaga                                                                        <  1:338753/149‑1 (MQ=255)                                                                            tCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACcaca                                                                          >  1:262411/1‑146 (MQ=255)                                                                               aTCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTcac                                                                                      >  8:100885/1‑131 (MQ=255)                                                                               aTCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTcac                                                                                      <  7:100885/131‑1 (MQ=255)                                                                                  cccACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTc                                                                  <  5:178919/148‑1 (MQ=255)                                                                                           gTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTc                                                                                                                 <  2:332043/92‑1 (MQ=255)                                                                                           gTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTc                                                                                                                 >  1:332043/1‑92 (MQ=255)                                                                                                             aaaaTATCTGGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTca                                                                               <  2:22326/109‑1 (MQ=255)                                                                                                             aaaaTATCTGGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTca                                                                               >  1:22326/1‑109 (MQ=255)                                                                                                             aaaaTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTc                                                                  >  7:89708/1‑122 (MQ=255)                                                                                                             aaaaTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTc                                                                  <  8:89708/122‑1 (MQ=255)                                                                                                                     tAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTatatttatgt                                >  1:177821/1‑142 (MQ=255)                                                                                                                     tAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAatat                                >  5:295936/1‑148 (MQ=255)                                                                                                                              ttcttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACt                      >  4:210090/1‑149 (MQ=255)                                                                                                                                cttcGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTATATTTATGTAAGAAGcttct                     <  2:262411/148‑5 (MQ=255)                                                                                                                                   cGGAGCATGCTTATTTAAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTATATTTATGTAAGAAGcttctg                    <  2:177821/146‑6 (MQ=255)                                                                                                                                        cATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTaaa            <  6:295936/149‑1 (MQ=255)                                                                                                                                         aTGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTaaaa           <  1:186448/149‑1 (MQ=255)                                                                                                                                          tGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTaaaaa          >  4:133511/1‑149 (MQ=255)                                                                                                                                             tGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCa       >  4:154545/1‑149 (MQ=255)                                                                                                                                                tGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTa    <  3:133511/149‑1 (MQ=255)                                                                                                                                                    aaaCATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATc  <  7:162238/147‑1 (MQ=255)                                                                                                                                                    aaaCATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATc  >  8:162238/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATC  >  NC_002947/4812038‑4812330

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑