Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,937,463 C→A intergenic (+104/‑5) fucI → / → fucK L‑fucose isomerase/L‑fuculokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,937,4630CA92.3% 80.7 / 1.6 26intergenic (+104/‑5)fucI/fucKL‑fucose isomerase/L‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (11/13);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

CGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAA  >  NC_000913/2937338‑2937598
                                                                                                                             |                                                                                                                                       
cGGTCCGTCGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACaag                                                                                                                           >  1:369964/1‑140 (MQ=255)
    ccGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGt                                                                                                                       <  2:30610/140‑1 (MQ=255)
         gTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATcc                                                                                                                  <  1:154995/140‑1 (MQ=255)
         gTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATcc                                                                                                                  <  1:274063/140‑1 (MQ=255)
                cGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACt                                                                                                           <  2:115035/140‑1 (MQ=255)
                    aaTACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGGTAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGAc                                                                                                       >  1:302584/1‑140 (MQ=255)
                    aaTACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTAAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGAc                                                                                                       <  1:152731/140‑1 (MQ=255)
                           tACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGgcg                                                                                                >  1:203360/1‑140 (MQ=255)
                             cTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGcgcg                                                                                              >  2:334419/1‑140 (MQ=255)
                             cTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGCACTCGACTGTGGcgcg                                                                                              >  1:441745/1‑140 (MQ=255)
                                  tgatgtgatgCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCa                                                                                          <  2:334421/139‑1 (MQ=255)
                                  tgatgtgatgCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCa                                                                                          >  1:334421/1‑139 (MQ=255)
                                       tgatgCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCa                                                                                    >  1:240526/1‑140 (MQ=255)
                                            cccGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGAcc                                                                                           <  2:18810/128‑1 (MQ=255)
                                            cccGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGAcc                                                                                           >  1:18810/1‑128 (MQ=255)
                                              cGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCAt                                                                             <  2:203360/140‑1 (MQ=255)
                                               ggTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATc                                                                            >  1:296116/1‑140 (MQ=255)
                                                            aCCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGc                                                               <  2:192535/140‑1 (MQ=255)
                                                                  ccTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCa                                                         <  1:334419/140‑1 (MQ=255)
                                                                           gcaagAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAAttgttg                                                >  1:275565/1‑140 (MQ=255)
                                                                                   ggTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGcc                                        >  1:49479/1‑140 (MQ=255)
                                                                                       aGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCaa                                    <  2:240710/140‑1 (MQ=255)
                                                                                             aTAAGCAATGTTATCCGGCTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTa                              <  2:296116/140‑1 (MQ=255)
                                                                                                               cTATATTGCAGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGa            >  2:179646/1‑140 (MQ=255)
                                                                                                                        aGGAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAaaca   <  2:302584/140‑1 (MQ=255)
                                                                                                                         ggAGAGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAaacaa  <  2:118292/140‑1 (MQ=255)
                                                                                                                             |                                                                                                                                       
CGGTCCGTTGTACAAGCGTTAATACATTACTCCCTGATGTGATGCCCGGTCGCTCCGGCTACCGGGCCTGAACAAGCAAGAGTGGTTAGCCGGATAAGCAATGTTATCCGGCTATATTGCAGGAGCGATTATGAAACAAGAAGTTATCCTGGTACTCGACTGTGGCGCGACCAATGTCAGGGCCATCGCGGTTAATCGGCAGGGCAAAATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAA  >  NC_000913/2937338‑2937598

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: