Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,348,184 T→C 8.1% intergenic (+91/‑123) npr → / → yrbL phosphohistidinoprotein‑hexose phosphotransferase component of N‑regulated PTS system (Npr)/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,348,1840TC 8.1% 313.3 / 16.7 111intergenic (+91/‑123)npr/yrbLphosphohistidinoprotein‑hexose phosphotransferase component of N‑regulated PTS system (Npr)/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (48/54);  new base C (5/4);  total (53/58)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 118 total aligned reads displayed.

CTCCCTTCAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTT  >  W3110S.gb/3348116‑3348254
                                                                    |                                                                      
cTCCCTTCAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTa                                                                       >  1:514320/1‑70 (MQ=255)
 tCCCTTCAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTAc                                                                      >  1:279411/1‑70 (MQ=255)
    cTTCAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACAttt                                                                  <  1:2379562/71‑1 (MQ=255)
    cTTCAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACAttt                                                                  <  1:3115446/71‑1 (MQ=255)
    cTTCAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACAttt                                                                  <  1:763840/71‑1 (MQ=255)
       cAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTtgt                                                               >  1:1233321/1‑71 (MQ=255)
       cAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTtgt                                                               >  1:2905943/1‑71 (MQ=255)
       cAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTtg                                                                <  1:4349657/70‑1 (MQ=255)
            tCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCt                                                          >  1:4603641/1‑71 (MQ=255)
              cccccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCttt                                                        >  1:5121720/1‑71 (MQ=255)
              cccccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCttt                                                        >  1:473304/1‑71 (MQ=255)
              cccccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCtt                                                         <  1:4437932/70‑1 (MQ=255)
              cccccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCtt                                                         <  1:125468/70‑1 (MQ=255)
              cccccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCtt                                                         <  1:4543032/70‑1 (MQ=255)
              cccccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCt                                                          <  1:2522778/69‑1 (MQ=255)
              cccccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGCACATTTTGTGCGCttt                                                        <  1:2865362/71‑1 (MQ=255)
               ccccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCtttt                                                       <  1:2828969/71‑1 (MQ=255)
                cccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCtttt                                                       >  1:4540041/1‑70 (MQ=255)
                cccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCtttt                                                       >  1:3097930/1‑70 (MQ=255)
                 ccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCtttt                                                       >  1:4181710/1‑69 (MQ=255)
                 ccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCtttt                                                       >  1:3947839/1‑69 (MQ=255)
                 ccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGCACATTTTGTGCGCTTTTCa                                                     >  1:603999/1‑71 (MQ=255)
                 ccGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGCACATTTTGTGCGCTTTTCa                                                     <  1:249187/71‑1 (MQ=255)
                   gACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCAc                                                    >  1:134129/1‑70 (MQ=255)
                   gACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACt                                                   >  1:1312030/1‑71 (MQ=255)
                   gACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACt                                                   <  1:4180748/71‑1 (MQ=255)
                        tAAGAAAATCACGTACTGAAATCTTTCTCAATCAACGTCATTTGt                                                                        <  1:2008052/45‑1 (MQ=255)
                        tAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGt                                                                        <  1:2967655/45‑1 (MQ=255)
                        tAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGt                                                                        <  1:1554487/45‑1 (MQ=255)
                        tAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGt                                                                        <  1:153672/45‑1 (MQ=255)
                        tAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGt                                                                        <  1:4423934/45‑1 (MQ=255)
                        tAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTATGCGCtt                                                         >  1:2496649/1‑60 (MQ=255)
                        tAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGCACATTTTGTGCGCTTTTCACttt                                                 >  1:1640631/1‑68 (MQ=255)
                        tAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGCACATTTTGTGCGCTTTTCACttt                                                 >  1:4805707/1‑68 (MQ=255)
                          aGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTa                                                                       <  1:5100297/44‑1 (MQ=255)
                          aGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTa                                                                       <  1:3889617/44‑1 (MQ=255)
                          aGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTa                                                                       <  1:1011099/44‑1 (MQ=255)
                          aGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTa                                                                       <  1:1153138/44‑1 (MQ=255)
                          aGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCaga                                             <  1:688290/70‑1 (MQ=255)
                          aGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCaga                                             <  1:208398/70‑1 (MQ=255)
                          aGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCag                                              >  1:3081188/1‑69 (MQ=255)
                          aGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCag                                              >  1:399296/1‑69 (MQ=255)
                          aGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCag                                              >  1:4188113/1‑69 (MQ=255)
                           gAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTc                                                <  1:2235686/66‑1 (MQ=255)
                              aaTCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAAc                                        >  1:2054919/1‑71 (MQ=255)
                                  aCGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCtt                                                         >  1:518315/1‑50 (MQ=255)
                                  aCGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCtt                                     >  1:3592566/1‑70 (MQ=255)
                                       cTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGaat                               <  1:1476831/71‑2 (MQ=255)
                                       cTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGaa                                <  1:4609752/70‑1 (MQ=255)
                                       cTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGaa                                >  1:3402392/1‑70 (MQ=255)
                                       cTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGaa                                >  1:733337/1‑70 (MQ=255)
                                         gAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCa                                                     >  1:2537461/1‑47 (MQ=255)
                                         gAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCag                                              >  1:4142771/1‑54 (MQ=255)
                                         gAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAAc                             >  1:878796/1‑71 (MQ=255)
                                          aaaTCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCtttt                                                       >  1:4879412/1‑44 (MQ=255)
                                          aaaTCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAAcc                            >  1:4332549/1‑71 (MQ=255)
                                           aaTCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCagaag                                           <  1:713047/55‑1 (MQ=255)
                                            aTCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCtt                          <  1:1830889/71‑1 (MQ=255)
                                            aTCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCtt                          <  1:641485/71‑1 (MQ=255)
                                            aTCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCtt                          <  1:1904567/71‑1 (MQ=255)
                                             tCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCtt                          <  1:619998/70‑1 (MQ=255)
                                             tCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCtt                          <  1:4493215/70‑1 (MQ=255)
                                             tCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCtt                          <  1:4291729/70‑1 (MQ=255)
                                               gTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTaa                        <  1:3513154/70‑1 (MQ=255)
                                               gTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTaa                        <  1:4046719/70‑1 (MQ=255)
                                               gTTCTCAATCAACGTCATTTGCACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTaa                        <  1:829735/70‑1 (MQ=255)
                                                ttCTCAATCAACGTCATTTCTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCtt                          <  1:4856024/67‑1 (MQ=255)
                                                 tctcAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCtt                          <  1:1009989/66‑1 (MQ=255)
                                                 tctcAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTaa                        <  1:3963119/68‑1 (MQ=255)
                                                  ctcAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAgg                    >  1:2196632/1‑71 (MQ=255)
                                                   tcaatcaaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTc                                                <  1:1885701/42‑1 (MQ=255)
                                                      atcaaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCAtt                >  1:3353483/1‑71 (MQ=255)
                                                       tcaaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTaa                        <  1:1281390/62‑1 (MQ=255)
                                                       tcaaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTg               >  1:1739381/1‑71 (MQ=255)
                                                         aaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTaa                                   >  1:1458081/1‑49 (MQ=255)
                                                         aaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTaa                        <  1:1662885/60‑1 (MQ=255)
                                                         aaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTaa                        >  1:2153139/1‑60 (MQ=255)
                                                         aaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGtt             <  1:389932/71‑1 (MQ=255)
                                                         aaCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGtt             <  1:150787/71‑1 (MQ=255)
                                                          aCGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGttt            <  1:505075/71‑1 (MQ=255)
                                                             tCATTTGTACATTTTGTGCGCTTTTCACTTTCAGATGAACCTTAAGAAAACCTTAAGAGGCATTGtt             <  1:4909449/67‑1 (MQ=255)
                                                             tCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAgg         <  1:1490944/71‑1 (MQ=255)
                                                             tCATTTGCACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAag                       >  1:2582262/1‑57 (MQ=255)
                                                             tCATTTGCACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAgg         >  1:1944119/1‑71 (MQ=255)
                                                               aTTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAagag                     >  1:893981/1‑57 (MQ=255)
                                                               aTTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGtt       >  1:2090194/1‑71 (MQ=255)
                                                               aTTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGtt       >  1:3250413/1‑71 (MQ=255)
                                                               aTTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGtt       >  1:3489954/1‑71 (MQ=255)
                                                               aTTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGtt       <  1:1524265/71‑1 (MQ=255)
                                                               aTTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGtt       >  1:4588128/1‑71 (MQ=255)
                                                                tttGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTaa                        <  1:2617930/53‑1 (MQ=255)
                                                                tttGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGttt      <  1:324549/71‑1 (MQ=255)
                                                                tttGCACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGttt      <  1:2676947/71‑1 (MQ=255)
                                                                  tGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTg    >  1:877768/1‑71 (MQ=255)
                                                                   gTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTaa                                   <  1:4202520/39‑1 (MQ=255)
                                                                   gTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTa           <  1:2079399/63‑1 (MQ=255)
                                                                   gTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTa           <  1:356732/63‑1 (MQ=255)
                                                                   gTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAgg         >  1:1559099/1‑65 (MQ=255)
                                                                    tACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGtt  >  1:483195/1‑71 (MQ=255)
                                                                    tACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGtt  <  1:2854466/71‑1 (MQ=255)
                                                                    |                                                                      
CTCCCTTCAAACTCCCCCCGACAATAAGAAAATCACGTACTGAAATCGTTCTCAATCAACGTCATTTGTACATTTTGTGCGCTTTTCACTTTCAGAAGAACCTTAAGAAAACCTTAAGAGGCATTGTTTAGGTTTTGTT  >  W3110S.gb/3348116‑3348254

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