Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,812,183	C→T	7.2%	intergenic (+103/‑290)	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,183	0	C	T	7.2%	161.2 / 2.0	56	intergenic (+103/‑290)	PP_4238/PP_5665	hypothetical protein/McbA family microcin B17‑like processing protein
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base T (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

AGACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTA  >  NC_002947/4812035‑4812328                                                                                                                                                      |                                                                                                                                                  agaCACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTc                                                                                                                                                   <  7:40669/149‑1 (MQ=255)   acacATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCa                                                                                                                                                  >  3:278526/1‑148 (MQ=255)    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                                                                                                                                       >  7:116103/1‑103 (MQ=255)                                                           caatcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCAt                                                                                                                                       <  8:116103/103‑1 (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                                                                                                <  7:158170/131‑1 (MQ=255)                                                                      cccGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACa                                                                                                >  8:158170/1‑131 (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                                                                             >  1:22326/1‑109 (MQ=255)                                                                                                                aaaaTATCTGGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTca                                                                             <  2:22326/109‑1 (MQ=255)                                                                                                                aaaaTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTc                                                                <  8:89708/122‑1 (MQ=255)                                                                                                                aaaaTATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTc                                                                >  7:89708/1‑122 (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)                                                                                                                                                      |                                                                                                                                                  AGACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTA  >  NC_002947/4812035‑4812328

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