Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,936,359	T→A	10.6%	V259D (GTT→GAT)	fucI →	L‑fucose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,936,359	0	T	A	10.6%	106.5 / 3.5	47	V259D (GTT→GAT)	fucI	L‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base T (14/25);  new base A (2/3);  total (17/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

CGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACT  >  NC_000913/2936226‑2936496                                                                                                                                      |                                                                                                                                          cGCCGTGTCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGAACTGcgcgg                                                                                                                                  <  2:177975/143‑2 (MQ=255) cGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGAACTGcgcgg                                                                                                                                  <  1:11297/143‑2 (MQ=255)   ccGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGaaa                                                                                                                                <  1:807655/143‑1 (MQ=255)     gTATCGATCAGAAGATTTACGACGAAGCCGAATTGTAAGTGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTCATGCCGAGCCACGCCGCGCACTCCTGCTCGACAGt                                                                                                                              >  1:494862/1‑143 (MQ=255)            tCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAAAAACAAACAGTATCAAAGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGa                                                                                                                       >  1:260480/1‑143 (MQ=255)             cAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGAt                                                                                                                      <  1:346623/143‑1 (MQ=255)                   aTTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGGTCTGCGCGAAAGTTTACGGATGGCggt                                                                                                                >  1:219382/1‑141 (MQ=255)                     ttACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGAtgt                                                                                                              <  2:767653/143‑1 (MQ=255)                     ttACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGAtgt                                                                                                              <  1:963645/143‑1 (MQ=255)                       acgacgAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATgtgt                                                                                                            >  1:1035738/1‑143 (MQ=255)                                ccGAATCGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGAc                                                                                                   <  2:460871/143‑1 (MQ=255)                                 cGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACa                                                                                                  >  1:602331/1‑143 (MQ=255)                                      tGGAAATGGCACTGGCCGCGGCTGATAAAAACTGCCGCTATGGCGAAGATGAAAATAACAAACAGTATCCACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatgatg                                                                                             <  2:845254/143‑1 (MQ=255)                                       ggAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGc                                                                                            <  2:894613/143‑1 (MQ=255)                                           aTGGCCCGGGCCTGGGCCGACAAAACCTTCCGCTATTGCGAAGATGAAACTACCATACAGTATCAACGTAATTCCGAGCAACGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCTACATGATGCAAgg                                                                                        <  2:600496/143‑1 (MQ=255)                                              gCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAAGCaa                                                                                     >  1:958358/1‑143 (MQ=255)                                              gCACTGGCCTGGGCGGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCaa                                                                                     <  1:424108/143‑1 (MQ=255)                                              gCACTGGCCTGGGCAGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGTTGTTTATCCGCGACATGATGCAAGGCaa                                                                                     <  2:129704/143‑1 (MQ=255)                                                     ccTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCaaa                                                                              <  2:674278/143‑1 (MQ=255)                                                         ggCAGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTgg                                                                          <  2:716880/143‑1 (MQ=255)                                                          gcggATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGc                                                                         <  2:755373/140‑1 (MQ=255)                                                                  aaaCTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTg                                                                 >  1:899788/1‑143 (MQ=255)                                                                      ttCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGATCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTcg                                                             >  1:692699/1‑143 (MQ=255)                                                                        gcGCTATGGCGGAAATGAAAGTAACAAGCAGTATCAACGTAATGCCGAACAAAGCCGCGCAGTTCTGCGGGAAAGTTTACTGATGGCGATGTGTAGCGGGGAGATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTcgcg                                                           <  2:989530/142‑1 (MQ=255)                                                                           cTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatgatg                                                                                             >  1:676040/1‑106 (MQ=255)                                                                           cTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACatgatg                                                                                             <  2:676040/106‑1 (MQ=255)                                                                            tATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGga                                                       >  2:404371/1‑143 (MQ=255)                                                                             aTGGCGAAGCTGAAAATACCCAACAGTACCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTCCTGATGGCGATGTGTATCCGCTACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaa                                                      <  1:643673/143‑1 (MQ=255)                                                                             aTGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaa                                                      <  1:783152/143‑1 (MQ=255)                                                                             aTGGAGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaa                                                      <  1:91113/143‑1 (MQ=255)                                                                                gCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGgaagaa                                                   >  1:833581/1‑143 (MQ=255)                                                                                 cGAAGATGAAAATAACATACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAAt                                                  <  2:203767/143‑1 (MQ=255)                                                                                 cGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCCCAGCAGCGAAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACCGCAAACTGGCCCACACTGAGCCCGTGgaggaaa                                                  >  2:73809/1‑142 (MQ=255)                                                                                        gAAACTAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTgg                                           <  2:30412/143‑1 (MQ=255)                                                                                        gAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTgg                                           >  1:985222/1‑143 (MQ=255)                                                                                        gAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTgg                                           >  2:1045702/1‑143 (MQ=255)                                                                                         aaaaTAACAACCAGTATCAACGTAAGGCCGAGCAAAGCCGCGCATTTCTGCGCGAAAGCTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGc                                          <  1:354780/143‑1 (MQ=255)                                                                                                 aaacaGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGCTGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGc                                  >  1:133815/1‑143 (MQ=255)                                                                                                     aGTATCAACGTAATGCCGAGCAAAGCCGCGCAGGTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAACATGGCCGCTATTTGTTGCGTGGAAGAATCACTTGGCTACAACGCCATc                              <  1:201941/143‑1 (MQ=255)                                                                                                           aaCGTAATCCCGAGCAAAGCCGCGCCGAACAGCCCGACAGCTTTCGCATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCg                        <  1:391318/143‑1 (MQ=255)                                                                                                                     cGACCAAAGCCGCCCAGCTCAGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCCACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAgg              <  1:30693/143‑1 (MQ=255)                                                                                                                      gggCAAAGGCGCGCAGTTCTGCGCGAAAGTTTACTGAGGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAggg             <  2:929883/141‑1 (MQ=255)                                                                                                                       aGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAgggg            >  1:934733/1‑143 (MQ=255)                                                                                                                          aaaGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCaa         <  1:89763/143‑1 (MQ=255)                                                                                                                             ttcgcgCAGAACTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACATGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGt      >  2:157313/3‑143 (MQ=255)                                                                                                                                gcgcAGTTCTGCGGGAAAGTTTACGGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCAc   <  2:493253/143‑1 (MQ=255)                                                                                                                                 cgaaGTTCTGGGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACt  <  2:695079/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                          CGCCGTATCGATCAGAAGATTTACGACGAAGCCGAATTGGAAATGGCACTGGCCTGGGCTGATAAAAACTTCCGCTATGGCGAAGATGAAAATAACAAACAGTATCAACGTAATGCCGAGCAAAGCCGCGCAGTTCTGCGCGAAAGTTTACTGATGGCGATGTGTATCCGCGACATGATGCAAGGCAACAGCAAACTGGCCGATATTGGTCGCGTGGAAGAATCACTTGGCTACAACGCCATCGCTGCGGGCTTCCAGGGGCAACGTCACT  >  NC_000913/2936226‑2936496

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