Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,755,639 G→A 100% intergenic (‑498/‑59) ydhZ ← / → pykF uncharacterized protein/pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,755,6390GA96.0% 226.5 / 0.4 74intergenic (‑498/‑59)ydhZ/pykFuncharacterized protein/pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (39/32);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAA  >  NC_000913/1755546‑1755735
                                                                                             |                                                                                                
gaCACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgc                                                                                            >  2:595260/1‑100 (MQ=255)
  cacaAGCACACATTCCTCTGCACGCTTTTTCGATTTCACCTATCCTGAGAGCGAGGCACCACCACTTTCTTAATACCGGATTCGCTTTCCGACAGTgcgc                                                                                          >  2:1246743/1‑100 (MQ=255)
  cacaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgcgc                                                                                          >  1:1238306/1‑100 (MQ=255)
  cacaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgcgc                                                                                          <  2:1548947/100‑1 (MQ=255)
     aaGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCa                                                                                       <  2:539711/100‑1 (MQ=255)
      aGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAg                                                                                      <  1:171293/100‑1 (MQ=255)
        cacacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGaa                                                                                    >  1:1978730/1‑100 (MQ=255)
          cacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAg                                                                                  <  1:78913/100‑1 (MQ=255)
                ccTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGtt                                                                            >  1:1721470/1‑100 (MQ=255)
                 ggtgcaaacaattttggtctttttcatgacagtcttagtctttaagttgagaaggatgggagaaacttGCTTTCCGACAGTGCGCCCAGAAAGCAAGttt                                                                           >  2:223725‑M1/69‑100 (MQ=255)
                  tctGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCTCCACTTTCGTAATACCGGATTCGCTTTCCGACAGTTCGCCCAGAAAGCAAGTTtc                                                                          >  1:100117/1‑100 (MQ=255)
                  tctGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTtc                                                                          <  2:1977434/100‑1 (MQ=255)
                     gCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       <  2:2233617/100‑1 (MQ=255)
                      cACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCTCCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTccc                                                                      <  2:1083855/100‑1 (MQ=255)
                       aacaattttggtctttttcatgacagtcttagtctttaagttgagaaggatgggagaaacttGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       >  1:406362‑M1/63‑98 (MQ=255)
                       aCGCTTTTTCGATGTCACCTATCCTTATAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       <  2:1959972/98‑1 (MQ=255)
                         caattttggtctttttcatgacagtcttagtctttaagttgagaaggatgggagaaacttGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATc                                                                   >  1:1347733‑M1/61‑100 (MQ=255)
                          cTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATcc                                                                  >  2:563913/1‑100 (MQ=255)
                             tttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtc                                                               <  1:1653563/100‑1 (MQ=255)
                             tgtCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtc                                                               <  2:1359326/98‑1 (MQ=255)
                               tCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtctc                                                             <  2:1721459/100‑1 (MQ=255)
                               tCGATGCCACCTATACATAAAGCGAGGCACCACCACTTTCGTAATACCGGATTCCCTTTCCGACAGTGCACCCAGAAAGCAAGTTTCTCCCATCCtt                                                                >  2:47404/1‑97 (MQ=255)
                                cGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCa                                                            >  2:1288926/1‑100 (MQ=255)
                                cGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCCGATTCGCTTTCCGACAGTGCGCTCAGAAAGCAAGTTTCTCCCATCCTtctt                                                             <  2:325390/99‑2 (MQ=255)
                                 gATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCaa                                                           >  1:488181/1‑100 (MQ=255)
                                  aTGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAAc                                                          <  2:1569290/100‑1 (MQ=255)
                                   tGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGTATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACt                                                         <  1:464630/100‑1 (MQ=255)
                                   tGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACt                                                         >  2:881434/1‑100 (MQ=255)
                                     tCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTa                                                       >  1:2156508/1‑100 (MQ=255)
                                     tCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTa                                                       >  1:1092007/1‑100 (MQ=255)
                                     tCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTa                                                       >  1:1216109/1‑100 (MQ=255)
                                      cACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaa                                                      >  2:2077039/1‑100 (MQ=255)
                                      cACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCATAAAGCAAGTTTCTCTCATCCTTCTCAACTTaa                                                      <  2:385545/100‑1 (MQ=255)
                                      cACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaa                                                      >  2:792939/1‑100 (MQ=255)
                                        ccTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaag                                                    <  1:1969162/100‑1 (MQ=255)
                                        ccTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaag                                                    >  2:1281203/1‑100 (MQ=255)
                                         cTATCCTTAGAGCGAGGTACCATCTCTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaaga                                                   <  2:2052236/100‑1 (MQ=255)
                                          tATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagac                                                  >  2:1963259/1‑100 (MQ=255)
                                            tCCTTAGAGCGAGGCACCACCCCTTTCTTAATACCGGATTCCCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagacta                                                <  2:1746812/100‑1 (MQ=255)
                                              cTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaag                                              >  2:1939754/1‑100 (MQ=255)
                                              cTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaag                                              >  1:1180820/1‑100 (MQ=255)
                                               ttAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaaga                                             <  1:792939/100‑1 (MQ=255)
                                                  gagCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaagact                                           >  1:85408/1‑99 (MQ=255)
                                                  gagCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTg                                          <  1:2109933/100‑1 (MQ=255)
                                                      gAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCAt                                      >  2:126530/1‑100 (MQ=255)
                                                       aGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATg                                     >  1:2005269/1‑100 (MQ=255)
                                                        ggCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGa                                    >  1:2190039/1‑100 (MQ=255)
                                                           accaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaa                                 <  1:648628/100‑1 (MQ=255)
                                                             caccacTTTCGTAATACCGGATTCGCTTTCGGAAAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGGCATGAAAAAg                               >  2:1475501/1‑100 (MQ=255)
                                                             caccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAg                               <  1:414096/100‑1 (MQ=255)
                                                             caccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAg                               >  1:1315328/1‑100 (MQ=255)
                                                             caccacTTTCGTAATACCGGAATCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATACTTCTCAACTTAAAGACTAAGACTGTCATGaa                                   >  2:1327875/1‑96 (MQ=255)
                                                              accacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGa                              >  2:244426/1‑100 (MQ=255)
                                                              accacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGa                              >  1:1790976/1‑100 (MQ=255)
                                                                   tttCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCaaa                         <  1:2077053/100‑1 (MQ=255)
                                                                     tCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTAAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAAt                       >  2:1190987/1‑100 (MQ=255)
                                                                      cGTAATTCCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTACCACTGTCATGAAAAAGACCAAAAtt                      <  2:168628/100‑1 (MQ=255)
                                                                       ntAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTg                     >  1:843663/2‑100 (MQ=255)
                                                                        tAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGt                    <  1:1708987/100‑1 (MQ=255)
                                                                           tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTg                 <  2:1525744/100‑1 (MQ=255)
                                                                            aCCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGc                >  2:1244208/1‑100 (MQ=255)
                                                                             ccGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGc                >  1:1904526/1‑99 (MQ=255)
                                                                             ccGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCa               <  1:1674880/100‑1 (MQ=255)
                                                                             ccGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCa               <  2:1040169/100‑1 (MQ=255)
                                                                              cGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCAc              <  2:1230425/100‑1 (MQ=255)
                                                                              cGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCAc              <  1:1356915/100‑1 (MQ=255)
                                                                               ggATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATTAAAAAGACCAAAATTGTTTGCAcc             >  2:305706/1‑100 (MQ=255)
                                                                                  ttCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATAGTTTGCACCATc          >  1:1706558/1‑100 (MQ=255)
                                                                                      CTTTCCGACCGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACGGTCATAAAAAAGACCAAAATTGTTGGCACCATCGGAc      <  2:1514238‑M1/100‑1 (MQ=255)
                                                                                      CTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAc      <  2:994605‑M1/100‑1 (MQ=255)
                                                                                       TTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAcc     >  2:1843449‑M1/1‑100 (MQ=255)
                                                                                       TTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAcc     >  2:1185914‑M1/1‑100 (MQ=255)
                                                                                        TTCCGACAGTGCGCCCAAAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGa       <  2:432109‑M1/97‑1 (MQ=255)
                                                                                          CCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAACAGACCAAAATTGTTTGCACCATCGGACAGaa  <  2:1027353‑M1/100‑1 (MQ=255)
                                                                                             |                                                                                                
GACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAA  >  NC_000913/1755546‑1755735

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