Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,142,512 C→A 10.1% intergenic (‑291/‑18) ptsA ← / → frwC putative PTS multiphosphoryl transfer protein PtsA/putative PTS enzyme IIC component FrwC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,142,5120CA10.1% 111.5 / 2.3 53intergenic (‑291/‑18)ptsA/frwCputative PTS multiphosphoryl transfer protein PtsA/putative PTS enzyme IIC component FrwC
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base A (4/2);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATG  >  NC_000913/4142370‑4142650
                                                                                                                                              |                                                                                                                                          
ccGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAACATATGGACATGGGCTTTCGc                                                                                                                                            >  1:12270/1‑143 (MQ=255)
        tCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCGACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATGAAAAATATGGACATGGGTCTTCGCAAGCCCgg                                                                                                                                    >  2:768190/1‑143 (MQ=255)
          gCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCTAGCCCGGGa                                                                                                                                  >  1:833197/1‑143 (MQ=255)
          gCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGa                                                                                                                                  >  1:259439/1‑143 (MQ=255)
           cAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTGGCAAGCCCGGGag                                                                                                                                 >  1:344115/1‑143 (MQ=255)
           cAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGaa                                                                                                                                 >  1:838677/1‑142 (MQ=255)
             aGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTGCGAAAGcccatgtcca                                                                                                                               <  2:302274/143‑8 (MQ=255)
              gTTATGAGGGGGATCGCATTGTTGTACTGATATTACAGAAGGCCAGTAAATGGCCTTTTTATCCTCTGTTTGCTCTGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGTCATGGGCTTTCGCAAGCCCGGGagcag                                                                                                                              <  2:453386/143‑1 (MQ=255)
               ttATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCATTTTTATCCACTGTTTGCTCAGCCTGCGATTGACTATTGTTTCGTCCATCAAAAATATGGACATGGGCTTTCTAAAGCCCGGGaggagc                                                                                                                             >  2:568095/1‑143 (MQ=255)
                   gAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGAGTTTTTATTCACCAAATACTGAATAATTGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGCGCAGCTTat                                                                                                                         >  2:891094/1‑143 (MQ=255)
                      gCGGATCGCATTTTTCTACTGCTATTACAAAAGTCCCGTAAATGGCCTTTCTATCCACCGTTCGCTACGTCTGTTTGTGCCTACTGACGCCTTCAGCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTatgaa                                                                                                                      <  2:647380/143‑1 (MQ=255)
                             gCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCtt                                                                                                                           <  2:88202/131‑1 (MQ=255)
                             gCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCtt                                                                                                                           >  1:88202/1‑131 (MQ=255)
                             gCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTg                                                                                                               >  2:285087/1‑143 (MQ=255)
                             gCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTGGCTCCGACTGCGATTGCCTAGTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTg                                                                                                               <  1:692680/143‑1 (MQ=255)
                                tttttGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTg                                                                                                            <  1:472642/143‑1 (MQ=255)
                                tttttGTACTGATATCACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTg                                                                                                            <  1:901018/143‑1 (MQ=255)
                                         tGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTg                                                                                                   <  1:845812/143‑1 (MQ=255)
                                          gATATTACATGAGTCCGTGTAATTGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTGCATCAAGAATATGGACATGGGCGTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGa                                                                                                  <  2:458731/143‑1 (MQ=255)
                                              ttACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGaaaaa                                                                                              >  1:434380/1‑143 (MQ=255)
                                                  aaaaGTCCAGTAAATGGCCATTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCGTATGAAAGAGTTGGTGCAGATCCTGAAAAATAcc                                                                                          >  2:715262/1‑143 (MQ=255)
                                                   aaaGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCGCCGCCTGCGATTGCCTACTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATAccc                                                                                         <  1:764357/143‑1 (MQ=255)
                                                                tGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCACTAATATGGACATGGGCTTTCGCAAGCCCGGTAGCAGCTTATGAATGAGTTGGTGCAGATCCTGCAAAATACCCGTCAGCATTTact                                                                            >  2:584193/1‑141 (MQ=255)
                                                                    cTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACg                                                                        >  2:369037/1‑143 (MQ=255)
                                                                     tttttATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGCGCAGCTTATGAATGAGTTGGTGCAGATCCTGAGAAATACTAGTCAGCATGTAAtaccgg                                                                       >  2:454905/1‑138 (MQ=255)
                                                                          aTCCACTGTTTTCTCCTCCTGTGATTGCCTTGTTTTTCCTCCATCAAAAATATGGACATGGGCTTTCGAAAGCCCGGGAGCAGCTTATGAATGAGTCGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGtt                                                                  <  1:41191/143‑1 (MQ=255)
                                                                           tCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGttt                                                                 >  1:545643/1‑143 (MQ=255)
                                                                              aCTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCGATCAAAAATATGGACATGGGCTTTCGAAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATACTGAAAAATAGACATCAGGAGTTAATGACGGGCGtataaa                                                              >  2:283299/1‑138 (MQ=255)
                                                                              aCTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCACGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTcac                                                              >  2:149816/1‑143 (MQ=255)
                                                                               cTGTTGGCTCCGCCTGCCGTTGTCTATTCTTATCTGCATCAAAAATATGGACATGCGCTTTAGCAACCCCGGTAACAGCTTATGAATGAGTTGGTAAAGAAAAAGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTcaca                                                             <  2:711745/143‑1 (MQ=255)
                                                                                      cTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATCAATGAGTTGGTGCATATCCTGACAAATACCCGACAGCACTTAACGACGCGCCGTTcacacacaatt                                                      >  2:118771/1‑138 (MQ=255)
                                                                                        ccGCCTGCGCTTGCCTATTGCTTCCTCCAGCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATTAATGAGTTGGTGCAGATCCTGAAAACTACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTcc                                                    <  2:745623/143‑1 (MQ=255)
                                                                                              gCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGt                                              <  2:128002/143‑1 (MQ=255)
                                                                                              gCGATTGCCTATTGTTGCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGt                                              <  2:557585/143‑1 (MQ=255)
                                                                                               cGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTg                                             <  1:632039/143‑1 (MQ=255)
                                                                                                       tctggtctcctccATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGTTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCggg                                     <  2:716615/136‑1 (MQ=255)
                                                                                                       tATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAACACACAACACCTGGGTGTGATGACGGGCGTTTCACACATGATTCCCTTCGTGGGATCggg                                     >  1:205596/1‑143 (MQ=255)
                                                                                                         ttGTTTCCGCCCTCAAGTATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCg                                   <  2:429328/143‑1 (MQ=255)
                                                                                                          tGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCgg                                  >  2:901702/1‑143 (MQ=255)
                                                                                                            tttCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTa                                >  1:843342/1‑143 (MQ=255)
                                                                                                                ctccATCACAACTATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTAtttt                            <  2:670434/143‑1 (MQ=255)
                                                                                                                ctccATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTAtttt                            >  1:646087/1‑143 (MQ=255)
                                                                                                                ctccATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAAAACCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTAtttt                            >  1:746616/1‑143 (MQ=255)
                                                                                                                   cATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCt                         >  1:383592/1‑143 (MQ=255)
                                                                                                                    aTCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGATTATTTTCCGTCAGCATTTAATGACGGGCGTTTCACACAGGATGGCCTTCGTGGTATCGGGCGGTATTTTGCTg                        <  1:169858/143‑1 (MQ=255)
                                                                                                                          actatGGACATGGGCTTTCGCGAGCCCGTGTGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGtt                  <  1:285087/141‑1 (MQ=255)
                                                                                                                             atGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTcc               >  1:12748/1‑143 (MQ=255)
                                                                                                                                  cccgggCTTGCGAAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTACGGCAt          >  1:873795/4‑143 (MQ=255)
                                                                                                                                    cgggCTTGCGAAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGGCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTACTTTGCTGGCGGTTCCGGTCgtgt        >  2:503760/2‑143 (MQ=255)
                                                                                                                                        cTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTa    >  2:842038/1‑143 (MQ=255)
                                                                                                                                        cTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTACCGGGCGGTAATTTGCTGGCGGGTTCCGTCATGTTGTa    >  2:42334/1‑143 (MQ=255)
                                                                                                                                         tttCGCAAGCCGGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGTTTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTAt   <  2:365114/143‑1 (MQ=255)
                                                                                                                                         tttCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGCTCCTGAAAAATACGCGTCAGCATTTAATGACGCGCGTTTCACACATGCTGCCCCTCGTGGGATGGGGGGGGAATTAGCTGACGGTATCCGACAcgttgcag   >  2:684840/1‑140 (MQ=255)
                                                                                                                                          ttCGCAAGCACGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATg  <  1:21131/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                          
CCGGTTTATCGCAAGTTATGAGGCGGATCGCATTTTTGTACTGATATTACAAAAGTCCAGTAAATGGCCTTTTTATCCACTGTTTGCTCCGCCTGCGATTGCCTATTGTTTCCTCCATCAAAAATATGGACATGGGCTTTCGCAAGCCCGGGAGCAGCTTATGAATGAGTTGGTGCAGATCCTGAAAAATACCCGTCAGCATTTAATGACGGGCGTTTCACACATGATTCCCTTCGTGGTATCGGGCGGTATTTTGCTGGCGGTTTCCGTCATGTTGTATG  >  NC_000913/4142370‑4142650

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