Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,930,677 +G intergenic (+175/‑264) kup → / → rbsD low affinity potassium transport system/RbsD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,930,6761.G94.6% 233.6 / 8.2 92intergenic (+174/‑265)kup/rbsDlow affinity potassium transport system/RbsD
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base G (34/53);  total (38/54)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCTGCGGTCTGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTG  >  NC_000913/3930601‑3930745
                                                                             |                                                                      
ccGCTGCGGTCTGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTc                                                                         <  2:4063774/76‑1 (MQ=255)
  gCTGCGGTCTGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGc                                                                     >  2:2326022/1‑76 (MQ=255)
  gCTGCGGTCTGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGc                                                                     >  2:179813/1‑76 (MQ=255)
   cTGCGGTCTGGAAGGATTAAGTTGCCGGATGATGGAAACCGCCTTCTTTATGTCGCATGAGTCGTTGATCCTCggc                                                                     >  2:4259801/1‑74 (MQ=255)
    tGCGGTCTGGAAGGGTTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCggca                                                                    >  1:525498/1‑73 (MQ=255)
    tGCGGTCTGGAAGGATTAAGTTGCCGGATGATGGAAACCCCCTTCTTTATGTCGCATGAGTCGTTGATCCCCGcca                                                                   >  2:310099/1‑74 (MQ=255)
       ggTCTGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAAc                                                                 <  1:637156/76‑1 (MQ=255)
       ggTCTGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAAc                                                                 <  2:2320550/76‑1 (MQ=255)
       ggTCTGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAAc                                                                 <  1:4387101/76‑1 (MQ=255)
         tCTGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGc                                                               >  1:34606/1‑76 (MQ=255)
          cTGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGcc                                                              <  2:4067982/76‑1 (MQ=255)
           tGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGccc                                                             <  2:4514749/76‑1 (MQ=255)
           tGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGccc                                                             >  1:1200473/1‑76 (MQ=255)
           tGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGccc                                                             >  2:3635199/1‑76 (MQ=255)
           tGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGccc                                                             <  1:3604817/76‑1 (MQ=255)
            ggAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCg                                                            <  2:648494/76‑1 (MQ=255)
             gAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGt                                                           <  1:3855528/76‑1 (MQ=255)
               aGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGtag                                                         >  1:1081413/1‑74 (MQ=255)
                ggATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGt                                                        >  1:4829369/1‑76 (MQ=255)
                  aTTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTAt                                                      >  1:3074738/1‑76 (MQ=255)
                  aTTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTAt                                                      <  1:4335102/76‑1 (MQ=255)
                   ttAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTAtt                                                     <  1:997054/76‑1 (MQ=255)
                    tAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTAttt                                                    >  1:1386778/1‑76 (MQ=255)
                     tctttGCCGGCTGCTGGACACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAGCGCCCGTGGTATTTg                                                   <  2:3082002/73‑1 (MQ=255)
                     aaGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTg                                                   <  2:349983/76‑1 (MQ=255)
                      aGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGc                                                  <  2:1900899/76‑1 (MQ=255)
                      aGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGc                                                  <  1:2326791/76‑1 (MQ=255)
                       gTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCg                                                 >  1:3104054/1‑76 (MQ=255)
                        ttGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGt                                                >  1:3667207/1‑76 (MQ=255)
                        ttGCCGGATGATGGAAACCTCCTTCTTTATGCCGCATGAGGCGGTGATCCTCCGCAAACCCCCGTGGTACTTGCgc                                                >  2:3814321/1‑75 (MQ=255)
                          gCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCt                                              <  2:2749568/76‑1 (MQ=255)
                           ccGGATGATGGAAACCTCCTTCTTTGTGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTg                                             >  1:1257251/1‑76 (MQ=255)
                           ccGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTg                                             >  1:578762/1‑76 (MQ=255)
                            cGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTgc                                            <  1:1749324/76‑1 (MQ=255)
                            cGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTgc                                            >  2:1841972/1‑76 (MQ=255)
                            cGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTgc                                            <  2:3006525/76‑1 (MQ=255)
                            cGGATGATGGAAACCTCCGTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTgc                                            <  2:1393879/76‑1 (MQ=255)
                             ggATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTgcg                                           <  1:3732304/76‑1 (MQ=255)
                              gatgatGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTgcgc                                          <  2:2482302/76‑1 (MQ=255)
                               gtgatgaGAACCTCCTTCTTTATTTGGCGTGAGTCGTTGATTGTCGGCAAACGCCCGTGGTATTTGGGTCTGcgcg                                         <  2:3437428/75‑1 (MQ=255)
                               atgatgGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGcgcg                                         <  2:4253151/76‑1 (MQ=255)
                               atgatgGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGcgcg                                         <  2:1994769/76‑1 (MQ=255)
                                   tgtaaaCCTCCTTCTGTGTGTCGCAGGAGTCGTGGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCaa                                     <  2:962977/73‑1 (MQ=255)
                                   tgGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCaa                                     <  1:768817/76‑1 (MQ=255)
                                    ggAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAg                                    >  1:3643724/1‑76 (MQ=255)
                                    ggAAACCTCCTTCTTTATGTCGCATGAGGCGTTGATCCTCGGCAAACGCCCGTGGCATTTGCGTCTGCGCGCCAAg                                    >  2:761220/1‑76 (MQ=255)
                                      aaaCCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCt                                  >  2:2758525/1‑76 (MQ=255)
                                       accctcctTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCtt                                 <  1:603142/74‑2 (MQ=255)
                                        aCCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGt                                >  1:3508980/1‑76 (MQ=255)
                                          ctcctTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTAc                              >  2:905267/1‑76 (MQ=255)
                                          ctcctTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTAc                              <  2:1951180/76‑1 (MQ=255)
                                          ctcctTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTAc                              >  2:3712719/1‑76 (MQ=255)
                                            cctTCTTTATGTCGCATGAGTCTTTGATCCTCGTCAATCGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACtt                            <  1:131845/76‑1 (MQ=255)
                                            cctTCTTTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACtt                            <  1:1250619/76‑1 (MQ=255)
                                             cttcttTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtg                           <  2:1644615/76‑1 (MQ=255)
                                             cttcttTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtg                           <  1:378579/76‑1 (MQ=255)
                                              ttcttTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCGGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgc                          <  1:3157505/76‑1 (MQ=255)
                                              ttcttTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgc                          <  2:2535923/76‑1 (MQ=255)
                                               tcttTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgct                         >  1:4803035/1‑76 (MQ=255)
                                               tcttTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgct                         >  1:1684428/1‑76 (MQ=255)
                                                cttTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgctg                        >  1:2605569/1‑76 (MQ=255)
                                                cttTATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgctg                        >  2:1731258/1‑76 (MQ=255)
                                                 tttATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgctgc                       >  1:4357761/1‑76 (MQ=255)
                                                  ttATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgctgct                      <  1:1379108/76‑1 (MQ=255)
                                                   tATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgctgctg                     <  1:677402/76‑1 (MQ=255)
                                                   tATGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgctgctg                     >  2:1769568/1‑76 (MQ=255)
                                                    aTGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgctgctgc                    >  2:2660900/1‑76 (MQ=255)
                                                    aTGTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTtgctgctgc                    >  2:613666/1‑76 (MQ=255)
                                                    aTGGCGCATGAGTCGTTGATCCTCGGCAATCGCCCGTGGTATTTGCGTCTGCGCGGCAATCTGTACTtgctgctgc                    <  2:3382428/76‑1 (MQ=255)
                                                      gTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCaa                  <  2:3166723/76‑1 (MQ=255)
                                                      gTCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCaa                  <  1:412867/76‑1 (MQ=255)
                                                       tCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAAc                 >  2:202372/1‑76 (MQ=255)
                                                       tCGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAAc                 <  1:1658053/76‑1 (MQ=255)
                                                        cGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCCACg                >  2:3920629/1‑76 (MQ=255)
                                                        cGCATGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACg                <  2:1847283/76‑1 (MQ=255)
                                                           aTGAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTaa             <  2:2821828/76‑1 (MQ=255)
                                                            tGAGGCGTTGATCCGCGGCGAAGGCCCGTGGTATTTGCGTCTGCGCGGCAAGGTGTACTTGCTGCTGCAACGTAAt            <  2:4857522/76‑1 (MQ=255)
                                                             gAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATg           <  2:2867383/76‑1 (MQ=255)
                                                             gAGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATg           <  1:1067720/76‑1 (MQ=255)
                                                              aGTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATgc          >  1:3958494/1‑76 (MQ=255)
                                                               ggcGTTGATCCTCGGCAAGCGCCCGTGGTGTTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACTTAATgcg         <  2:2851837/74‑1 (MQ=255)
                                                               gTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATgcg         <  2:2811896/76‑1 (MQ=255)
                                                               gTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATgcg         >  1:2225626/1‑76 (MQ=255)
                                                               gTCGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATgcg         <  1:4420773/76‑1 (MQ=255)
                                                               gTCGTGGATCCACGGCACACGCCCGTGGTCTTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATgcg         <  2:3183119/76‑1 (MQ=255)
                                                                 cGTTGATCCTCGGCAATCGGCCGTGGTATTTGCGACTGCGTGGCAAGCTGTACTTGCTGCTGCGACGTAATGCGCt       <  2:4402820/76‑1 (MQ=255)
                                                                 cGTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCt       <  1:606053/76‑1 (MQ=255)
                                                                  gTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTg      >  1:210971/1‑76 (MQ=255)
                                                                  gTTGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTg      >  1:3888964/1‑76 (MQ=255)
                                                                   ttGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGc     <  1:2517725/76‑1 (MQ=255)
                                                                   ttGATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGc     <  2:1734164/76‑1 (MQ=255)
                                                                     gATCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCg  >  2:4652471/1‑75 (MQ=255)
                                                                      aTCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTg  <  1:1731259/76‑1 (MQ=255)
                                                                      aTCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTg  <  2:3580419/76‑1 (MQ=255)
                                                                      aTCCTCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTg  >  1:4196014/1‑76 (MQ=255)
                                                                             |                                                                      
CCGCTGCGGTCTGGAAGGATTAAGTTGCCGGATGATGGAAACCTCCTTCTTTATGTCGCATGAGTCGTTGATCCTCGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTG  >  NC_000913/3930601‑3930745

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