Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,458,117	T→G	8.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	8.7%	129.9 / 4.8	46	intergenic (+17/+23)	PP_1276/algA	hypothetical protein/mannose‑6‑phosphate isomerase/mannose‑1‑phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base G (2/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

AAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATCTC  >  NC_002947/1457969‑1458265                                                                                                                                                     |                                                                                                                                                     aaaCTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGt                                                                                                                                                      <  3:192683/149‑1 (MQ=255) aaaCTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTAGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGt                                                                                                                                                      >  1:265808/1‑149 (MQ=255)       gAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGcc                                                                                                                                                 >  8:267056/1‑148 (MQ=255)          aTGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAgg                                                                                                                                             <  3:103346/149‑1 (MQ=255)             aCTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGccc                                                                                                                                                >  1:180991/1‑143 (MQ=255)             aCTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGccc                                                                                                                                                <  2:180991/143‑1 (MQ=255)              cTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGGCTGACCAGGGGc                                                                                                                                          <  7:224677/148‑1 (MQ=255)               tGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCTGACTCGCCCCTGGGCTGACCAGGGGc                                                                                                                                          >  3:156909/1‑147 (MQ=255)                   acaAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCa                                                                                                                                    <  4:113182/149‑1 (MQ=255)                    caAGCATCAGGTCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAg                                                                                                                                   <  8:265098/149‑1 (MQ=255)                         aTCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCAttt                                                                                                                              <  8:189683/149‑1 (MQ=255)                          tCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTa                                                                                                                             <  7:267056/149‑1 (MQ=255)                            aGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTGGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACt                                                                                                                           >  8:248180/1‑149 (MQ=255)                            aGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACt                                                                                                                           >  8:34276/1‑149 (MQ=255)                              ggCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTgc                                                                                                                         >  5:193817/1‑149 (MQ=255)                               gCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCCGCGACCTTTCAGGGTGGTTGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTgcg                                                                                                                        <  2:265808/149‑1 (MQ=255)                                        ttCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCt                                                                                                               >  4:322321/1‑149 (MQ=255)                                           gAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCa                                                                                                            <  4:311163/149‑1 (MQ=255)                                            aaTTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTTGTCAGCCCAGGGGCGCGTCAGCAGTTACTGCGCGATGGACTTCAc                                                                                                           >  8:108635/1‑149 (MQ=255)                                                               tttGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATg                                                                                        >  5:134632/1‑149 (MQ=255)                                                                     gCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTg                                                                                  >  3:36520/1‑149 (MQ=255)                                                                       gACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGa                                                                                >  6:129629/1‑149 (MQ=255)                                                                           tttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATgggg                                                                                     <  8:219480/140‑1 (MQ=255)                                                                           tttCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATgggg                                                                                     >  7:219480/1‑140 (MQ=255)                                                                               aGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTgcgc                                                                         <  6:193817/148‑1 (MQ=255)                                                                                      tgggtggCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGACCATaca                                                                  >  7:135980/1‑148 (MQ=255)                                                                                       gggtggCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATAcaca                                                                >  4:132138/1‑149 (MQ=255)                                                                                               cATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                        >  4:60283/1‑149 (MQ=255)                                                                                               cATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGaa                                                        <  3:322321/149‑1 (MQ=255)                                                                                                            cGTCATAGCCCCGAGCGTTGAATAAGCAGGCTGGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGAt                                               <  7:248180/145‑1 (MQ=255)                                                                                                            cGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGAt                                               <  7:34276/145‑1 (MQ=255)                                                                                                            cGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATg                                              <  3:60283/146‑1 (MQ=255)                                                                                                                 tAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTc                                      <  6:134632/149‑1 (MQ=255)                                                                                                                       cGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCACCATTTACTGCGAGATCGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGAGATCGATGTCATCCTCGCCCAg                                >  7:131218/1‑149 (MQ=255)                                                                                                                                 aTAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCg                      <  4:36520/149‑1 (MQ=255)                                                                                                                                     gCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCgcgc                                                                                              >  4:37010/1‑73 (MQ=255)                                                                                                                                     gCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCgcgc                                                                                              <  3:37010/73‑1 (MQ=255)                                                                                                                                      tgggCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAATCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTg                 >  6:283504/3‑149 (MQ=255)                                                                                                                                         gCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATacac                                                                 >  7:194397/1‑98 (MQ=255)                                                                                                                                         gCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATacac                                                                 <  8:194397/98‑1 (MQ=255)                                                                                                                                          cTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCAcc             >  5:12898/1‑149 (MQ=255)                                                                                                                                          cTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCAcc             >  4:40562/1‑149 (MQ=255)                                                                                                                                           tCGCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTc                                                           >  8:228036/1‑102 (MQ=255)                                                                                                                                           tCGCCCCTGGGCTGACCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTc                                                           <  7:228036/102‑1 (MQ=255)                                                                                                                                             gCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTc           >  1:9003/1‑148 (MQ=255)                                                                                                                                              ccccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGAAGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCga         <  5:283504/149‑1 (MQ=255)                                                                                                                                               acatGGTCAGCCCTGGGGCGAGTCAGAATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgat        <  8:135980/146‑1 (MQ=255)                                                                                                                                                ccTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCgatg       <  6:12898/149‑1 (MQ=255)                                                                                                                                                     tCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGACTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGAtctc  <  7:108635/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AAACTCGAAATGACTGCCACAAGCATCAGGGCTGAATCCTTCGAATTTGAAGGGGAGGATGATTTGCAGCGACCTTTCAGGGTGGTGGGTGGCCCATTCGACATCTACGTCATAGCCCCGAGCGTTGAATAAGCAGGCTCGCCCCTGGTCAGCCCAGGGGCGAGTCAGCATTTACTGCGCGATGGTCTTCACCGAAACGCCGCGCTCGATGGGGGTGGACGTGCGCCCATACACATCCTCGAAGCGCTCGATGTCATCCTCGCCCAGGTAGCTCCCGGACTGCACCTCGATGATCTC  >  NC_002947/1457969‑1458265

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