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_0009132,935,9170GA97.7% 134.9 / NA 44G112S (GGC→AGC) fucIL‑fucose isomerase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (22/21);  total (23/21)
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

AAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCG  >  NC_000913/2935801‑2936043
                                                                                                                    |                                                                                                                              
aaaaaTTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGt                                                                                                          <  2:831562/139‑1 (MQ=255)
                             aCCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTcaca                                                                             >  1:328134/1‑139 (MQ=255)
                               cATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGc                                                                           <  2:1034459/139‑1 (MQ=255)
                                      gTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAg                                                                    >  2:1312262/1‑139 (MQ=255)
                                       tAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg                                                                   <  1:1741802/139‑1 (MQ=255)
                                         aCGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGaa                                                                                                                        >  2:86911/1‑84 (MQ=255)
                                         aCGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGaa                                                                                                                        <  1:86911/84‑1 (MQ=255)
                                         aCGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgcg                                                                                                              >  2:1779798/1‑94 (MQ=255)
                                         aCGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgcg                                                                                                              <  1:1779798/94‑1 (MQ=255)
                                         aCGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCa                                                                 >  1:1126234/1‑139 (MQ=255)
                                                gCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGc                                                          >  2:1083014/1‑139 (MQ=255)
                                                 cTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCa                                                         <  1:1083014/139‑1 (MQ=255)
                                                 cTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCa                                                         >  1:1505612/1‑139 (MQ=255)
                                                    gTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATtc                                                      <  1:1312262/139‑1 (MQ=255)
                                                      gCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATtctc                                                    <  2:1505612/139‑1 (MQ=255)
                                                           ggCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgcg                                                                                                              <  1:209406/76‑1 (MQ=255)
                                                           ggCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgcg                                                                                                              <  1:1084312/76‑1 (MQ=255)
                                                           ggCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgcg                                                                                                              >  2:209406/1‑76 (MQ=255)
                                                           ggCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgcg                                                                                                              >  2:1084312/1‑76 (MQ=255)
                                                                tGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGt                                          >  1:266416/1‑139 (MQ=255)
                                                                   aaCCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCTCTGaa                                                                                                                        >  2:1212744/1‑58 (MQ=255)
                                                                   aaCCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGaa                                                                                                                        <  1:1212744/58‑1 (MQ=255)
                                                                         cGACATGGATCCAACCCGCCCGAAGGCCACTTCGGGCTTTAACGGCACTGAACGCGACGGAGCTGTCTACCTGGCAGCGGTTCTGGCAGCTCACAGCCAGAAGCGCATCCCAGCATGCCCCATTTACGGTCATGACGtt                                 >  1:437366/1‑139 (MQ=255)
                                                                             aTGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAgg                             <  2:328134/139‑1 (MQ=255)
                                                                             aTGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAgg                             <  1:996892/139‑1 (MQ=255)
                                                                                   ccAACCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCg                       >  2:895137/1‑139 (MQ=255)
                                                                                      aCCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg                                                                   >  2:1316185/1‑92 (MQ=255)
                                                                                      aCCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAgg                                                                   <  1:1316185/92‑1 (MQ=255)
                                                                                      aCCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATg                    >  1:733397/1‑139 (MQ=255)
                                                                                      aCCCGCCCGAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATg                    >  1:267755/1‑139 (MQ=255)
                                                                                           cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgc                                                                                                               >  2:1860/1‑43 (MQ=255)
                                                                                           cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGgc                                                                                                               <  1:1860/43‑1 (MQ=255)
                                                                                           cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTcac                                                                              >  2:1004187/1‑76 (MQ=255)
                                                                                           cccgAAGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTcac                                                                              <  1:1004187/76‑1 (MQ=255)
                                                                                                aGGCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGAt          >  1:567669/1‑139 (MQ=255)
                                                                                                  gCCATTTGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTc        <  2:567669/139‑1 (MQ=255)
                                                                                                      tttGGGGCTTTAACAGCACTGAACGCCCCGGCGCCGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGc    >  1:127572/1‑139 (MQ=255)
                                                                                                        tGGGGCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCg  <  2:733397/139‑1 (MQ=255)
                                                                                                            gCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCa                                                         >  2:1750548/1‑80 (MQ=255)
                                                                                                            gCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCa                                                         <  1:1750548/80‑1 (MQ=255)
                                                                                                            gCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGa                   <  1:739786/118‑1 (MQ=255)
                                                                                                            gCTTTAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGa                   >  2:739786/1‑118 (MQ=255)
                                                                                                               ttAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGtt                                 <  2:1146203/101‑1 (MQ=255)
                                                                                                               ttAACAGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGtt                                 >  1:1146203/1‑101 (MQ=255)
                                                                                                                    |                                                                                                                              
AAAAATTCAGCAGTCAGAATGTAGGCCTCACCATTACGGTAACGCCTTGCTGGTGCTATGGCAGTGAAACCATCGACATGGATCCAACCCGCCCGAAGGCCATTTGGGGCTTTAACGGCACTGAACGCCCCGGCGCTGTTTACCTGGCAGCGGCTCTGGCAGCTCACAGCCAGAAAGGCATCCCAGCATTCTCCATTTACGGTCATGACGTTCAGGATGCCGATGACACATCGATTCCTGCCG  >  NC_000913/2935801‑2936043

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