Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,442,209	G→A	9.6%	R108W (CGG→TGG) ‡	nicP‑II ←	porin‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,442,209	0	G	A	9.6%	137.6 / 4.3	52	R108W (CGG→TGG) ‡	nicP‑II	porin‑like protein
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base A (3/2);  total (29/23)
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.24e-01

AATCGCTGAACAGCAGCACCGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACTGGGTGCTGGC  >  NC_002947/4442062‑4442353                                                                                                                                                    |                                                                                                                                                 aaTCGCTGAACAGCAGCACCGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCgg                                                                                                                                                 >  7:202729/1‑149 (MQ=255)  atcccagAACACCAGCACCGCCAGGTTCGGCATCACCTCGCCCACTTTCAGCCCCGCTCGCGACACCTCGAACTTCATCGCCGCACCCCCCCGGCTGTACTCGCCCCCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCCGACCggg                                                                                                                                                <  8:119189/143‑1 (MQ=255)       tGAACTGCAGCACGGGTAGGTTTGGCATGAACTCGCCCTCTTTCAGCTCGGTTTGCGACACCTTGATCTTTTTCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCGGGACCGGGCAt                                                                                                                                             <  8:183403/147‑1 (MQ=255)       tGAACAGCAGCACCGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGt                                                                                                                                           >  6:184052/1‑149 (MQ=255)                    cGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAGCCCGGAACATGCCCGGTCCGGGttgctgc                                                                                                                              <  2:69455/149‑2 (MQ=38)                    cGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAGCCCGGAACATGCCCGGTCCGGGttgctgc                                                                                                                              <  6:429345/149‑2 (MQ=38)                    cGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGt                                                                                                                              <  8:411535/149‑1 (MQ=255)                      gCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTc                                                                                                                             >  4:317817/1‑148 (MQ=255)                          gTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCt                                                                                                                        >  7:164594/1‑149 (MQ=255)                                cATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAAc                                                                                                                  >  6:285748/1‑149 (MQ=255)                                    aaCTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCcag                                                                                                              <  8:88765/149‑1 (MQ=255)                                      cTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCcag                                                                                                              >  7:89132/1‑147 (MQ=255)                                           ccACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCg                                                                                                       >  2:79039/1‑149 (MQ=255)                                           ccACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCg                                                                                                       <  6:7676/149‑1 (MQ=255)                                               tttCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTaa                                                                                                   >  1:41939/1‑149 (MQ=255)                                                    gCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATc                                                                                               >  1:168505/1‑148 (MQ=255)                                                        ggTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCcac                                                                                            <  8:89132/147‑1 (MQ=255)                                                          tttGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACAccc                                                                                        >  8:275562/1‑149 (MQ=255)                                                              cGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCacac                                                                                          <  6:250563/143‑1 (MQ=255)                                                              cGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCacac                                                                                          >  5:250563/1‑143 (MQ=255)                                                               gACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTAGAAACCCAGCAGGCCGGTAACATCCACACCCACGc                                                                                    >  5:238844/1‑148 (MQ=255)                                                                   ccTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAgg                                                                               <  8:202729/149‑1 (MQ=255)                                                                   ccTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAgg                                                                               <  2:168505/149‑1 (MQ=255)                                                                   ccTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAgg                                                                               <  3:317817/149‑1 (MQ=255)                                                                    cTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAAcc                                                                                                                 <  8:212902/114‑1 (MQ=255)                                                                    cTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAAcc                                                                                                                 >  7:212902/1‑114 (MQ=255)                                                                       gAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGt                                                                             <  7:77302/147‑1 (MQ=255)                                                                          cTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGCCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCt                                                                        >  4:10633/1‑149 (MQ=255)                                                                             gATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTgg                                                                      <  5:285748/148‑1 (MQ=255)                                                                               tCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAATCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGta                                                                   >  1:49282/1‑149 (MQ=255)                                                                                  ccgcACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATaa                                                                >  2:135917/1‑149 (MQ=255)                                                                                     cACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCg                                                             <  6:267839/149‑1 (MQ=255)                                                                                       cccAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGa                                                           <  8:164594/149‑1 (MQ=255)                                                                                        ccAGGCGGCTGTATTCGTCCGCAGACTTGCCGCTGTCCAGGCTGGGAAGCAACCCGGACCGGGCATCTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAg                                                          >  2:162598/1‑149 (MQ=255)                                                                                              ggcTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGACTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCa                                                    >  5:204063/1‑149 (MQ=255)                                                                                                   tATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCCCGCCAAGGGTGACCTGGGTATAACCGGAGCTGGCATGGAg                                               >  4:291702/1‑149 (MQ=255)                                                                                                     ttCGTCCGCCGACTTGCCGCTGTCCAGGCTTGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAgg                                              <  8:173443/148‑1 (MQ=255)                                                                                                     ttCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAATCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAgg                                              <  2:49282/148‑1 (MQ=255)                                                                                                     ttCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAgg                                              <  7:275562/148‑1 (MQ=255)                                                                                                     ttCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACACACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGa                                             >  5:408172/1‑149 (MQ=255)                                                                                                           cgccgACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAAcc                                       >  8:208633/1‑149 (MQ=255)                                                                                                            gccgACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAAccc                                      >  5:29855/1‑149 (MQ=255)                                                                                                                aCTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTgc                                  <  6:408172/149‑1 (MQ=255)                                                                                                                      cGCCGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCAt                            >  4:339817/1‑149 (MQ=255)                                                                                                                        cTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTc                          >  8:408767/1‑149 (MQ=255)                                                                                                                         tGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGGCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATTCCCGGAGCTGGCAGTGAGGATAAACCCCTGCGCCCCTTcc                         >  4:195263/1‑149 (MQ=255)                                                                                                                             cAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTgg                      <  7:287477/148‑1 (MQ=255)                                                                                                                             cAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGCCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTgg                      <  3:10633/148‑1 (MQ=255)                                                                                                                                gctcgacagcagCCCGGAACATGCCCGGTCCGGGTTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTgc                  >  5:162961/7‑149 (MQ=37)                                                                                                                                   cgacagcagCCCGGAACATGCCCGGTCCGGGTTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGcgcg               >  8:343464/4‑149 (MQ=38)                                                                                                                                      cagcagCCCGGAACATGCCCGGTCCGGGTTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACt            >  8:188386/1‑149 (MQ=255)                                                                                                                                        gcaACCCGGACCGTGCATGTTCCTGGCTGCTGGCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACTgg          <  6:238844/149‑1 (MQ=255)                                                                                                                                             ccGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACTGGGTGCt     >  5:220714/1‑149 (MQ=255)                                                                                                                                                gACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACTGGGTGCTGgc  >  1:293644/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                 AATCGCTGAACAGCAGCACCGGCAGGTTGGGCATCAACTCGCCCACTTTCAGCTCGGTTTGCGACACCTTGAACTTGATCGCCGCACCCAGGCGGCTGTATTCGTCCGCCGACTTGCCGCTGTCCAGGCTGGGCAGCAACCCGGACCGGGCATGTTCCGGGCTGCTGTCGAGCTTGAAACCCAGCAGGCCGGTAACATCCACACCCACGCCAAGGGTGCCCTGGGTATAACCGGAGCTGGCATTGAGGATAAACCCCTGCGCCCATTCCTGGGTGCGCGACTGGGTGCTGGC  >  NC_002947/4442062‑4442353

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