Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,542,318	A→G	9.8%	Y135Y (TAT→TAC)	galE ←	UDP‑glucose 4‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,542,318	0	A	G	9.8%	136.6 / 5.3	51	Y135Y (TAT→TAC)	galE	UDP‑glucose 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base G (2/3);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

GATGGGTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGACGT  >  NC_002947/3542173‑3542463                                                                                                                                                  |                                                                                                                                                  gATGGGTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTat                                                                                                                                                <  2:250801/149‑2 (MQ=255)     ggTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACggg                                                                                                                                             <  6:288121/148‑1 (MQ=255)     ggTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGt                                                                                                                                            <  3:201136/149‑1 (MQ=255)       tttCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTg                                                                                                                                          <  5:45010/149‑1 (MQ=255)           cGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCg                                                                                                                                      <  4:59841/149‑1 (MQ=255)           cGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCg                                                                                                                                      <  7:298912/149‑1 (MQ=255)                     cGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgt                                                                                                                            <  1:207603/149‑1 (MQ=255)                      gCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgtg                                                                                                                           >  2:148228/1‑149 (MQ=255)                            cTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCt                                                                                                                     >  2:233044/1‑149 (MQ=255)                                   gCCTTCCGGGTAAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGCTGTCCTGAAGAATCTGCTCGGCCATCCATTTGCTCAACCCATACGGGGTGATCGGCCCCTTGGTGTGTGCGTCctcga                                                                                                              >  8:297946/1‑149 (MQ=255)                                     cTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcgatc                                                                                                            >  1:284304/1‑149 (MQ=255)                                      ttCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcgatcg                                                                                                           >  3:284673/1‑149 (MQ=255)                                      ttCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcgatcg                                                                                                           >  1:233140/1‑149 (MQ=255)                                           ggTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACg                                                                                                      <  2:111666/149‑1 (MQ=255)                                                 cgcCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGc                                                                                                >  7:71848/1‑149 (MQ=255)                                                  gcCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCg                                                                                               >  1:318605/1‑149 (MQ=255)                                                           aTTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGc                                                                                                         <  1:69992/130‑1 (MQ=255)                                                           aTTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGc                                                                                                         >  2:69992/1‑130 (MQ=255)                                                           aTTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCg                                                                                       >  7:301430/1‑148 (MQ=255)                                                               aaaTAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatc                                                                                                                >  8:227673/1‑119 (MQ=255)                                                               aaaTAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatc                                                                                                                <  7:227673/119‑1 (MQ=255)                                                               aaaTAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGtata                                                                                   <  1:136309/148‑1 (MQ=255)                                                                         ggCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTc                                                                        <  1:233044/149‑1 (MQ=255)                                                                            ttACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGTTGagc                                                                                                                                       >  3:230526/1‑81 (MQ=255)                                                                            ttACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGTTGagc                                                                                                                                       <  4:230526/83‑3 (MQ=255)                                                                            ttACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTCATCGGCCCCTTGGTGTGTGCTTCATCGATAGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCgaa                                                                     >  5:166436/1‑149 (MQ=255)                                                                              aCCGACTTCAACCCCTAGGCCCGGTGGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGTTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCgaa                                                                     <  3:52246/147‑1 (MQ=255)                                                                                    ttCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAacac                                                             >  8:121931/1‑149 (MQ=255)                                                                                        aaCCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCTgg                                                         >  7:213320/1‑149 (MQ=255)                                                                                              aGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCt                                                    >  8:15949/1‑148 (MQ=255)                                                                                               ggCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTg                                                   <  8:71848/148‑1 (MQ=255)                                                                                                       cGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCa                                          >  4:139229/1‑149 (MQ=255)                                                                                                        gAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCAt                                         >  5:276852/1‑149 (MQ=255)                                                                                                        gAAGTCCTCAAGAATCTGCTCGACCATACATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGCCACGTACTGCGGGTTGCCGTATACGGCGGCACTCCAAGAAAACACCAGGTGCCTGATGCCGGCAt                                         >  4:5303/1‑149 (MQ=255)                                                                                                          aGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTg                                       <  8:301430/149‑1 (MQ=255)                                                                                                            tCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGAc                                     <  7:297946/149‑1 (MQ=255)                                                                                                              cTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCa                                                           >  8:301233/1‑125 (MQ=255)                                                                                                              cTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCa                                                           <  7:301233/125‑1 (MQ=255)                                                                                                               tCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCAt                                  >  7:274169/1‑149 (MQ=255)                                                                                                                   gAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGcc                              >  2:303024/1‑149 (MQ=255)                                                                                                                     aTCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTg                            <  3:139229/149‑1 (MQ=255)                                                                                                                          cTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTgca                          <  7:121931/146‑1 (MQ=255)                                                                                                                           tCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGcagcag                      >  5:335135/1‑149 (MQ=255)                                                                                                                              aCCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAggg                    >  2:18644/1‑148 (MQ=255)                                                                                                                                catccatTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGaaa                 <  2:284304/149‑1 (MQ=255)                                                                                                                                catccatTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGaaa                 <  1:148228/149‑1 (MQ=255)                                                                                                                                     atTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTg            >  8:164596/1‑149 (MQ=255)                                                                                                                                           tCAACCCGTATGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTg        >  3:328680/1‑147 (MQ=255)                                                                                                                                           tCAACCCGTATGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTg        <  4:328680/147‑1 (MQ=255)                                                                                                                                           tCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGc       >  2:214879/1‑148 (MQ=255)                                                                                                                                               cccATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGACgt  >  1:185430/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GATGGGTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGACGT  >  NC_002947/3542173‑3542463

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