Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,423,375 A→G N177S (AAT→AGT)  ulaE → L‑ribulose‑5‑phosphate 3‑epimerase UlaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,423,3750AG95.7% 144.7 / ‑2.6 47N177S (AAT→AGT) ulaEL‑ribulose‑5‑phosphate 3‑epimerase UlaE
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (24/21);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GTCGTTTCCGTGACGGCCTGAAAGAGAGCGTTGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAATCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACATGTGAAAGACACCAAACCTGGCGTCTTCA  >  NC_000913/4423240‑4423506
                                                                                                                                       |                                                                                                                                   
gtcgtTTCCGTGACGGCCTGAAAGAGAGCGTTGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCg                                                                                                                                 <  2:804244/140‑1 (MQ=255)
gtcgtTTCCGTGACGGCCTGAAAGAGAGCGTTGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCg                                                                                                                                 <  2:690286/140‑1 (MQ=255)
         gTGACGGCCTGAAAGAGAGCGTTGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAg                                                                                                                        >  1:750852/1‑140 (MQ=255)
                          agCGTTGAGATGGCAAGCCGCGTGCTGGTTACGCTGGTGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCgg                                                                                                       <  1:306672/140‑1 (MQ=255)
                          agCGTTGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATc                                                                                                         >  2:407446/1‑138 (MQ=255)
                          agCGTTGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATc                                                                                                         <  1:407446/138‑1 (MQ=255)
                          agCGTTGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTAAACAGTCCGTGGTTCCAACTCTACCCGGATATCgg                                                                                                       >  2:309982/1‑140 (MQ=255)
                            cGTTGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCa                                                                                                     >  1:132271/1‑140 (MQ=255)
                             gTTGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCaa                                                                                                    >  2:445488/1‑140 (MQ=255)
                              ttGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAAc                                                                                                   >  1:808397/1‑140 (MQ=255)
                              ttGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCACACTATCTTAACAGTCTGTGGTTCCAGCTCTACCCGGATATCGGCAAc                                                                                                   >  2:306672/1‑140 (MQ=255)
                                      gCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGc                                                                                           >  2:586395/1‑140 (MQ=255)
                                       cAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCg                                                                                          <  2:132271/140‑1 (MQ=255)
                                        aaGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGt                                                                                         <  1:445488/140‑1 (MQ=255)
                                            cgcgcgCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGa                                                                                     <  2:750852/140‑1 (MQ=255)
                                              cgcgcAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACa                                                                                   <  2:600990/140‑1 (MQ=255)
                                              cgcgcAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCACTAGCCGGATATCGGCAAGCTGTCGGGGTGGGACa                                                                                   >  2:393054/1‑140 (MQ=255)
                                                   aGGTGACGCTGGCGACGGAGATCATGGATTATCCGTTAATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGAt                                                                              <  1:372406/140‑1 (MQ=255)
                                                                              aTTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATc                                                   >  1:527537/1‑140 (MQ=255)
                                                                              aTTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCATGGTAGATCGATGTGCAGATGTAGTTGCAGGCCGGAATc                                                   >  2:370681/1‑140 (MQ=255)
                                                                                 aTCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCggg                                                >  1:386/1‑140 (MQ=255)
                                                                                  tCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGGCGGCGTGGGACAGCGAGGTGCAGATTGAGTGGCAGGCCGGAATCGGGc                                               >  2:223644/1‑140 (MQ=255)
                                                                                       tgatgaGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGc                                                          <  2:586462/124‑1 (MQ=255)
                                                                                       tgatgaGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGc                                                          >  1:586462/1‑124 (MQ=255)
                                                                                         atgaGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAAt                                                    <  1:436817/128‑1 (MQ=255)
                                                                                         atgaGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAAt                                                    >  2:436817/1‑128 (MQ=255)
                                                                                         atgaGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATAtcgt                                        >  2:551701/1‑140 (MQ=255)
                                                                                           gaGCTCCATCAGCAATGCGCTGTGATACGAGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATAtcgtcg                                      <  1:223644/140‑1 (MQ=255)
                                                                                           gaGCTCCATCAGCAATGCGCTGTGAAACGCGCACTATTTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATAtcgtcg                                      <  1:393054/140‑1 (MQ=255)
                                                                                           gaGCTCCATCAGCAAGGTGCTGTGATACGCGCACTATCTTAACAATCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATAtcgtcg                                      <  2:265133/140‑1 (MQ=255)
                                                                                              cTCCATCAGCAAGGCACTGGGATACTCGCTCTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCgg                                   <  1:309982/140‑1 (MQ=255)
                                                                                              cTCCATCAGCAAGGCACTGGGACACTCGCCCTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCgg                                   <  1:378875/140‑1 (MQ=255)
                                                                                               tCCATCAGCAAGGCGATGGGATTCGGGCATTATCTCATCAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGt                                  <  1:370681/140‑1 (MQ=255)
                                                                                               tCCATCAGAAAGGCGCTGGGATGCGCGCCCGATTTTAACGGTCCGTGGCTCCAGCTCTCCCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCGGGAATCGGGCATATCGCCGCGGt                                  <  2:403350/140‑1 (MQ=255)
                                                                                                 cATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTAc                                >  1:564253/1‑140 (MQ=255)
                                                                                                 cATCAGCAAGGCGCTGGGATACCCGCACTATGTCACCAGCCCGTGGTTGCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTAc                                <  2:295454/140‑1 (MQ=255)
                                                                                                      gCAAGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACAtgtg                           <  1:586395/140‑1 (MQ=255)
                                                                                                        aaGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACATGTGaa                         >  2:143706/1‑140 (MQ=255)
                                                                                                         aGGCGCTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACATGTGaaa                        <  1:68423/140‑1 (MQ=255)
                                                                                                            cgcTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACATGTGaaa                        >  1:487279/1‑137 (MQ=255)
                                                                                                            cgcTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACATGTGaaa                        <  2:487279/137‑1 (MQ=255)
                                                                                                              cTGGGATACGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACATGTGAAAGacac                   >  2:743864/1‑140 (MQ=255)
                                                                                                                     aCGCGCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACATGTTAAAGACACCAAACCt            >  2:54357/1‑140 (MQ=255)
                                                                                                                      cgcgCACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACATGTGAAAGACACCAAACCTg           <  2:564253/140‑1 (MQ=255)
                                                                                                                         gcACTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCAGATCGTCGAGGTTCAAGTGAAAGACACCAAACCTGGCg        >  1:289990/1‑140 (MQ=255)
                                                                                                                           aCTATCTCAACAATCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACATGTGAAAGACACCAAACCTGGCGTc      >  1:720431/1‑140 (MQ=255)
                                                                                                                            cTATCTCAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACATGTGAAAGACACCAAACCTGGCGTCt     >  1:621188/1‑140 (MQ=255)
                                                                                                                               tctcAACAGTCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTTGCGGTACATGTCAACGACACCTAACCTGGCGTCTTCa  >  1:376878/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                   
GTCGTTTCCGTGACGGCCTGAAAGAGAGCGTTGAGATGGCAAGCCGCGCGCAGGTGACGCTGGCGATGGAGATCATGGATTATCCGTTGATGAGCTCCATCAGCAAGGCGCTGGGATACGCGCACTATCTCAACAATCCGTGGTTCCAGCTCTACCCGGATATCGGCAACCTGTCGGCGTGGGACAACGATGTGCAGATGGAGTTGCAGGCCGGAATCGGGCATATCGTCGCGGTACATGTGAAAGACACCAAACCTGGCGTCTTCA  >  NC_000913/4423240‑4423506

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