Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,068,332	C→G	9.8%	G325R (GGT→CGT) ‡	PP_4465 ←	porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,068,332	0	C	G	9.8%	140.7 / 6.3	51	G325R (GGT→CGT) ‡	PP_4465	porin
Reads supporting (aligned to +/- strand):  ref base C (22/24);  new base G (3/2);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

CCCGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCA  >  NC_002947/5068184‑5068477                                                                                                                                                     |                                                                                                                                                  cccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacc                                                                                                                                                   >  1:313257/1‑149 (MQ=255) cccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTAACCTTCTGGtcgacc                                                                                                                                                   >  6:156733/1‑149 (MQ=255) accGACTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACAATAACGTCGGTAGATAGATATAGGGATTGGAGGGCCCCTTCCCCTTATGGtcgccc                                                                                                                                                   >  6:263095/2‑149 (MQ=255)  ccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAATAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacct                                                                                                                                                  >  5:301744/1‑149 (MQ=255)  ccGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGtcgacct                                                                                                                                                  <  2:313257/149‑1 (MQ=255)           ggCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCAGCTTCTGGTCGACCTCGACCATc                                                                                                                                          <  6:88292/148‑1 (MQ=255)            gCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGAGGTCGACCagaag                                                                                                                                        <  5:364424/149‑5 (MQ=255)             ccTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAAc                                                                                                                                       >  7:179003/1‑149 (MQ=255)             ccTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAAc                                                                                                                                       <  2:246387/149‑1 (MQ=255)                   acacCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGtt                                                                                                                                 <  4:140465/149‑1 (MQ=255)                      ccTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGCCCTCGACCATCAGCAACGTTGc                                                                                                                               >  6:390491/1‑148 (MQ=255)                           ggcggcGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGtt                                                                                                                                 <  6:112162/141‑1 (MQ=255)                           ggcggcGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGtt                                                                                                                                 >  5:112162/1‑141 (MQ=255)                                 gAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACgg                                                                                                                   >  8:20525/1‑149 (MQ=255)                                        tAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCa                                                                                                              >  1:249845/1‑147 (MQ=255)                                             gTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGc                                                                                                          >  8:332134/1‑146 (MQ=255)                                              tAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCa                                                                                                        <  7:363175/147‑1 (MQ=255)                                               aaaGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGGCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGa                                                                                                      <  2:422942/148‑1 (MQ=255)                                                aaGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGAc                                                                                                     <  1:333311/148‑1 (MQ=255)                                                   gCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTa                                                                                                  >  8:132480/1‑148 (MQ=255)                                                   gCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGAATATAGATATAGGGATTGGACGCCCCCTTCACCTACTGGTCGACCACGACCATCAACCACGATGCTACTATTGACGGTAGCAAAGCCAGACTTaa                                                                                                 >  3:181741/1‑148 (MQ=255)                                                     cAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCa                                                                                               <  7:20525/149‑1 (MQ=255)                                                     cAAAAGTAGAATTTTCGCACGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTAtaa                                                                                               >  6:423199/1‑147 (MQ=255)                                                          gTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGAc                                                                                          <  7:55283/149‑1 (MQ=255)                                                                ttttCGCCCTCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACttgttg                                                                                    <  6:301744/149‑1 (MQ=255)                                                                ttttCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGCCTTATCACTGACttgtt                                                                                     >  8:364267/1‑148 (MQ=255)                                                                     gCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGCTGCAAAGCCAGACTTATCACTGACTTGTTGATGcc                                                                               <  6:173666/149‑1 (MQ=255)                                                                      cccGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGccc                                                                              >  8:259034/1‑149 (MQ=255)                                                                       ccGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGccc                                                                              >  7:264359/1‑148 (MQ=255)                                                                                tGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAgc                                                                    <  5:156733/149‑1 (MQ=255)                                                                                tGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAgc                                                                    <  6:4552/149‑1 (MQ=255)                                                                                  aaGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAgcgc                                                                  <  7:132480/149‑1 (MQ=255)                                                                                   aGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAgcgc                                                                  <  8:179003/148‑1 (MQ=255)                                                                                             cATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCg                                                       >  3:145987/1‑149 (MQ=255)                                                                                              aTAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAgc                                                                    >  4:350607/1‑135 (MQ=255)                                                                                              aTAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAgc                                                                    <  3:350607/135‑1 (MQ=255)                                                                                              aTAAAGTCGGTATATAGATTTAGGGATTGGACGCCCACTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGc                                                      <  8:379296/149‑1 (MQ=255)                                                                                                 aCGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCtt                                                    <  4:145987/148‑1 (MQ=255)                                                                                                   gTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTta                                                  >  4:195750/1‑148 (MQ=255)                                                                                                      ggTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAgg                                              >  5:359130/1‑149 (MQ=255)                                                                                                        tatataGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTg                                            <  7:259034/149‑1 (MQ=255)                                                                                                           atagatatagGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGa                                           >  7:404471/1‑147 (MQ=255)                                                                                                           atagatatagGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCAATGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGa                                           >  8:165244/1‑147 (MQ=255)                                                                                                             agatatagGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGc                                       >  3:11804/1‑149 (MQ=255)                                                                                                                       aTTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTTTAGGTGAAAGCCAAAGACCa                              <  5:243554/148‑1 (MQ=255)                                                                                                                        ttGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAgg                            <  8:404471/149‑1 (MQ=255)                                                                                                                             cGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGta                       <  7:279262/149‑1 (MQ=255)                                                                                                                                ccccTTCACCTTCTGGTCGACCTCGCCCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtatta                    >  7:121211/1‑149 (MQ=255)                                                                                                                                ccccTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGtatta                    >  5:75331/1‑149 (MQ=255)                                                                                                                                    ttCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTa                >  7:22549/1‑149 (MQ=255)                                                                                                                                       ggttgatGGTCGAGGTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCGAAGACCAGGTGGTATTATCTAttt             >  6:193661/7‑149 (MQ=255)                                                                                                                                        gttgatGGTCGAGGTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTTTCTAtttt            >  6:409110/6‑149 (MQ=255)                                                                                                                                         ttgatGGTCGAGGTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTc           <  5:409110/145‑1 (MQ=255)                                                                                                                                            cTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTccc         <  5:423199/148‑1 (MQ=255)                                                                                                                                            cTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTcc          <  7:332134/147‑1 (MQ=255)                                                                                                                                            atGGTCGAGGTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCa        >  8:222839/2‑149 (MQ=255)                                                                                                                                                 cgacctcgaccATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCAtttt    >  2:182755/1‑148 (MQ=255)                                                                                                                                                  gacctcgaccATCAACCCCGTTTCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCa  <  5:263095/149‑1 (MQ=255)                                                                                                                                                  gacctcgaccATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGAcc                               <  2:369214/120‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCCGCCTTTAGGCCTGGAACACCTAAGGCGGCGAAATCGTAGCTGTAAAGGCCAAAAGTAGAATTTTCGCCCGCGCGCGTGAAGCTATTCACCATAACGTCGGTATATAGATATAGGGATTGGACGCCCCCTTCACCTTCTGGTCGACCTCGACCATCAACCACGTTGCTACTATTGACGGTTGCAAAGCCAGACTTATCACTGACTTGTTGATGCCCTAAAGTCAGCGCATGACCACCCGCTTTATAGGTGAAAGCCAAAGACCAGGTGGTATTATCTATTTTCCCATTTTCA  >  NC_002947/5068184‑5068477

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