Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	146,128	T→A	14.0%	intergenic (+33/‑115)	gltP → / → PP_0138	glutamate/aspartate‑proton DAACS transporter/membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	146,128	0	T	A	14.0%	109.6 / 11.8	43	intergenic (+33/‑115)	gltP/PP_0138	glutamate/aspartate‑proton DAACS transporter/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base A (4/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01

TCATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTCGGCTTTTTCCTG  >  NC_002947/145983‑146275                                                                                                                                                  |                                                                                                                                                    tCATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCg                                                                                                                                                 >  1:221360/1‑150 (MQ=255)            tGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGt                                                                                                                                      <  1:333902/150‑1 (MQ=255)               gAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTAt                                                                                                                                   <  1:595921/150‑1 (MQ=255)                            acgCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTg                                                                                                                      >  1:606309/1‑150 (MQ=255)                                  tCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCcgc                                                                                                                <  2:606309/150‑1 (MQ=255)                                    aaGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCgcgc                                                                                                              >  2:715746/1‑150 (MQ=255)                                      gggCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGTGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCt                                                                                                            >  1:572984/1‑150 (MQ=255)                                      gggCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCggg                                                                                                                                               <  2:433760/115‑1 (MQ=255)                                      gggCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCggg                                                                                                                                               >  1:433760/1‑115 (MQ=255)                                       ggCGAGCCGTACTGCGCCTCGCCAATGTTAGAGAAGCAGGGAGCTGTATTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTa                                                                                                           <  2:536406/150‑1 (MQ=255)                                          gAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTCGTGGTTGGCGAAACCGCGAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTGGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACaa                                                                                                        >  1:368830/1‑150 (MQ=255)                                                  ctacGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGTGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTcgc                                                                                                >  1:141805/1‑150 (MQ=255)                                                        cTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGc                                                                                          <  2:243908/150‑1 (MQ=255)                                                        cTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGATTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGc                                                                                          >  1:623676/1‑150 (MQ=255)                                                                  gCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAgg                                                                                >  1:437448/1‑150 (MQ=255)                                                                      agaAGCAGGGAGCTGTAGTGGTTGGCGAAACCGGCTAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGc                                                                            <  2:623676/150‑1 (MQ=255)                                                                         ggCAGGGACCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCGTGATGGCTGTTTGACCATGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGGGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGAt                                                                         <  2:368830/149‑1 (MQ=255)                                                                             gggAGCTGTAGTGGTTGGCGAAACCGCCAAGTGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGc                                                                     <  2:572984/150‑1 (MQ=255)                                                                              ggAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAg                                                                                                      <  2:576749/116‑1 (MQ=255)                                                                              ggAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAg                                                                                                      >  1:576749/1‑116 (MQ=255)                                                                              ggAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAACCTGATGGCTGTGTGACCAAGAGCACCGACTTCTCGGGGCTTTTAGTTAACTGCGAGATATTGAGCCGCGCCTACAAAGGATCGAGATCGGCAGTCACAGGTTGCTCTACgac                                                                    >  2:128405/1‑148 (MQ=255)                                                                                   tGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGc                                                                                                    >  1:282506/1‑113 (MQ=255)                                                                                   tGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGc                                                                                                    <  2:282506/113‑1 (MQ=255)                                                                                     tAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTc                                                             >  1:420868/1‑150 (MQ=255)                                                                                          gTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCg                                                            >  1:14958/1‑146 (MQ=255)                                                                                          gTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCg                                                            <  2:14958/146‑1 (MQ=255)                                                                                               cGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCt                                                                    <  1:109894/133‑1 (MQ=255)                                                                                               cGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCt                                                                    >  2:109894/1‑133 (MQ=255)                                                                                               cGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCt                                                                    >  2:109753/1‑133 (MQ=255)                                                                                               cGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCt                                                                    <  1:109753/133‑1 (MQ=255)                                                                                                 aaagtaccAAGCGCTGACCCTGATGGCTGTTAGACCACGAGCCCCGACTTGCCTGGGATTGGAGTTATCTGCGAGAGCTTTATCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACACTCTCATTCTGGCACTTTTCa                                                 >  2:290154/7‑150 (MQ=255)                                                                                                   aCCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAgg                                               >  2:657031/1‑150 (MQ=255)                                                                                                                   ccTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGc                               >  1:363178/1‑150 (MQ=255)                                                                                                                               ttGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCttt                   <  2:646512/150‑1 (MQ=255)                                                                                                                                   ccAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCt               <  2:363178/150‑1 (MQ=255)                                                                                                                                    cAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTc              <  2:420868/150‑1 (MQ=255)                                                                                                                                    aaaaagCCCCGACAAGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGa                                  <  2:397371/126‑1 (MQ=255)                                                                                                                                    aaaaagCCCCGACAAGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGa                                  >  1:397371/5‑130 (MQ=255)                                                                                                                                    aaaaagCCCCGACAAGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTc              >  1:556426/5‑150 (MQ=255)                                                                                                                                          ccccGACAAGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACAt                                       >  2:343799/1‑119 (MQ=255)                                                                                                                                          ccccGACAAGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACAt                                       <  1:343799/119‑1 (MQ=255)                                                                                                                                           cccGACAAGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTCGGCtttt       >  1:443752/1‑150 (MQ=255)                                                                                                                                                cTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAGCATACGGATGCTCAACGGCCTTTGGCTCGGCTTTTTCCtg  <  1:715746/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TCATTGCTCGCTGGGAAGGCATGTACGACGCGGTCAAGGGCGAGCAGTACTACGCCTCGCTGATGGCCGAGAAGCAGGGAGCTGTAGTGGTTGGCGAAACCGCCAAGCGCTGAGCCTGATGGCTGTTTGACCAAGAGCCCCGACTTGTCGGGGCTTTTTGTTATCTGCGAGATCTTGAGCCGCGCCTACAAAGGCTCGCGATCGCCAGTCACAGGCTGCGATACGCTATTATTCGGGCACTTTTCAGGGGGAACATACGGATGCTCAACGGCCTTTGGCTCGGCTTTTTCCTG  >  NC_002947/145983‑146275

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