Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,897,086 T→C C105R (TGT→CGT)  yieH → phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,897,0860TC95.7% 158.0 / ‑1.7 46C105R (TGT→CGT) yieHphosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (20/24);  total (21/25)
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.20e-01

CTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCTACGATATTCAGCGCTGGAAGCCAGACCCGGC  >  NC_000913/3896951‑3897210
                                                                                                                                       |                                                                                                                            
cTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgcgtg                                                                                                                           >  2:343315/1‑139 (MQ=255)
cTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgcgtg                                                                                                                           <  1:554605/139‑1 (MQ=255)
cTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgcgtg                                                                                                                           <  2:512660/139‑1 (MQ=255)
           tgtTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAAc                                                                                                                >  1:89541/1‑139 (MQ=255)
                cGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGccc                                                                                                           <  2:253682/139‑1 (MQ=255)
                   tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCaaa                                                                                                        >  1:311751/1‑139 (MQ=255)
                    aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAAt                                                                                                       >  2:585013/1‑139 (MQ=255)
                     gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATa                                                                                                      >  2:946979/1‑139 (MQ=255)
                        aaaaCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACa                                                                                                   <  1:585013/139‑1 (MQ=255)
                        aaaaCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACa                                                                                                   <  2:596554/139‑1 (MQ=255)
                             aGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATg                                                                                              <  2:973695/139‑1 (MQ=255)
                             aGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATg                                                                                              <  1:946979/139‑1 (MQ=255)
                                 gCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                                                                          <  1:343315/139‑1 (MQ=255)
                                     aaCACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTc                                                                                      >  2:301211/1‑139 (MQ=255)
                                      acacGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                                                                          >  1:390569/1‑134 (MQ=255)
                                      acacGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                                                                          <  2:390569/134‑1 (MQ=255)
                                         cGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATg                                                                                  <  1:117630/139‑1 (MQ=255)
                                           tttACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgcgtgtg                                                                                                                         >  1:180622/1‑98 (MQ=255)
                                           tttACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgcgtgtg                                                                                                                         <  2:180622/98‑1 (MQ=255)
                                            ttACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                                                               <  1:1137779/139‑1 (MQ=255)
                                               ccGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAg                                                                            <  1:261418/139‑1 (MQ=255)
                                                  tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTg                                                                         >  1:518818/1‑139 (MQ=255)
                                                   gCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGc                                                                                                             <  1:1167880/102‑1 (MQ=255)
                                                   gCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGc                                                                                                             >  2:1167880/1‑102 (MQ=255)
                                                             cTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCa                                                              <  2:89541/139‑1 (MQ=255)
                                                                 gCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                                                               <  1:1133959/118‑1 (MQ=255)
                                                                 gCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                                                               <  2:930221/118‑1 (MQ=255)
                                                                 gCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                                                               >  2:1133959/1‑118 (MQ=255)
                                                                 gCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                                                               >  1:930221/1‑118 (MQ=255)
                                                                   tGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACtt                                                        <  2:400080/139‑1 (MQ=255)
                                                                   tGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACtt                                                        <  2:518818/139‑1 (MQ=255)
                                                                          ttCAGAACTGGAAGCCGTCGAAGGGGTTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGAt                                                 >  2:72758/1‑139 (MQ=255)
                                                                              gAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAAc                                             <  2:942401/139‑1 (MQ=255)
                                                                               aaCTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACt                                            >  2:205619/1‑139 (MQ=255)
                                                                               aaCTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACt                                            >  2:741832/1‑139 (MQ=255)
                                                                                    gAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCa                                       >  1:638204/1‑139 (MQ=255)
                                                                                        ccATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCgg                                   <  1:205619/139‑1 (MQ=255)
                                                                                        ccATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCgg                                   <  2:638204/139‑1 (MQ=255)
                                                                                         cATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCa                                       <  2:1092902/134‑1 (MQ=255)
                                                                                         cATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCa                                       >  1:1092902/1‑134 (MQ=255)
                                                                                                     tGGAGCGCTCCTGTCAGCGATCACTGGGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAGATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCTACGATATTCAg                      <  1:301211/139‑1 (MQ=255)
                                                                                                            cTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCTa                                >  2:55101/1‑122 (MQ=255)
                                                                                                            cTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCTa                                <  1:55101/122‑1 (MQ=255)
                                                                                                            cTCCTGTCAGCGATCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCTACGATATTCAGCGCTGGa               >  1:280577/1‑139 (MQ=255)
                                                                                                                         tCACTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCTACGATATTCAGCGCTGGAAGCCAGACCCGGc  >  1:1162560/1‑139 (MQ=255)
                                                                                                                           aCTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCTACGATATTCAGTGCTGGAAGCCAGACCCgg   >  1:469462/1‑136 (MQ=255)
                                                                                                                           aCTGCGCCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCTACGATATTCAGTGCTGGAAGCCAGACCCgg   <  2:469462/136‑1 (MQ=255)
                                                                                                                                gcCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATggg                                                                                >  2:802846/1‑54 (MQ=255)
                                                                                                                                gcCAATGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATggg                                                                                <  1:802846/54‑1 (MQ=255)
                                                                                                                                    aTGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCTACGATATTCAGCGCTGGAAGCCAGAcc      >  1:48617/1‑124 (MQ=255)
                                                                                                                                    aTGCGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCTACGATATTCAGCGCTGGAAGCCAGAcc      <  2:48617/124‑1 (MQ=255)
                                                                                                                                       |                                                                                                                            
CTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCTACGATATTCAGCGCTGGAAGCCAGACCCGGC  >  NC_000913/3896951‑3897210

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