Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,590,330	G→C	10.0%	G284A (GGC→GCC)	nicP‑I →	porin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,590,330	0	G	C	10.0%	169.2 / 8.0	58	G284A (GGC→GCC)	nicP‑I	porin‑like protein
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base C (2/4);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

GTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCGTACCAGGTA  >  NC_002947/3590188‑3590469                                                                                                                                               |                                                                                                                                            gttgGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGTGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGc                                                                                                                                        <  7:137796/148‑1 (MQ=255) gttgGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCggg                                                                                                                                         <  7:301552/147‑1 (MQ=255)  ttgGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGCCATGCCGGGca                                                                                                                                       <  2:177949/148‑1 (MQ=255)    ggCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGCCATGCCGGGcagt                                                                                                                                     <  2:257347/148‑3 (MQ=255)        tGAAGCACAGCGAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTc                                                                                                                                >  2:184927/1‑149 (MQ=255)        tGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGCCATGCCGGGcagt                                                                                                                                     >  4:305181/1‑142 (MQ=255)        tGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGCCATGCCGGGcagt                                                                                                                                     <  3:305181/144‑3 (MQ=255)             cacaGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGCCATGCCGGGcagt                                                                                                                                     <  1:188022/139‑3 (MQ=255)             cacaGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGCCATGCCGGGcagt                                                                                                                                     >  2:188022/1‑137 (MQ=255)              acaGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGcc                                                                                                                          >  1:308747/1‑149 (MQ=255)               caGCAAACCCGTCGGCGACTGGATACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCg                                                                                                                         >  2:196803/1‑149 (MQ=255)                     aCCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTAc                                                                                                                    >  8:282566/1‑148 (MQ=255)                        cGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAgg                                                                                                                <  3:72386/149‑1 (MQ=255)                                 cTGGATACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATg                                                                                                       <  1:196803/149‑1 (MQ=255)                                  tGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATgg                                                                                                      <  7:282566/149‑1 (MQ=255)                                    ggTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGc                                                                                                     <  6:277131/148‑1 (MQ=255)                                         tCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAgg                                                                                                                >  1:63984/1‑132 (MQ=255)                                         tCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAgg                                                                                                                <  2:63984/132‑1 (MQ=255)                                         tCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGc                                                                                               >  6:214720/1‑149 (MQ=255)                                           aGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGc                                                                                               <  7:132201/147‑1 (MQ=255)                                            gTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCg                                                                                            >  5:8927/1‑149 (MQ=255)                                            gTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCg                                                                                            >  5:84336/1‑149 (MQ=255)                                                 aaCCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTcc                                                                                        >  7:277613/1‑148 (MQ=255)                                                   ccTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCAc                                                                                     >  1:379794/1‑149 (MQ=255)                                                       ggTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACgg                                                                                   <  8:124689/147‑1 (MQ=255)                                                                  cgagcgTGAGACAGGCGAAGCGCAGCTGGGCAACATCGACAACCACGCGTTCTTCTCGCTGTTCACTGCCCGCCATGGCGGGCACAGATTCCATGCCGGCTACCACGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACa                                                                      <  6:243031/146‑1 (MQ=255)                                                                   gcgcgTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGATCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACAt                                                                     >  1:91931/1‑149 (MQ=255)                                                                            cAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCTGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGtt                                                            <  5:4232/149‑1 (MQ=255)                                                                            cAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGtt                                                            <  6:84336/149‑1 (MQ=255)                                                                             aaGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTg                                                           <  8:277613/149‑1 (MQ=255)                                                                              aGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTgg                                                          <  5:134963/149‑1 (MQ=255)                                                                               ggCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCAt                                                                                                            <  6:369675/98‑1 (MQ=255)                                                                               ggCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCAt                                                                                                            >  5:369675/1‑98 (MQ=255)                                                                                   aTGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAAc                                                     >  2:278677/1‑149 (MQ=255)                                                                                   aTGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAAc                                                     <  5:214720/149‑1 (MQ=255)                                                                                     gcgcTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAAc                                                     >  4:392886/1‑147 (MQ=255)                                                                                       gctgctGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGc                                                        <  3:290247/142‑1 (MQ=255)                                                                                       gctgctGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGc                                                        >  4:290247/1‑142 (MQ=255)                                                                                             gggCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGa                                           >  1:246204/1‑149 (MQ=255)                                                                                               gCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGt                                                                                 <  8:192650/109‑1 (MQ=255)                                                                                               gCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGt                                                                                 >  7:192650/1‑109 (MQ=255)                                                                                                 aaCATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTa                                       >  5:327080/1‑149 (MQ=255)                                                                                                  aCATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTAt                                      <  2:308747/149‑1 (MQ=255)                                                                                                    aTCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTAt                                      >  3:117163/1‑147 (MQ=255)                                                                                                          aaCCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGc                              <  1:278677/149‑1 (MQ=255)                                                                                                                gCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTacac                        <  3:392886/149‑1 (MQ=255)                                                                                                                        ctcGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGaa                 <  6:8927/148‑1 (MQ=255)                                                                                                                        ctcGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAAc                >  2:321271/1‑149 (MQ=255)                                                                                                                           gCTGTTCTCTGCCCGGCATGGCGGGCACAGCTTCCATGGCGGCTACGAGGGCATTTATGGCGTTAGCCCGTTGCCACGGGTATTTGCCATCATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTTCCTACACCGACGAACGCt             <  2:91931/149‑1 (MQ=255)                                                                                                                           gCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCt             <  1:321271/149‑1 (MQ=255)                                                                                                                            cTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTc            <  2:379794/149‑1 (MQ=255)                                                                                                                              gTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCg           >  3:3207/1‑148 (MQ=255)                                                                                                                               ttCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCGTa         >  4:343050/1‑149 (MQ=255)                                                                                                                                 cACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCGTAc        <  1:184927/148‑1 (MQ=255)                                                                                                                                   cTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCGTACCAg     <  5:343106/149‑1 (MQ=255)                                                                                                                                      cccGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCg                                                              >  7:324401/1‑89 (MQ=255)                                                                                                                                      cccGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCg                                                              <  8:324401/89‑1 (MQ=255)                                                                                                                                      cccGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCGTACCAGGTa  >  5:101602/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                            GTTGGCCTGAAGCACAGCAAACCCGTCGGCGACTGGGTACTCAGTGCCAACCTGGGTTATTTCACCGCGCGTGAGCAAGGCGATGCGCTGCTGGGCAACATCGACAACCAGGCGTTCTTCTCGCTGTTCACTGCCCGGCATGGCGGGCACAGCTTCCATGCCGGCTACCAGGGCATTTATGGCGATAGCCCGTTCCCACGGGTATTTGCCAACATCTCGCCGTTGGGCAACGAGGTGCCGACCTATGAGTTTGCCTACACCGACGAACGCTCGTACCAGGTA  >  NC_002947/3590188‑3590469

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