Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,261,044 C→T 9.0% A123T (GCT→ACT)  fruA ← fused fructose‑specific PTS enzymes: IIBcomponent/IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,261,0440CT 9.0% 127.8 / 2.8 44A123T (GCT→ACT) fruAfused fructose‑specific PTS enzymes: IIBcomponent/IIC components
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base T (2/2);  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 = 5.77e-01

GATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGCTTCACTCAGGAACAGCTCAGG  >  NC_000913/2260897‑2261191
                                                                                                                                                   |                                                                                                                                                   
gATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc                                                                                                                                                    >  1:387159/1‑149 (MQ=255)
               cTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtgtgt                                                                                                                                     <  2:26378/149‑1 (MQ=255)
                 gccgccACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCGCCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGc                                                                                                                                   <  1:330048/149‑1 (MQ=255)
                    gccACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTAc                                                                                                                                <  1:342248/149‑1 (MQ=255)
                              ccAGATTCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCggg                                                                                                                      >  2:289104/1‑149 (MQ=255)
                                aGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCa                                                                                                                    >  2:61080/1‑149 (MQ=255)
                                 gATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCaa                                                                                                                   <  1:340618/149‑1 (MQ=255)
                                              gACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTAc                                                                                                                                >  1:398073/1‑123 (MQ=255)
                                              gACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTAc                                                                                                                                <  2:398073/123‑1 (MQ=255)
                                                     tcGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTAcgc                                                                                                >  3:9316/1‑148 (MQ=255)
                                                        ggAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGttt                                                                                            >  1:328657/1‑149 (MQ=255)
                                                             gATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGa                                                                                        >  3:174380/1‑148 (MQ=255)
                                                              aTTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGAcc                                                                                      >  2:265836/1‑149 (MQ=255)
                                                               ttGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGTCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCg                                                                                     >  2:98156/1‑149 (MQ=255)
                                                                tGCATTACCCGCACCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATAGCTTCAGTCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCCGACCg                                                                                     >  4:262677/1‑148 (MQ=255)
                                                                        ccGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCa                                                                            >  1:29977/1‑149 (MQ=255)
                                                                         cgcgCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCa                                                                            <  4:191686/148‑1 (MQ=255)
                                                                                 cAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGc                                                                    <  4:65113/148‑1 (MQ=255)
                                                                                 cAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGATTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTTCGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGc                                                                    >  3:59055/1‑148 (MQ=255)
                                                                                         cGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGtgtgtg                                                                                                                                    <  3:222942/76‑1 (MQ=255)
                                                                                         cGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGtgtgtg                                                                                                                                    >  4:222942/1‑76 (MQ=255)
                                                                                             tttCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGAc                                                        >  4:321673/1‑148 (MQ=255)
                                                                                                    tttCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTc                                                                                           <  1:201941/106‑1 (MQ=255)
                                                                                                    tttCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTc                                                                                           >  2:201941/1‑106 (MQ=255)
                                                                                                        acccaccAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTa                                            >  2:186672/1‑149 (MQ=255)
                                                                                                          ccaccaGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTaa                                           <  3:151760/148‑1 (MQ=255)
                                                                                                                   aCGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCa                                  >  3:79250/1‑148 (MQ=255)
                                                                                                                   aCGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCAt                                 >  1:354409/1‑149 (MQ=255)
                                                                                                                    cGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATg                                <  2:29977/149‑1 (MQ=255)
                                                                                                                      ttttttCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCg                                                      <  1:105680/125‑1 (MQ=255)
                                                                                                                      ttttttCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCg                                                      >  2:105680/1‑125 (MQ=255)
                                                                                                                      ttttttCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGc                              <  1:289104/149‑1 (MQ=255)
                                                                                                                           tCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTg                         <  2:33443/149‑1 (MQ=255)
                                                                                                                            cGCTTCGGTTTTAATGGCTTCAGTCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTg                         <  3:262677/148‑1 (MQ=255)
                                                                                                                            cGCTTCGGTTTCAATGGCTTCAGTCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTgg                        <  1:98156/149‑1 (MQ=255)
                                                                                                                            cGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTg                         <  4:271470/148‑1 (MQ=255)
                                                                                                                             gCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGc                       <  2:204688/149‑1 (MQ=255)
                                                                                                                               ttCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGCtt                     >  2:15575/1‑149 (MQ=255)
                                                                                                                               ttCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGCtt                     >  2:55238/1‑149 (MQ=255)
                                                                                                                                   gTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGCTTCACt                 >  1:381183/1‑149 (MQ=255)
                                                                                                                                     ttCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGCTTCACTCa               <  2:268117/149‑1 (MQ=255)
                                                                                                                                     ttCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGCTTCACTCa               <  1:306885/149‑1 (MQ=255)
                                                                                                                                         aTGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGCTTCACTCAGGaa           >  2:201555/1‑149 (MQ=255)
                                                                                                                                               tCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGCTTCACTCAGGAACAGCt      >  3:21650/1‑148 (MQ=255)
                                                                                                                                                cAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGCTTCACTCAGGAACAGCTc     <  4:9316/148‑1 (MQ=255)
                                                                                                                                                  gccgccATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGCTTCACTCAGGAACAGCTCAgg  <  2:381183/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGCTACGCCAGTCGGGCAAGCAGTCACCGCAACTACGCGTTTCGGACCGCTGGCGGCAACTGGTGCTGTCGCAGCGACCGGCGCAGTGTAAGGTTTCGCATGGCCTTTGGCTTCACTCAGGAACAGCTCAGG  >  NC_000913/2260897‑2261191

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