Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,933,574 G→A 8.9% L39L (CTG→TTG)  fucA ← L‑fuculose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,933,5740GA 8.9% 135.8 / 3.5 45L39L (CTG→TTG) fucAL‑fuculose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (19/22);  new base A (2/2);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

CGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAA  >  NC_000913/2933434‑2933710
                                                                                                                                            |                                                                                                                                        
cGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGcatccc                                                                                                                                    <  3:114639/147‑1 (MQ=255)
      gCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGcatcccatc                                                                                                                                 >  5:145941/1‑144 (MQ=255)
            aTAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGc                                                                                                                                         >  4:46584/1‑130 (MQ=255)
            aTAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGc                                                                                                                                         <  3:46584/130‑1 (MQ=255)
               ggCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTac                                                                                                                     <  6:88177/147‑1 (MQ=255)
               ggCTGCAATATGGAAACGCCATTCGCTTGAGGGGATCTTTCCTTCATCATGTTTACCGTTGCCATAAATAAAGACAATATGCGAATCCGTCAGTTTTTGATATTGATTGCCTGTTGGCGTAATCAGCATCCCATCCTGATAACGTaca                                                                                                                    <  2:265813/148‑1 (MQ=255)
                        atGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAAc                                                                                                                         <  4:350022/134‑1 (MQ=255)
                        atGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAAc                                                                                                                         >  3:350022/1‑134 (MQ=255)
                                  cATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTcc                                                                                                 <  6:321008/148‑1 (MQ=255)
                                    ttCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGACATAAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCACATCATGATAAAGTACACTGACGTTCCCCGATGTccac                                                                                               >  3:137069/1‑146 (MQ=255)
                                     tCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTac                                                                                                                     >  3:116331/1‑125 (MQ=255)
                                     tCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTac                                                                                                                     <  4:116331/125‑1 (MQ=255)
                                       gCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTgg                                                                                            >  7:338746/1‑148 (MQ=255)
                                           gAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTc                                                                                         <  5:135306/147‑1 (MQ=255)
                                               ggAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTcc                                                                                    >  4:251003/1‑148 (MQ=255)
                                                gAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAACATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTcc                                                                                                 >  5:287941/1‑134 (MQ=255)
                                                gAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAACATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTcc                                                                                                 >  5:163165/1‑134 (MQ=255)
                                                gAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAACATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTcc                                                                                                 <  6:163165/134‑1 (MQ=255)
                                                gAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAACATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTcc                                                                                                 <  6:287941/134‑1 (MQ=255)
                                                      tccttcctCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCgg                                                                             >  4:335033/1‑148 (MQ=255)
                                                            ctCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCAtt                                                                       >  1:309287/1‑148 (MQ=255)
                                                              cATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTc                                                                     >  5:166734/1‑148 (MQ=255)
                                                               aTGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTcc                                                                    <  3:335033/148‑1 (MQ=255)
                                                               aTGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTc                                                                     >  2:48801/1‑147 (MQ=255)
                                                                tggttACCGTTTCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTGCCCCGCTGTCCCCTGGGTCAGTCCCAGGCGGGTCATTTCCa                                                                   <  4:137069/145‑1 (MQ=255)
                                                                 gTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAg                                                                  >  2:200029/1‑148 (MQ=255)
                                                                    tACCGTTGCCATCAATAAAGACAATATGCGCCTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGc                                                                >  4:315897/1‑147 (MQ=255)
                                                                    tACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGc                                                                >  8:355969/1‑147 (MQ=255)
                                                                             cATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCaat                                                      <  6:145941/148‑1 (MQ=255)
                                                                                cAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCaataat                                                   >  5:55643/1‑148 (MQ=255)
                                                                                 aaTAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATc                                                  >  5:294244/1‑148 (MQ=255)
                                                                                                  aCTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTAtt                                 <  1:48801/148‑1 (MQ=255)
                                                                                                         cAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTc                           >  3:318666/1‑147 (MQ=255)
                                                                                                          aGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTcc                          <  8:338746/147‑1 (MQ=255)
                                                                                                               ttCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGc                    >  4:225119/1‑148 (MQ=255)
                                                                                                                tCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGGTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGc                    <  6:55643/147‑1 (MQ=255)
                                                                                                                 cATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCt                   <  4:318666/147‑1 (MQ=255)
                                                                                                                 cATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCt                   <  1:200029/147‑1 (MQ=255)
                                                                                                                  atatGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTa                  <  3:225119/147‑1 (MQ=255)
                                                                                                                    atGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTTTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCt                   >  7:27636/1‑144 (MQ=255)
                                                                                                                    atGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTTTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCt                   <  8:27636/144‑1 (MQ=255)
                                                                                                                       gAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACctct             <  6:166734/147‑1 (MQ=255)
                                                                                                                        aaTGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACctctct           <  2:309287/148‑1 (MQ=255)
                                                                                                                            ccTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGGGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTg        <  3:251003/147‑1 (MQ=255)
                                                                                                                                  ggCTTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGTCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCaa  <  6:294244/147‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                        
CGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGTTTACCGTTGCCATCAATAAAGACAATATGCGACTCCGTCAGTTTTTCATATGGAATGCCTGTAGGCGTAATCAGCATCCCATCCTGATAACGTACACTGACGTTCCCCGCTGTCCCCTGGTTCAGTCCCAGGCGGGTCATTTCCAGGCAAGTGTCAATAATCTGACGAGCAAGTTTATTTCGTTCCATTAGCTACCTCTCTCTGATTCAA  >  NC_000913/2933434‑2933710

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