Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,308,169 G→A R104C (CGT→TGT)  alsE ← D‑allulose‑6‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,308,1690GA90.0% 21.3 / NA 10R104C (CGT→TGT) alsED‑allulose‑6‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (3/6);  minor base C (1/0);  total (4/6)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACC  >  NC_000913/4307946‑4308373
                                                                                                                                                                                                                               |                                                                                                                                                                                                            
tGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACAGCGGATTTCATCAATCAGGCGGAACGc                                                                                                                                                                                   >  1:8289/1‑251 (MQ=255)
  gTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACAGCGGATTTCATCAATCAGGCGGAACGCCt                                                                                                                                                                                 >  1:2214/1‑251 (MQ=255)
                                                    gCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACAGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCAcgn                                                                                                                               <  1:55829/251‑2 (MQ=255)
                                                                      cAAGTTTATCCAGCATTTCAGGAATGNACGGTTNNCCNGCAAAGCCGGGATCGACAGNCANNNCNNNANTTTTANCGNNCTTATGGATNTAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTNATGTCATGANNGCGGCTTtcatca                                                                                                                                                                                                 >  2:49145/1‑167 (MQ=255)
                                                                                          ggAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACAGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCAt                                                                                         >  1:168577/1‑251 (MQ=255)
                                                                                                       gTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACAGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATg                                                                                                                <  1:65774/215‑1 (MQ=255)
                                                                                                          cGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACAGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAgn                                                                         <  1:43451/251‑2 (MQ=255)
                                                                                                                              gTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTACATGGCCTAAAATGACGTCCCCGAATTACGGATCAACCCCCCTTTCTTGTCATGACCACGGATTTCATAAATCAGGAGGCACGCCTGGACGTTGATgg                                                                                                                                                                     >  1:144825/1‑139 (MQ=255)
                                                                                                                                                                        cnatncTCAANTGGCNCGTGCGGGTTAAGGANCAGCCCCACTTTCATGTCANNANAGCGGNTTtca                                                                                                                                                                                                    <  2:86204/61‑1 (MQ=255)
                                                                                                                                                                                 aCTGGCGTCTCCGGGGTAAGGATCAGCCCCACTTTCATGGCATGACAGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACTCGCCAGGTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAGCCTGACTAAATAAGAACGGTGAGAGTGGCAGATTGGGGACAAAGTGAcc  <  1:188717/251‑1 (MQ=255)
                                                                                                                                                                                                  anNGATCAGNNCCACTTTCATGTCATGACAGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTGTCCGGATGCAGAGTGATGAAATCTGc                                                                                                                                         <  2:2214/99‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                            
TGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCGGCCTTATGGATATAGTATTTCATGGCCTCAACTGGCGTCTCCGGGTTAAGGATCAGCCCCACTTTCATGTCATGACGGCGGATTTCATCAATCAGGCGGAACGCCTGGCCGTTGATGGTTTCCGGATGCAGAGTGATGAAATCTGCTCCCGCACGCGCCAGTTGAGCAATGTAATCCTGCGGCCGCGTCACCATCAGATGACAGTCGAGCGGTTTAGTTGCCAGTTTTTTAACCTGACTTACGAAGAACGGTGAGAGTGTCAGATTGGGGACAAAGTGACC  >  NC_000913/4307946‑4308373

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