Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,295,031:1	+T	15.0%	intergenic (+41/+29)	slyB → / ← nhaA‑I	outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,295,031	1	.	T	15.0%	123.7 / 19.4	40	intergenic (+41/+29)	slyB/nhaA‑I	outer membrane lipoprotein/Na(+)/H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base . (12/22);  new base T (3/3);  total (15/25)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCTACCC  >  NC_002947/1294882‑1295177                                                                                                                                                       |                                                                                                                                                   tgGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGTGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCa                                                                                                                                                     >  2:474898/1‑150 (MQ=255)  ggCNGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑g                                                                                                                                                   >  2:701310/1‑150 (MQ=255)     agcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATGTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCg                                                                                                                                                >  1:41964/1‑150 (MQ=255)      gcaTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAACCCCCGACATGGCTGCTGCCA‑GTCgg                                                                                                                                               <  2:1123622/150‑1 (MQ=255)       caTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCggg                                                                                                                                              <  2:807678/150‑1 (MQ=255)        aTGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCgggg                                                                                                                                             >  1:188575/1‑150 (MQ=255)          gcgcgcCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGtt                                                                                                                                           <  1:145931/150‑1 (MQ=255)                cTANGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCa                                                                                                                                     >  2:702747/1‑150 (MQ=255)                 tATGTGCAGGCCGTGCAGGAGAACGAAATGTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCAt                                                                                                                                    <  2:41964/150‑1 (MQ=255)                      gCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTc                                                                                                                               >  2:35263/1‑150 (MQ=255)                      gCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTc                                                                                                                               >  1:507685/1‑150 (MQ=255)                        aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAg                                                                                                                             <  1:108319/150‑1 (MQ=255)                        aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCGGGGtttt                                                                                                                                         <  2:624351/138‑1 (MQ=38)                        aGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGCTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGAC‑TGGCAGCAGCCATGTCGGGGtttt                                                                                                                                         >  1:624351/1‑138 (MQ=38)                                  ggAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGAAGTAATAATAACCGCCAACATGACACCTCCCT‑GTCAGGATTTGTTCAGGCCTCAGCCAAGTCAgg                                                                                                                   <  2:852206/150‑1 (MQ=255)                                          aaaTCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTg                                                                                                           <  1:701310/150‑1 (MQ=255)                                          aaaTCTTCCGTGTTGGCGACCGAGTGCGAATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACACGGGCGCTTGCA‑GTCGGGTTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTg                                                                                                           <  1:1041467/150‑1 (MQ=255)                                                 ccGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGc                                                                                                    <  2:771968/150‑1 (MQ=255)                                                                   gcgcATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGa                                                                                  <  1:853249/150‑1 (MQ=255)                                                                     gcATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGACCATTGCGGTTAct                                                                                <  2:545450/150‑2 (MQ=255)                                                                            gACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGc                                                                         <  2:188575/150‑1 (MQ=255)                                                                                            aGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCaaa                                                         >  1:649650/1‑150 (MQ=255)                                                                                                   tttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaa                                                  <  2:314872/150‑1 (MQ=255)                                                                                                    ttCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAgaacgaac                                                 <  1:35263/150‑1 (MQ=255)                                                                                                      cGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCa                                               <  1:474898/150‑1 (MQ=255)                                                                                                       gCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAg                                              <  1:24699/150‑1 (MQ=255)                                                                                                           tGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgcag                                                             >  2:725489/1‑131 (MQ=255)                                                                                                           tGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCgcag                                                             <  1:725489/131‑1 (MQ=255)                                                                                                                gTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGc                                     <  2:649650/150‑1 (MQ=255)                                                                                                                         taataaAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAcc                                                                 >  2:109905/1‑113 (MQ=255)                                                                                                                         taataaAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCAcc                                                                 <  1:109905/113‑1 (MQ=255)                                                                                                                           ataaAAAACCCCGACATGGCTGCTGCCA‑GTCgggg                                                                                                                                             <  1:543343/35‑1 (MQ=255)                                                                                                                           ataaAAAACCCCGACATGGCTGCTGCCA‑GTCgggg                                                                                                                                             >  2:543343/1‑35 (MQ=255)                                                                                                                            taaAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCa                         <  2:507685/150‑1 (MQ=255)                                                                                                                               aaaaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGc                      >  2:725001/1‑150 (MQ=255)                                                                                                                                 aaCCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCg                    <  1:725001/150‑1 (MQ=255)                                                                                                                                  aCCCCGAC‑TGGCAGCAGCCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCgg                   >  1:482410/1‑150 (MQ=255)                                                                                                                                      cGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCACGCCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTa               <  2:760243/150‑1 (MQ=255)                                                                                                                                      cGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTa               >  1:40942/1‑150 (MQ=255)                                                                                                                                               tgctgcCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGa                                                                      <  1:917511/87‑1 (MQ=255)                                                                                                                                               tgctgcCATGTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGa                                                                      >  2:917511/1‑87 (MQ=255)                                                                                                                                                   gcCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCTGGGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTCAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCTAccc  >  2:252232/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                   TGGCAGCATGCGCGCCTATGTGCAGGCCGTGCAGGAGAACGAAATCTTCCGTGTTGGCGACCGAGTGCGCATCATGACTGTAGACGGTACCAGTCGCGTTTCGCACTGACCGTAGTGCAGTAATAAAAAACCCCGACATGGCTGCTGCCA‑GTCGGGGTTTTTTCATGCCTCAGCCAAGTCAGGCAATGCTGGTTTTGCGGCTGGCCATTGCGGTGACGGCGTAGCCGATCACCGCAGCAAAGAACGAACCAGTAAGAATGCCCATGCGGTCCATGCCGGCGTACTCACTGCTACCC  >  NC_002947/1294882‑1295177

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