Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,297,722 A→T 8.7% F192I (TTT→ATT)  alsE ← allulose‑6‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,297,7220AT 8.7% 126.0 / 5.4 46F192I (TTT→ATT) alsEallulose‑6‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (25/17);  new base T (2/2);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

CCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGG  >  CP009273/4297574‑4297863
                                                                                                                                                     |                                                                                                                                             
cccccATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATaaa                                                                                                                                               <  2:49367/149‑1 (MQ=255)
   ccATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGAc                                                                                                                                            >  1:251163/1‑149 (MQ=255)
    cATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCAGGAGGCTGTACCTCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACa                                                                                                                                           <  2:723241/149‑1 (MQ=255)
         ccACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGGAGTG‑CCGACGATAAAGACATCCGc                                                                                                                                      <  2:647565/149‑1 (MQ=255)
           aCGCCCGCTACGACGGTATGCTGTTTTTGAATGAGGCTGTGCCGCGCTTTTTGCAGCCAGACTCTGCGCGGTCATAATTCTCCATGCTTCGTCGAGATTTTCCGCATCGTTAACCGCGCCAGCAGTG‑CCGACGAGAAAGTCGGCCGccc                                                                                                                                    >  1:666022/1‑149 (MQ=255)
                cGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCGCTTTTTGCATCCATAATCTGCGCGGTCATAATTCTCCATCCTTCGTCGATATTTTCCGCATGATTAAACACGCCGCAAGTG‑CCGACGATCAAGACATCCGCCCCTgcc                                                                                                                               >  1:814540/1‑149 (MQ=255)
                              gCTGTTTTTGCATGAGGCTGTACCTCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTc                                                                                                                 >  1:676105/1‑149 (MQ=255)
                              gCTGTTTTTGCACGAGGCTGTTCCTCGCTTTTTGCATCCAGAACCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAATTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTc                                                                                                                 <  1:776043/149‑1 (MQ=255)
                               cTGTTTTTGCATGAGGCTGTACCTCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCGTAAGTTTTTCg                                                                                                                >  1:807018/1‑149 (MQ=255)
                               cTGTTTTTGCATGAGGCTGTACCTCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCCGCAGTTTTTCg                                                                                                                >  1:319750/1‑149 (MQ=255)
                                        cATGAGGCTGGACCTCGCTTTTTGCAGCCAGAATCTGCGCGGGCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGG‑AGTGACCGGCGGTCAAGACATCCGCCCCTGCCACCCGCAGCTTTTCGtgaattacc                                                                                                       >  1:245720/1‑141 (MQ=255)
                                             ggCTGTTCCTCGCTTTTTGCAGCCAGAATCTGCGCTGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGtt                                                                                                  <  1:700619/149‑1 (MQ=255)
                                                      tCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAATGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAAcc                                                                                         >  2:694235/1‑149 (MQ=255)
                                                      tCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAATGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAAcc                                                                                         >  2:274299/1‑149 (MQ=255)
                                                      tCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGATGCACAAAcc                                                                                         >  1:394814/1‑149 (MQ=255)
                                                      tCGCTTTTTGCAGCCAGAATCTGCGCCGTCATAACTCTCCATTCTTCGTCGATATTTGCCGCATGATTAACCAGGCCGGAAGTG‑CCGACGATAAAGTCATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCAGGTTGCAGGAAcc                                                                                         <  1:530050/149‑1 (MQ=255)
                                                      tCGCTGTTCGCAGCCAGAATCTGCGCGGTCATAATGCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAAcc                                                                                         <  1:609692/149‑1 (MQ=255)
                                                       cGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGCTATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGCCGATAAAGACATCCGCCCCTGCCGCCATCAGTGCTTCGTACGTTGCCTGGTTGAAGGAACCg                                                                                        >  1:554022/1‑149 (MQ=255)
                                                         cTTTTTGCGGCCAGACTCTGCGCGGTCATAATTCTCCCTGCTTCGTTGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCATGTCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTc                                                                                      <  1:267926/149‑1 (MQ=255)
                                                                  gCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAAt                                                                             >  2:823216/1‑149 (MQ=255)
                                                                     aGAATCTGCGCGGTCATACTTCCTCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACCTCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTc                                                                          <  2:815551/149‑1 (MQ=255)
                                                                               cgGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGa                                                                <  2:251163/149‑1 (MQ=255)
                                                                                  tCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCt                                                             >  1:649915/1‑149 (MQ=255)
                                                                                   cATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCtt                                                            >  1:455611/1‑149 (MQ=255)
                                                                                    aTAATTCTCCATGCTTCGTCGGTATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTAGCAGGAACCGTCCACCTCAATTTCGTACTCCCGACCTtc                                                           >  2:615367/1‑149 (MQ=255)
                                                                                                 gttCGTCGATATTTTGCGCATGATTAAACAGACGGGAAGTG‑CCGACGATAAAGACATGCGCAACTGCCGCGAAAAGTTTTTGGTAAGTTGCTGGGTGGTGGGAACGGTGTTCCTCAATTTCATATTCCAGACCTTCTCTTGCTCGCCAt                                              >  1:298347/2‑149 (MQ=255)
                                                                                                 cTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCAt                                              <  2:339459/149‑1 (MQ=255)
                                                                                                 cTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCAt                                              <  2:705716/149‑1 (MQ=255)
                                                                                                 cTTCGTCGAGATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGccct                                              >  1:281318/1‑147 (MQ=255)
                                                                                                      tcgATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCtt                                         >  2:14207/1‑149 (MQ=255)
                                                                                                        gATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACAGCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCGCGCCATGCCTTCa                                       >  2:622544/1‑149 (MQ=255)
                                                                                                            ttttCCGCATGATTACACAGGCCGGAAGTG‑CCGACGATAAATACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTTGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTc                                   <  1:458924/149‑1 (MQ=255)
                                                                                                                     tGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGttt                          >  2:349628/1‑149 (MQ=255)
                                                                                                                       aTTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTCCTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTAt                        >  2:719952/1‑149 (MQ=255)
                                                                                                                        ttAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATc                       >  2:475640/1‑149 (MQ=255)
                                                                                                                        ttAAACAGGCCGGAAGTG‑CCGACGATAAAGACAGCCGCCCCTGCCGCTATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATc                       >  2:561465/1‑149 (MQ=255)
                                                                                                                        atAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATc                       >  2:443649/2‑149 (MQ=255)
                                                                                                                                  ccgAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCa             <  1:781927/147‑1 (MQ=255)
                                                                                                                                  cGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCGGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCa             >  2:609468/1‑149 (MQ=255)
                                                                                                                                  cGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCa             <  2:729368/149‑1 (MQ=255)
                                                                                                                                   ggAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAg            <  1:823759/149‑1 (MQ=255)
                                                                                                                                      aGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTACGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGaa         <  2:814540/149‑1 (MQ=255)
                                                                                                                                            cgacggTAAAGCCCTCCTCCCCCGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAAc    <  1:253591/149‑1 (MQ=255)
                                                                                                                                            cgacgaTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAAc    >  1:794235/1‑149 (MQ=255)
                                                                                                                                             gacgaTAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACg   >  2:624403/1‑149 (MQ=255)
                                                                                                                                              acgaTAATGACATCCGCCCCTGCCGCCATCCGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACgg  <  1:274299/149‑1 (MQ=255)
                                                                                                                                              acgaTAATGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACgg  <  1:694235/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
CCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCGCTTTTTGCAGCCAGAATCTGCGCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTG‑CCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGG  >  CP009273/4297574‑4297863

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