Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,474,160	G→A	12.6%	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	12.6%	98.0 / 5.7	40	intergenic (‑65/+115)	PP_2167/tal	PhoD‑like phosphatase/transaldolase
Reads supporting (aligned to +/- strand):  ref base G (17/17);  new base A (3/2);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

GAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCTTT  >  NC_002947/2474012‑2474308                                                                                                                                                     |                                                                                                                                                     gAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACTGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAg                                                                                                                                                      <  5:155966/149‑1 (MQ=255)         cggccagacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGc                                                                                                                                              <  1:83581/149‑1 (MQ=255)               gacgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGccc                                                                                                                                        <  5:256391/149‑1 (MQ=255)                 cgCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCcgc                                                                                                                                      <  8:105840/149‑1 (MQ=255)                                                           aGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGcaaa                                                                                            <  2:116710/149‑1 (MQ=255)                                                                cAACGGGTTAGGGGAGGTCATCGACAGCTCGGGGGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAAt                                                                                                                 <  5:221509/123‑1 (MQ=255)                                                                      gTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCGGGCTTGCCcgc                                                                                                                                      <  8:182023/96‑1 (MQ=14)                                                                      gTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCTCCTACACTGATCGATGTAAGGGCGGGCTTGCCcgc                                                                                                                                      >  7:182023/1‑96 (MQ=14)                                                                             gAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTAGCGGGCAAGCCCGCCCGTACATCGATAAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAAATCAATAAACAAATAAAGCGTGGCAAAAAAAACGGGGGGaaaaa                                                                           >  4:243090/1‑147 (MQ=38)                                                                                   aTCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGaa                                                                     <  3:243090/148‑1 (MQ=255)                                                                                    tCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCg                                                                                                   <  7:157477/117‑1 (MQ=255)                                                                                       aCAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCg                                                                                 <  1:133129/132‑1 (MQ=255)                                                                                                     tATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGt                                                                                                  <  3:370821/101‑1 (MQ=255)                                                                                                     tATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAgg                                                   >  7:182197/1‑148 (MQ=255)                                                                                                           ggcCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAacac                                            >  1:12237/1‑149 (MQ=255)                                                                                                               tGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAAc                                        >  2:55460/1‑149 (MQ=255)                                                                                                               tGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGCCGTGGCAAACAAAACGGGCGGAAAATACTGAAAACTTACCGCCCGTAAGGCAAACACCAAc                                        >  8:140541/1‑149 (MQ=39)                                                                                                                gTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACt                                       >  4:269015/1‑149 (MQ=255)                                                                                                                 ttCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAAc                                        <  8:182197/147‑1 (MQ=255)                                                                                                                 ttCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTg                                      >  4:331041/1‑149 (MQ=255)                                                                                                                 ttCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGCAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAAc                                        <  8:43366/147‑1 (MQ=255)                                                                                                                  tCGCGGGAAAGCCCGCCCTTACATCGATAAGTGTATGAGCGGGCTTGCCCGCGAACAGGCCAGAGAATTGAATACACAAATAAAGCCTGGCAAACAAAACGCGCCGAAAATACTTAACACTTACCCCCCGTAAGGAAAAAACGAAcag                                      >  2:95412/1‑146 (MQ=255)                                                                                                                    gcgGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTcccgc                                                           >  4:223734/1‑125 (MQ=255)                                                                                                                    gcgGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTcccgc                                                           <  3:223734/125‑1 (MQ=255)                                                                                                                        gCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAAt                                                                                                                 >  3:319219/1‑67 (MQ=255)                                                                                                                        gCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAAt                                                                                                                 <  4:319219/67‑1 (MQ=255)                                                                                                                         cAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGc                              <  2:12237/149‑1 (MQ=255)                                                                                                                            gcccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCa                                 >  3:47690/1‑143 (MQ=37)                                                                                                                            gcccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCa                                 <  4:47690/143‑1 (MQ=37)                                                                                                                            gcccgctcCTACACTGATCGATGTAAGGGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtt                           >  3:131378/1‑149 (MQ=37)                                                                                                                               cgcccTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCtttt                         >  3:44678/1‑148 (MQ=255)                                                                                                                                    ttACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTa                   >  1:241178/1‑149 (MQ=255)                                                                                                                                        aTCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCa                       >  7:246680/1‑141 (MQ=255)                                                                                                                                        aTCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCa                       <  8:246680/141‑1 (MQ=255)                                                                                                                                        aTCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCaa               >  7:31084/1‑149 (MQ=255)                                                                                                                                        aTCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCa                >  6:254451/1‑148 (MQ=255)                                                                                                                                            aTCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACtt            >  2:84993/1‑148 (MQ=255)                                                                                                                                             tCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTtct          <  3:331041/149‑1 (MQ=255)                                                                                                                                                gtgtAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTtctc         >  3:57484/1‑147 (MQ=255)                                                                                                                                                gtgtAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTtctc         >  7:238358/1‑147 (MQ=255)                                                                                                                                                 tgtAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAg      <  5:254451/149‑1 (MQ=255)                                                                                                                                                     ggAGCGGGCTTGCCCGCGAACAGGCCCGATAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCttt  <  7:140541/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GAGCCAGACGGCCAGACGCTGGGGTTCCAGGCGGCGCAGTACCGGGCCGGCGAGAACGAGGGGCAACGGGTTAGGGGAGGTCATCGACAGCTCGGGTGGCTATGGCGGCCTGTTCGCGGGCAAGCCCGCCCTTACATCGATCAGTGTAGGAGCGGGCTTGCCCGCGAACAGGCCCGAGAATTCAATACACAACTAAGGCGTGGCAAACAAAACGGGCGGAAAATACTGAACACTTCCCGCCCGTCAGGCAAACACGAACTGTATCAGGCCTTTTCAGCTATCAACTTCTCAAGCTTT  >  NC_002947/2474012‑2474308

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