Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,812,184	A→T	7.5%	intergenic (+104/‑289)	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,184	0	A	T	7.5%	153.0 / 3.0	53	intergenic (+104/‑289)	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 T (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.35e-01

ACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATC  >  NC_002947/4812037‑4812330                                                                                                                                                     |                                                                                                                                                   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                                                                                                  >  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)                                                                                                                                                     |                                                                                                                                                   ACACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACC‑CTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAACCATTATC  >  NC_002947/4812037‑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: ‑