Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,307,957 A→G G174G (GGT→GGC)  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,307,9570AG97.6% 138.9 / ‑4.1 41G174G (GGT→GGC) alsED‑allulose‑6‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (22/18);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

GCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCG  >  NC_000913/4307835‑4308092
                                                                                                                          |                                                                                                                                       
gcgGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTc                                                                                                                        <  2:607233/140‑1 (MQ=255)
 cgGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCg                                                                                                                       <  1:802701/140‑1 (MQ=255)
   gTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTa                                                                                                                     <  1:821635/140‑1 (MQ=255)
       tAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTcc                                                                                                                 <  2:41793/140‑1 (MQ=255)
             tcCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCt                                                                                                           >  2:879508/1‑140 (MQ=255)
              ccATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCtt                                                                                                          >  2:574234/1‑140 (MQ=255)
                  gCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCg                                                                                                      >  1:60087/1‑140 (MQ=255)
                   cTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGt                                                                                                     >  2:888517/1‑140 (MQ=255)
                    ttCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGtt                                                                                                    >  1:922469/1‑140 (MQ=255)
                     tcgtcgATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTc                                                                                                   >  2:774261/1‑140 (MQ=255)
                       gtcgATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCAc                                                                                                 >  2:906377/1‑140 (MQ=255)
                       gtcgATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCAc                                                                                                 <  1:1131163/140‑1 (MQ=255)
                       gtcgATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCAc                                                                                                 <  2:60087/140‑1 (MQ=255)
                        tcgATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCtt                                                                                                          >  1:834403/1‑130 (MQ=255)
                        tcgATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCtt                                                                                                          <  2:834403/130‑1 (MQ=255)
                         cgATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGc                                                                                               <  2:922469/140‑1 (MQ=255)
                          catatTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACTTCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGcc                                                                                              <  2:158004/139‑1 (MQ=255)
                                     cATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGt                                                                                   >  2:235077/1‑140 (MQ=255)
                                         aTTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGTCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCgg                                                                               >  2:1000490/1‑140 (MQ=255)
                                           tAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCa                                                                             >  2:244258/1‑140 (MQ=255)
                                                  gCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTAt                                                                      >  2:849273/1‑140 (MQ=255)
                                                           gCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCAttt                                                             <  1:849273/140‑1 (MQ=255)
                                                            ccGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTc                                                            >  2:818530/1‑140 (MQ=255)
                                                            ccGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTc                                                            >  2:461558/1‑140 (MQ=255)
                                                            ccGACAATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTc                                                            >  2:956737/1‑140 (MQ=255)
                                                                   tAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATg                                                     >  2:1119898/1‑140 (MQ=255)
                                                                   tAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATg                                                     <  1:879508/140‑1 (MQ=255)
                                                                        aCATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACgg                                                <  1:818530/140‑1 (MQ=255)
                                                                           tCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTg                                             >  2:1079986/1‑140 (MQ=255)
                                                                           tCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTg                                             <  1:1062726/140‑1 (MQ=255)
                                                                             cGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTc                                           >  2:1068710/1‑140 (MQ=255)
                                                                                 ccTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGc                                       <  1:1029066/140‑1 (MQ=255)
                                                                                     ccgccATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAg                                   <  1:1068710/140‑1 (MQ=255)
                                                                                            cAGTTTTTCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGAt                            <  1:774261/140‑1 (MQ=255)
                                                                                               tttttCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGa                         <  1:574234/140‑1 (MQ=255)
                                                                                                   tCGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGt                     >  1:23666/1‑140 (MQ=255)
                                                                                                    cGTAAGTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTc                    >  2:94960/1‑140 (MQ=255)
                                                                                                         gTTGCCTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACACCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGAc               >  2:531148/1‑140 (MQ=255)
                                                                                                             ccTGGTTGCAGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTa           >  2:484468/1‑140 (MQ=255)
                                                                                                                      aGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCg  <  2:172969/140‑1 (MQ=255)
                                                                                                                      aGGAGCCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACACCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCg  <  1:531148/140‑1 (MQ=255)
                                                                                                                          |                                                                                                                                       
GCGGTCATAATTCTCCATGCTTCGTCGATATTTTCCGCATGATTAAACAGGCCGGAAGTGCCGACGATAAAGACATCCGCCCCTGCCGCCATCAGTTTTTCGTAAGTTGCCTGGTTGCAGGAACCGTCCACCTCAATTTCGTACTCCAGACCTTCTCGTTCACGCCATGCCTTCAGTTCGGCAAGTTTATCCAGCATTTCAGGAATGAACGGTTGTCCGGCAAAGCCGGGATCGACAGTCATGACCGTAATTTTATCG  >  NC_000913/4307835‑4308092

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