Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,159,336 A→C 15.8% intergenic (+34/+26) ptsG → / ← fhuE fused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,3360AC15.8% 163.6 / 22.6 68intergenic (+34/+26)ptsG/fhuEfused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base C (3/8);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGT  >  NC_000913/1159238‑1159433
                                                                                                  |                                                                                                 
cGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                  >  1:938935/1‑100 (MQ=255)
     ttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                             <  1:1332728/100‑1 (MQ=255)
     ttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                             <  1:677754/100‑1 (MQ=255)
     ttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                             >  1:1601927/1‑100 (MQ=255)
       cGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGcc                                                                                           <  1:984680/100‑1 (MQ=255)
        ggTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCt                                                                                          >  1:900091/1‑100 (MQ=255)
          tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                        <  1:1415401/100‑2 (MQ=255)
           aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                       <  2:1544961/100‑3 (MQ=255)
           aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                       <  2:1537738/100‑3 (MQ=255)
           aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                       <  1:2037160/100‑3 (MQ=255)
           aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtttt                                                                                       >  1:851361/1‑100 (MQ=255)
           aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtttt                                                                                       >  1:1589277/1‑100 (MQ=255)
            cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                      <  2:1589265/100‑1 (MQ=255)
            cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                      >  1:1698292/1‑100 (MQ=255)
            cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                      <  2:562668/100‑1 (MQ=255)
             tAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttagtc                                                                                     <  2:1487958/100‑5 (MQ=39)
             tAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtttttt                                                                                     <  1:1240503/100‑1 (MQ=255)
             tAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtttttt                                                                                     <  1:112717/100‑1 (MQ=255)
              aaaTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGTCAGCCATCTGGCTGCCttagtct                                                                                    <  2:148341/100‑6 (MQ=255)
              aaaTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCATCCATCTGGCTGCCttagtct                                                                                    <  2:149225/100‑6 (MQ=255)
              aaaTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTa                                                                                    >  2:655366/1‑100 (MQ=255)
                aTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACa                                                                                  >  2:1491929/1‑100 (MQ=255)
                aTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACa                                                                                  <  1:365130/100‑1 (MQ=255)
                aTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACa                                                                                  <  1:2124854/100‑1 (MQ=255)
                aTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACa                                                                                  <  2:274288/100‑1 (MQ=255)
                 tCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTACTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACa                                                                                  <  1:1370463/99‑1 (MQ=255)
                    gATAACCTGAAAACCGAGATGGATGAGTATATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACGAAGGCAGCCAGATGGCTGCCTTTTTTACAGgtg                                                                              <  2:1429536/100‑1 (MQ=255)
                     aTAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgtgt                                                                             <  2:1782401/100‑1 (MQ=255)
                      tAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGtt                                                                            >  1:2076130/1‑100 (MQ=255)
                           tGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTCAGGCAGCCAGCTGGATGCCTTTTTTACAGGTGTTCTTCa                                                                       >  2:1977312/1‑100 (MQ=255)
                             aaaaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGACTGCCTTTTTTACAGGTGTTATTCAGa                                                                     <  2:87391/100‑1 (MQ=255)
                                aCCGAGATGGATGGGTACATCCGTAACCACCAATCTGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTGTTTTACAGGTGTTATTCAGAAtt                                                                  <  2:1753340/100‑1 (MQ=255)
                                 ccGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTg                                                                 <  2:1659827/100‑1 (MQ=255)
                                   gagaTGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGAt                                                               <  1:2228452/100‑1 (MQ=255)
                                   gagaTGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGAt                                                               <  2:792480/100‑1 (MQ=255)
                                       tggatgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGCTACGt                                                           >  2:1380365/1‑100 (MQ=255)
                                       tggatgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGt                                                           >  1:1964194/1‑100 (MQ=255)
                                          atgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGcc                                                        <  2:961071/100‑1 (MQ=255)
                                            gAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCgg                                                      <  1:601469/100‑1 (MQ=255)
                                                 cATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTTCAGGTGTTATTCAGAATTGATACGTGCCGGTAAt                                                  >  2:968111/1‑99 (MQ=255)
                                                   tCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCt                                               >  2:497144/1‑100 (MQ=255)
                                                      gTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGaa                                            <  1:638464/100‑1 (MQ=255)
                                                       tAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGaaa                                           <  2:681619/100‑1 (MQ=255)
                                                           cACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAc                                       >  2:21948/1‑100 (MQ=255)
                                                               aaTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCgg                                   >  2:1037874/1‑100 (MQ=255)
                                                               aaTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCgg                                   >  1:1158200/1‑100 (MQ=255)
                                                                         aCGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAg                         <  1:1530393/100‑1 (MQ=255)
                                                                          cGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGa                        >  2:2142914/1‑100 (MQ=255)
                                                                          cGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGa                        >  2:2305584/1‑100 (MQ=255)
                                                                               gggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATa                   >  1:640894/1‑100 (MQ=255)
                                                                               gggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATa                   <  1:1891829/100‑1 (MQ=255)
                                                                                ggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAg                  >  1:135475/1‑100 (MQ=255)
                                                                                 gagaCTAAGGCAGCCAGATGGCTGGCTTTTTTACAGGTGTTAGTCAGAATTGATACGTGCCTGTAAGGCTGGAATTACGCGGTGTGCCGTAGACGATAGa                 >  2:704586/1‑100 (MQ=255)
                                                                                  agaCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGa                 <  1:2060976/99‑1 (MQ=255)
                                                                                   gaCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAc               <  1:700365/100‑1 (MQ=255)
                                                                                   aaaaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAc               >  1:2009047/5‑100 (MQ=39)
                                                                                    aaaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAcc              >  1:338925/4‑100 (MQ=255)
                                                                                     cTAAGGCAGCCAGATGGCTGCCTTTTTTATTTTTTTTTTTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCt             >  2:2173426/1‑100 (MQ=255)
                                                                                     cTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCt             >  1:1412318/1‑100 (MQ=255)
                                                                                     cTAAGGCAGCCAGATGGCGGCCTTTTTTACAGGTGTTATTCAGACTTGATACGTGACGCTAATGCTGAAATTCCACAGCGTCCCGTAGACGATAGAACCt             >  2:2125683/1‑100 (MQ=255)
                                                                                       aaGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTc           <  1:1595647/100‑1 (MQ=255)
                                                                                       aaGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTc           >  1:1582752/1‑100 (MQ=255)
                                                                                       aaGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTc           <  2:1698279/100‑1 (MQ=255)
                                                                                         ggaaGCCAGATGGCCGCCTTTTTTACAGGTGTTATCAAGAATTGATACGTACCGGTAATGCGGAAATTACCCGGGTTGCCGAAGAc                       >  2:1425629/4‑86 (MQ=255)
                                                                                          gCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCAc        >  2:344255/1‑100 (MQ=255)
                                                                                               cAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTgg   <  2:1615142/100‑1 (MQ=255)
                                                                                               cAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTgg   <  1:1297096/100‑1 (MQ=255)
                                                                                                aGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGt  >  2:508444/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
CGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGT  >  NC_000913/1159238‑1159433

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: