Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,474,160	G→A	13.9%	intergenic (‑65/+115)	PP_2167 ← / ← tal	PhoD‑like phosphatase/transaldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,474,160	0	G	A	13.9%	82.7 / 5.8	36	intergenic (‑65/+115)	PP_2167/tal	PhoD‑like phosphatase/transaldolase
Reads supporting (aligned to +/- strand):  ref base G (14/17);  new base A (3/2);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

GAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCTTT  >  NC_002947/2474012‑2474308                                                                                                                                                      |                                                                                                                                                     gAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAg                                                                                                                                                      <  7:43744/149‑1 (MQ=255)     cagacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGc                                                                                                                                                  <  1:294043/149‑1 (MQ=255)      agacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGA‑TCGATGTAAGGGCg                                                                                                                                                 <  4:378701/149‑1 (MQ=37)       gacggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCgg                                                                                                                                                <  3:126041/149‑1 (MQ=255)         cggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCATCGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGA‑TCGATGTAAgggcgggc                                                                                                                                              <  1:83511/149‑1 (MQ=37)                                                 ggCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAg                                                                                                      <  1:345480/149‑1 (MQ=255)                                                 ggCGAGAACGAGGGGCAACGGGTCAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCTA‑TCAGTGTAGCAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAg                                                                                                      <  4:362800/149‑1 (MQ=255)                                                             gggCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTAGCGGGCAAGCCCGCCCGTTCCTCGA‑TGAGTGTAGGCGCGGGCTTGGCCGCGGACAGGCCCGCGAATTGAATACACAACTAAGGCGGGGgaaac                                                                                           >  1:109178/1‑148 (MQ=255)                                                                  aCGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACg                                                                                     <  3:188199/149‑1 (MQ=255)                                                                    gggTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACggg                                                                                   <  1:20165/149‑1 (MQ=255)                                                                               ggTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACt                                                                        <  1:309360/149‑1 (MQ=255)                                                                                          gCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTccc                                                             <  6:323823/149‑1 (MQ=255)                                                                                           cTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATAATGAACACTTcccg                                                            >  5:193981/1‑149 (MQ=255)                                                                                                 tGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTc                                                      <  7:103302/149‑1 (MQ=255)                                                                                                        ggcggcCGGTTCGCGGGCACGCCCCCCATAACATCTA‑ACAGTGTAGCCCCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCaaa                                               <  1:180739/149‑1 (MQ=255)                                                                                                          cggcCTGTTCGCGGGCAAGCCCGCCCGTACATCGC‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACCGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAACCAAAACGGGCGGGAAAGACTGACCACTTCCCGCCCGGCAGGCAAaca                                             >  5:384290/1‑149 (MQ=38)                                                                                                            gcCTGTTCGCGGGCAAGCCCGCCGTTACAGAGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAAGACTTCCAGCCCGTCAGGCAAaaac                                            >  6:317729/1‑148 (MQ=38)                                                                                                              cTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGaa                                         <  6:5215/149‑1 (MQ=255)                                                                                                              cTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGGGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGa                                          >  5:5215/1‑148 (MQ=255)                                                                                                               tGTTCGCGGGCAAGCCCGCCCTTACATCGAGTCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGaa                                         >  1:233111/1‑149 (MQ=255)                                                                                                                 ttCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGTAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACt                                       <  7:218647/148‑1 (MQ=255)                                                                                                                  tCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCATACACGAACTGt                                     >  4:253487/1‑149 (MQ=255)                                                                                                                      gggCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACATACCCGAGAATTCAATACACAACTAAGGCGGGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTAt                                   >  1:296640/1‑147 (MQ=255)                                                                                                                          aaGCCCGCCCTTACATCGA‑TCAGTGTAGGCGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAACCGGGCGGAAAATACTGCACACTTCCCGCCCGTCTGGCACACCCGAACTGTATCAGGcc                             >  4:4497/1‑149 (MQ=255)                                                                                                                            gcccgctcCTACACTGA‑TCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtt                           >  2:209458/1‑149 (MQ=37)                                                                                                                              ccgctcCTACACTGA‑TCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtttt                         >  1:246719/1‑149 (MQ=37)                                                                                                                              ccgctcCTACACTGA‑TCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCttt                          >  3:213482/1‑148 (MQ=37)                                                                                                                                    ttACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTa                   >  7:352824/1‑149 (MQ=255)                                                                                                                                        aTCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtttt                         <  1:298703/139‑1 (MQ=255)                                                                                                                                        aTCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATc                 <  7:347070/147‑1 (MQ=255)                                                                                                                                        aTCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATc                 >  8:347070/1‑147 (MQ=255)                                                                                                                                        aTCGA‑TCAGGGGAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAACCAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCGAACACGAACTGTATCAGGCCtgtg                         >  2:298703/1‑136 (MQ=255)                                                                                                                                         tCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAAc              >  8:310638/1‑149 (MQ=255)                                                                                                                                         atcgtccAGGGTCGGAGCGGCCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTACGGCGTGGCAAACACACCGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCC‑TTTCAGCTATCAAc              <  6:384290/143‑1 (MQ=255)                                                                                                                                                aGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCa        <  5:390731/149‑1 (MQ=255)                                                                                                                                                 gtgtAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTtctc         >  5:173377/1‑147 (MQ=255)                                                                                                                                                 gtgtAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCaa       <  3:366062/149‑1 (MQ=255)                                                                                                                                                      ggAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCttt  <  2:233111/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGA‑TCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCTTT  >  NC_002947/2474012‑2474308

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