Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,935,917	G→A	G112S (GGC→AGC)	fucI →	L‑fucose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,935,917	0	G	A	97.7%	137.0 / NA	44	G112S (GGC→AGC)	fucI	L‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/22);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATGTTG  >  NC_000913/2935780‑2936049                                                                                                                                          |                                                                                                                                     cTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTGACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGATGGCCATTTGGGGCTTTAACgg                                                                                                                                     <  1:967597/139‑1 (MQ=255)    aaGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCAc                                                                                                                                  >  1:1453348/1‑139 (MQ=255)     aGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCATGCTGGTGCTATGGCAGTGAAACCATCGACAGGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACt                                                                                                                                 >  1:1046325/1‑139 (MQ=255)        cGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGaa                                                                                                                              <  2:639771/139‑1 (MQ=255)             ctTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGccc                                                                                                                         <  2:1453348/139‑1 (MQ=255)                                             ggCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAg                                                                                         >  1:778796/1‑139 (MQ=255)                                                 tCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTca                                                                                     >  1:1280047/1‑139 (MQ=255)                                                  cACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTcac                                                                                    <  2:839643/139‑1 (MQ=255)                                                   aCCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTcaca                                                                                   <  1:1006369/139‑1 (MQ=255)                                                   aCCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTca                                                                                     <  2:50086/137‑1 (MQ=255)                                                   aCCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTca                                                                                     >  1:50086/1‑137 (MQ=255)                                                        tACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGAACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAg                                                                              >  2:875596/1‑139 (MQ=255)                                                             tAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg                                                                         <  1:1427056/139‑1 (MQ=255)                                                             tAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg                                                                         <  2:1305323/139‑1 (MQ=255)                                                               aCGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCa                                                                       >  2:1553108/1‑139 (MQ=255)                                                               aCGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCa                                                                       >  1:769701/1‑139 (MQ=255)                                                                     tGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTcac                                                                                    <  1:1200675/120‑1 (MQ=255)                                                                     tGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTcac                                                                                    >  2:1200675/1‑120 (MQ=255)                                                                       cTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCa                                                               <  1:1553108/139‑1 (MQ=255)                                                                           tGCTATGGCAGTGAAACCATCGACATGGATCCAAACCGCCAGAAGGCAGTTTGGGGCGGTAACGGCACTGAACGCCCCGAGGCTGTTTACCTGGCAGCGGCTGTGGCAGCTCACAGCCAGAAAGGCATCCCAGCAttat                                                           >  1:1192841/1‑137 (MQ=255)                                                                                     gTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACgg                                                 >  1:195031/1‑139 (MQ=255)                                                                                      tGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCAtt                                                             <  2:397918/126‑1 (MQ=255)                                                                                      tGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCAtt                                                             >  1:397918/1‑126 (MQ=255)                                                                                       gAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGt                                                                                                                <  1:162536/74‑1 (MQ=255)                                                                                       gAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGt                                                                                                                >  2:162536/1‑74 (MQ=255)                                                                                         aaCCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCAt                                             <  1:280955/139‑1 (MQ=255)                                                                                            cATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGAc                                          >  2:1472941/1‑139 (MQ=255)                                                                                                     ggATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg                                                                         >  2:1134297/1‑99 (MQ=255)                                                                                                     ggATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg                                                                         <  1:1134297/99‑1 (MQ=255)                                                                                                          cAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGa                            >  1:1443416/1‑139 (MQ=255)                                                                                                           aaCCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGAt                           <  2:1368338/139‑1 (MQ=255)                                                                                                            aCCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATg                          >  2:770430/1‑139 (MQ=255)                                                                                                             cccgcccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgcg                                                                                                                    <  1:1070309/48‑1 (MQ=255)                                                                                                             cccgcccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgcg                                                                                                                    >  2:1070309/1‑48 (MQ=255)                                                                                                                 cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATc                                                                     <  1:38704/91‑1 (MQ=255)                                                                                                                 cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATc                                                                     >  2:38704/1‑91 (MQ=255)                                                                                                                 cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCa                                                               >  2:1603386/1‑97 (MQ=255)                                                                                                                 cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCa                                                               <  1:1603386/97‑1 (MQ=255)                                                                                                                 cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAgg                                   <  2:1023999/125‑1 (MQ=255)                                                                                                                 cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAgg                                   >  1:1023999/1‑125 (MQ=255)                                                                                                                 cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGAcaca                     >  2:677930/1‑139 (MQ=255)                                                                                                                 cccgAAGGCCACTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGAcaca                     <  2:1046325/139‑1 (MQ=255)                                                                                                                                  gCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATGt    <  2:1443416/139‑1 (MQ=255)                                                                                                                                   cTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATGtt   <  1:677930/139‑1 (MQ=255)                                                                                                                                    tttAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATGTTg  <  2:769701/139‑1 (MQ=255)                                                                                                                                          |                                                                                                                                     CTGAAGCCGCTGCTTGCGAAGAAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCGATGTTG  >  NC_000913/2935780‑2936049

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