Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,458,117	T→G	15.7%	intergenic (+17/+23)	PP_1276 → / ← algA	hypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,458,117	0	T	G	15.7%	93.4 / 14.7	45	intergenic (+17/+23)	PP_1276/algA	hypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (15/23);  new base G (2/5);  total (17/28)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

AAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATC  >  NC_002947/1457969‑1458263                                                                                                                                                     |                                                                                                                                                   aaaCTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGt                                                                                                                                                    >  5:290667/1‑149 (MQ=255)      cGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGcc                                                                                                                                               >  3:175395/1‑149 (MQ=255)        aaaTGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGAGATGAATAAGCAGGCTCGCCCCTGGTCAGcc                                                                                                                                               >  7:276492/1‑147 (MQ=255)          aTGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAgg                                                                                                                                           <  6:290667/149‑1 (MQ=255)           tGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGGCTGACCAggg                                                                                                                                          <  1:219063/149‑1 (MQ=255)               tGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGa                                                                                                                                      <  1:344246/149‑1 (MQ=255)                 ccACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGa                                                                                                                                      <  8:36737/147‑1 (MQ=255)                 ccACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGGCTGACCAGGGGCGa                                                                                                                                      <  6:216848/147‑1 (MQ=255)                  cacaAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATACGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTc                                                                                                                                   <  1:120574/149‑1 (MQ=255)                   acaAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGt                                                                                                                                    >  8:11796/1‑147 (MQ=255)                    caAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAg                                                                                                                                 >  5:144821/1‑149 (MQ=255)                      aGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCa                                                                                                                               >  5:123783/1‑149 (MQ=255)                          tCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGGCTGACCAGGGGCGa                                                                                                                                      <  7:40269/138‑1 (MQ=255)                          tCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGGCTGACCAGGGGCGa                                                                                                                                      >  8:40269/1‑138 (MQ=255)                                    aTCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATg                                                                                                                 <  1:35290/149‑1 (MQ=255)                                       cTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGc                                                                                                                                        >  3:180966/1‑123 (MQ=255)                                       cTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGc                                                                                                                                        <  4:180966/123‑1 (MQ=255)                                        ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCt                                                                                                             >  7:152566/1‑149 (MQ=255)                                        ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCt                                                                                                             >  1:181411/1‑149 (MQ=255)                                          cGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCtt                                                                                                            <  7:11796/148‑1 (MQ=255)                                            aaTTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCAc                                                                                                         <  6:144821/149‑1 (MQ=255)                                             aTTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCAATTACTGCGCGATGGTCTTCAcc                                                                                                        >  5:163739/1‑149 (MQ=255)                                             aTTAGAAGGCGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATATACTGCGCCATGGTAGTCAcc                                                                                                        >  6:172616/1‑149 (MQ=255)                                                     gggAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACGTCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAAcgc                                                                                                 <  2:179455/148‑1 (MQ=255)                                                              atTTGCAGCGACCTTTCAGGGTGGTGTGTGTGCCATTCGACATCTACGTCATAGCCCGGAGCGTTTTATATGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGAt                                                                                       <  5:172616/149‑1 (MQ=255)                                                                 tGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGTGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATgg                                                                                     <  1:298670/148‑1 (MQ=255)                                                                   cAGCGACCTTTCAGGGTTGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATgggg                                                                                   <  2:301700/148‑1 (MQ=255)                                                                   cAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATggggg                                                                                  <  6:123783/149‑1 (MQ=255)                                                                   cAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATgggg                                                                                   <  6:332391/148‑1 (MQ=255)                                                                   cAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATgggg                                                                                   <  1:128938/148‑1 (MQ=255)                                                                        aCCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGAc                                                                             <  6:163739/149‑1 (MQ=255)                                                                                 tgttgtgggtggCACATTCGACATCTACGTCATAGCCCCGAGAGTTGAATAAGCAGGCTCGCCCCTGTTCAGCCCAGGGGCGAGTCAGCATTTACTGCTCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCcccc                                                                     <  8:276492/144‑1 (MQ=255)                                                                                          tggCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATcc                                                           <  8:152566/149‑1 (MQ=255)                                                                                              ccATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGa                                                       <  2:181411/149‑1 (MQ=255)                                                                                                                  aGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCg                                   >  7:131376/1‑149 (MQ=255)                                                                                                                       cGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAg                              >  7:155054/1‑149 (MQ=255)                                                                                                                                         acTCGCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTgcgc                                                                       <  2:47751/89‑1 (MQ=38)                                                                                                                                          cTCGCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTc                                                         <  2:59902/103‑1 (MQ=255)                                                                                                                                          cTCGCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTc                                                         >  1:59902/1‑103 (MQ=255)                                                                                                                                              ccccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCga       <  4:175395/149‑1 (MQ=255)                                                                                                                                               cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCTGGACTGCACCTCgat      >  2:122998/1‑149 (MQ=255)                                                                                                                                               cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTc         >  6:35093/1‑146 (MQ=255)                                                                                                                                               cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTc         <  5:35093/146‑1 (MQ=255)                                                                                                                                               cccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgat      >  6:60548/1‑149 (MQ=255)                                                                                                                                                ccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgatg     <  8:131376/149‑1 (MQ=255)                                                                                                                                                    gTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGAtc  <  5:60548/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   AAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATC  >  NC_002947/1457969‑1458263

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